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1.  SOLID  STATE 


1.1  SEMICONDUCTORS 

1.1.1  Electron  Probe  Analysis  of  Semiconductors 

AF-AFOSR-68- 1414,  Air  Force  Office  of  Scientific  Research 
D.  B.  Wittry,  D.  F.  Kyser,  J.  H.  McCoy,  H.  C.  Marciniak, 

T.  Rao-Sahib  and  A.  Van  Couvering 

Investigations'  ire  being  conducted  in  the  following  areas: 
cathodoluminescence,  electron  beam  modulated  reflectance,  surface 
state  densities,  non-linear  effects  in  cathodoluminescence  electron 
probe  x-ray  microanalysis  and  selected-area  electron  spectrometry.  A 
brief  report  of  the  present  status  and  recent  progress  in  these  investi¬ 
gations  follows: 

1.  1. 1.  1  Cathodoluminescence  of  Semiconductors 

D.  F.  Kyser,  H.  Marciniak  and  D.  B.  Wittry 

Work  on  the  ultrahigh  vacuum  chamber  and  beam  column  has 
continued  since  the  last  semiannual  report.  With  this  system,  we  hope 
to  provide  a  clean  "atmosphere"  in  which  to  study  semiconductor  speci¬ 
mens  with  minimum  effects  due  to  contamination.  The  vacuum  chamber, 
manufactured  by  Ion  Dynamics  Corporation  has  been  attached  to  an 

i 
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Ultek  TNB  ultrahigh  vacuum  system.  A  pressure  of  3  x  10  Torr  was 

achieved  after  only  a  few  hours  of  pumping  and  it  appears  possible  to 

.  9 

reach  the  goal  of  better  than  10  Torr. 

With  this  system,  it  will  be  possible  to  vary  the  specimen 
o  o 

temperature  from  20  K  to  500  K  using  an  Air  Products  and  Systems 
Cryo  Tip  refrigerator  and  an  electric  heater.  A  single  magnetic  lens 
will  be  used  to  focus  an  electron  beam  on  the  specimen.  A  standard 
ARL  microprobe  light  optics  assembly  will  be  used  to  focus  the  lumines¬ 
cence  on  a  high  resolution  grating  spectrometer.  The  first  application 
of  this  system  will  be  the  study  of  the  temperature  dependence  of  cathod- 
oluminescence  of  semiconductors.  Additional  experiments  planned 
include  electron  beam  modulated  reflectance  at  low  temperatures  and 
studies  of  the  changes  in  surface  recombination  velocities  with  surface 
conditions. 


1.1.  1.2  Electron  Beam  Modulated  Reflectance  of  Sf  liconductors 
J.  McCoy  and  D.  B.  Wittry 

As  reported  previously*,  a  focused  elecf^on  beam  has  been 
used  to  modulate  the  reflectivity  of  germanium.  Further  experiments 
have  been  made  since  the  last  semiannual  report  to  determine  the  rela¬ 
tive  importance  of  excess  carrier  density  and  thermal  effects.  For 
example,  the  amplitude  of  the  modulated  response  has  been  measured 
as  a  function  of  the  beam  current  density.  If  the  carrier  lifetime  is 
constant,  the  amplitude  of  the  modulated  reflectance  should  vary  as  the 
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one  third  power  of  the  beam  density  for  an  effect  caused  by  excess 
carriers.  On  the  other  hand,  if  thermal  effects  dominate,  the  modulated 
reflectance  should  vary  as  the  one  half  power  of  the  beam  density.  A 

one  third  power  dependence  was  found  over  four  orders  of  magnitude 

\ 

change  in  the  beam  density. 

The  heat  transport  equation  was  solved  for  the  case  of  a 

transient  disk-shaped  heat  source  on  the  surface  of  a  semi-infinite 
2 

medium.  This  model  approximates  the  case  of  heating  by  a  defocussed 
electron  beam  which  is  pulsed  as  in  the  studies  of  electron  beam  modulated 
reflectance.  The  time  response  was  obtained  from  a  digital  computation. 
These  calculations  have  not  ruled  out  the  possibility  of  thermal  effects 
as  a  cause  of  the  modulated  reflectance  since  the  response  time  was 
found  to  be  about  the  same  as  the  expected  carrier  lifetime  in  germanium. 
Further  experiments  are  in  progress  to  measure  the  actual  excess 
carrier  densities  obtained  during  electron  bombardment.  Hopefully 
these  experiments  will  give  more  insight  into  the  physical  basis  for  the 
electron  beam  modulation  of  reflectance. 
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1.  1,  1.3  Non-linear  Cathodoluminescence  of  Semiconductors 
T.  Rao- Sahib  and  D.  B.  Wittry 

The  voltage  dependence  of  cathodoluminescence  has  been 
previously  used*  to  measure  diffusion  lengths  and  surface  recombination 
velocities  of  n-type  gallium  arsenide.  However,  for  p-type  specimens, 
the  experimental  results  do  not  agree  with  theory,  or  appear  to  indicate 
unreasonably  low  surface  recombination  velocities  (e.g.  ,  S  =  0)  or 
unreasonably  large  space  charge  depletion  layer  widths  (i.  e.  ,  d  =  .  4u 
in  heavily  doped  material). 

Most  p-type  gallium  arsenide  has  exhibited  a  non-linear 

relation  between  luminescence  intensity  and  electron  beam  current  at 

constant  voltage.  Modifications  to  the  theory  for  a  one  dimensional  case 

can  be  made  to  account  for  this  non-linearity,  but  we  have  not  been 

able  to  construct  a  satisfactory  theory  foi  the  three  dimensional  case. 

For  this  reason,  we  have  attempted  to  use  experimental  conditions  that 

2 

more  ciosely  approximate  the  one  dimensional  case.  This  is  done  by 
irradiating  a  large  area  of  the  specimen  and  using  a  selected  area  aper¬ 
ture  in  the  plane  of  the  image  formed  by  the  objective  lens;  that  is  at 
the  entrance  slit  of  the  monochromator  Two  cases  have  been  examined, 
namely,  a)  a  beam  diameter  of  180  u  and  a  selected  area  corresponding 
to  27  by  27  u  ,  and  b)  a  beam  diameter  of  72  u  and  a  selected  area 

corresponding  to  5  ^  by  5  p .  The  experimental  results  obtained  on  a 

18  3 

Zn  doped  GaAs  specimen  (carrier  concentration  10  /cm  )  are  summar¬ 
ised  in  Table  I. 
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Table  I 


power  (mW) 

Average  exponent  m 

S 

d/L 

Mu) 

Case  a 

5 

1.0 

10 

0.033 

0.8 

50 

1.2 

10 

0.  033 

0.  75 

Case  b 

15 

1.2 

50 

0.033 

0.  96 

40 

1.4 

50 

0.033 

0.8 

It  can  be  seen  from  Table  I  that  the  values  obtained  with 
various  power  levels  are  consistent.  This  is  due  to  the  use  of  calculated 
curves  that  take  account  of  the  appropriate  exponent  m  in  the  relation 
I  =  im  where  I  is  the  intensity  and  i  is  the  current  at  constant  voltage. 

In  each  case,  the  average  value  of  the  exponent  m  is  determined  for  the 
power  level  actually  used  in  the  voltage  dependence  curves.  The  fit  of 
the  experimental  curves  to  the  calculated  curves  is  within  experimental 
error  for  the  cases  shown  in  Table  I,  indicating  some  progress  in  our 
understanding  of  the  voltage  dependence  of  cathodoluminescence  of  p-type 
GaAs.  Further  experiments  are  planned  on  p-type  specimens  of  GaAs 
with  various  doping  levels. 
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1 .  1 .  1 . 4  Surface  State  Density  Studies 

H.  Marciniak  and  D.  B.  Wittry 

Experiments  are  being  planned  to  measure  the  surface  state 
density,  width  of  the  space  charge  depletion  region,  and  surface  recom¬ 
bination  velocity  for  oxidized  Si  and  for  CdSe.  Preliminary  experiments 
are  being  conducted  in  the  ARL  microprobe,  but  because  of  the  lack  of  a 
controlled  environment,  sane  of  these  experiments  must  be  deferred 
until  the  new  electron  probe  column  is  finished. 

1.  1.  1.  5  Quantitative  Electron  Probe  Microanalysis 
D.  F.  Kyser  and  T.  Rao-Sahib 

Measurement  of  x-ray  intensity  ratios  from  pairs  of  semi¬ 
conductor  compounds  were  repeated  on  the  best  available  material  and 
the  results  compared  with  the  calculations  of  Brown  et  al.  *  The  ratio 

of  the  AsK  intensities  from  GaAs  and  InAs  was  measured  to  be  1.  173 

a 

(30  kV  beam  potential,  52.  5  degrees  takeoff  angle).  Transport  equation 

program  calculations  gave  1.209  without  a  fluorescence  correction  and 

1.  170  with  the  inclusion  of  a  fluorescence  correction.  For  the  case  of 

the  GaK  line  the  numbers  were  1. 828  (meas.  ),  1.  791  (calc.  ),  and  1. 843 
CL 

(calc,  with  fluor.  ).  For  the  case  of  the  InL  line  in  InP  and  InAs  the 

a 

numbers  were  1.  326  (meas.  )  and  1.  340  (calc.  ).  There  was  no  significant 
correction  for  fluorescence  in  this  third  case.  In  all  three  cases  the 
agreement  between  theory  and  experiment  was  1%  or  better. 
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1.  Brown,  D.  B.  ,  Wittry,  D.  B.  ,  and  Kyser,  D.F.  (to  be  published). 

1 .  1 .  1 . 6  Selec ted  A  rea  Electron  Spectrometry  in  the 
Transmission  Electron  Microscope 
D.  B.  Wittry 

An  electron  spectrometer  system  for  use  with  the  trans¬ 
mission  electron  microscope  was  constructed  at  the  Cavendish  Laboratory, 
University  of  Cambridge  (during  sabbatical  leave).  This  system  shown  in 
Fig,  1  can  be  used  for  the  study  of  the  energy  spectrum  of  electrons  from 
a  selected  area  of  the  image  or  diffraction  pattern  formed  in  normal 
operation  of  the  electron  microscope.  A  similar  system  is  now  being 
designed  £or  use  with  the  Hitachi  HU  12  5  at  the  University  of  Southern 
California, 

1  o 

Preliminary  results  obtained  on  a  1,000  A  selected  area 
of  i  polycrystalline  BeO  film  are  shown  in  Fig.  2.  In  this  figure,  the 
energy  loss  corresponding  to  oxygen  K  excitation  can  be  seen  as  well  as 
some  fine  structure  similar  to  that  observed  in  soft  x-ray  absorption 
spectra.  The  prelim  nary  results  indicate  that  the  x-ray  energy  loss 
spectra  can  be  used  for  microanalysis  in  the  transmission  electron 
microscope  with  sensitivity  significantly  higher  than  that  of  electron 
microprobe  x-ray  analysis  of  thin  specimens  and  light  elements. 

In  addition  to  practical  applications  to  microanalysis,  the 
x  ay  energy  loss  spectra  are  of  considerable  fundamental  interest. 
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Fundamental  investigations  of  the  fine  structure  in  the  energy  loss 
spectrum  can  be  performed  using  the  full  capabilitier  of  the  transmission 
electron  microscope  to  determine  the  crystal  structure  and  orientation  of 
the  specimen. 
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1.1.2  Site  Distribution  of  Silicon  in  Silicon-Doped  Gallium 
Arsenide 

AF-AFOSR -496-67,  Joint  Services  Electronics  Program 
Fi9628-6S-C-D169,  Air  Force  Cambridge  Research 
Labora  tories 

NGR-05-018-083,  National  Aeronautics  and  Space 
Administration 

W.  P.  Allred,  G,  Cumn  ing,  J.  Kung  and  W.  G.  Spitzer 


Introduction 

Gallium  arsenide  which  is  grown  from  a  stochiometric  melt 

and  doped  with  eilicon  is  always  an  n-type  semiconductor.  However,  it 

U  2)  18-3 

is  known  ’  that  at  silicon  concentrations  of  [Si]  >  10  cm  ,  the 

carrier  concentration  n,  in  the  exhaustion  range,  may  be  much  less  than 

[  Si  ]  .  The  difference  between  the  values  of  n  and  [Si  ]  is  related  to  the 

amphoteric  nature  of  silicon  in  GaAs.  Silicon  on  a  gallium  site,  Si^, 

is  a  donor  while  Si  is  expected  to  be  an  acceptor.  There  is  also  some 

A  8 
(3) 

evidence  from  photoluminescence  studies  for  the  existence  of  a  sub¬ 
stantial  concentration  of  (^Si^-  Si^  J  nearest  neighbor  pairs  when  the 

r,,  ,  ~  ..19  -3 

[  Si  ]  >  10  cm 

There  have  been  several  recent  studies  of  the  infrared 

absorption  associated  with  localized  vibrational  modes  of  impurities  in 

(4,  5,  6) 

semiconductors.  In  particular,  three  of  these  studies  '  ’  have  dealt 
with  silicon  impurities  in  GaAs.  It  was  found  that  one  could  correlate 
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many  of  the  observed  absorption  bands  with  the  presence  of  different 
point  defects  involving  silicon.  The  present  work  uses  this  identification 
to  investigate  the  influence  of  the  other  impurities,  specifically  tellurium 
and  zinc,  on  the  site  distribution  of  silicon. 

Review 

There  has  been  considerable  progress  in  the  understanding 

of  the  localized  vibrational  mode  absorption  in  silicon-doped  GaAs. 

Because  of  the  importance  of  previous  studies  for  the  interpretation  of 

the  present  results,  a  brief  review  will  be  presented  here.  This  review 

is  particularly  important  since,  at  the  time  of  this  writing,  one  of  the 

(6) 

more  important  papers  is  not  yet  published. 

Silicon,  as  an  isolated  substitutional  impurity  on  either  the 

gallium  or  arsenic  sublattice,  has  tetrahedral  point  group  symmetry. 

To  terms  quadratic  ui  the  displacement  of  the  silicon,  the  potential  energy 

is  that  of  a  spherically  symmetric,  harmonic  oscillator.  Because  the 

silicon  is  lighter  than  either  the  gallium  or  arsenic,  one  can  expect 

spatially  localized  vibrational  modes  associated  primarily  with  the 

motion  of  the  silicon  impurity  at  frequencies  above  the  highest  phonon 

(7) 

frequency  of  the  unperturbed  lattice.  A  single  silicon  defect  will  give 

three  vibrational  modes,  and  because  of  the  symmetry  of  the  potential, 

the  modes  will  be  degenerate.  Thus  a  single  infrared  absorption  baud 

is  exDected  for  Si_  and  another  one  for  Si.  .  The  two  sites  can  have 
Ga  As 

different  frequencies  because  of  a  slightly  different  mass  defect  and 
also  because  the  different  nearest  neighbors  in  the  two  cases  should 
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influence  the  effective  force  constants.  In  Table  I  it  is  seen  that  the 

observed  local  mode  frequencies  are  u>  =  384  cm  *  for  Si_  and  399 

ua 

cm  *  for  Si.  . 

As 


In  some  cases  the  silicon  impurity  is  no  longer  isolated 
from  ether  defects  but  .<>  paired  vith  another  silicon  or  another  impurity 
on  a  neighboring  site.  There  can  be  different  reasons  for  the  existence 
of  these  ion  rairs.  There  is  evidence  already  mentioned'  *  that 
f  Si  -  Si  j  pairs  tend  to  form  at  high  [  Si  ]  .  Moreover,  the  samples 
are  strongly  n-type  and  hence  the  resulting  free  carrier  absorption  must 
be  reduced  by  the  introduction  of  an  electrically  compensating  impurity. 
Lithium  or  copper  is  diffused  at  elevated  temperature  to  produce  the 
compensation.  Since  the  material  is  initially  n-type,  the  compensating 
diffusant  must  produce  acceptors  which  are  generally  assumed  to  be  Li^ 
or  Cu„  .  There  is  also  a  tendency  for  the  Si„  donors  and  the  Li„  or 
Cu^a  acceptors  to  interact  resulting  in  (^Si^-  Li_a^/  or  C^Ga*  ^UGa^ 
pairs.  Clearly,  in  these  cases,  the  silicon  is  no  longer  in  a  site  with 
tetrahedral  symmetry. 

-n 

For  the  (  Si_  -  Li_  )  pairs  where  the  defect  involves 
\  Ga  Gay 

impurities  on  second  neighbor  sites,  a  weak  coupling  model  ... 

(6) 

give  reasonable  agreement  with  experiment.  The  effect  of  one  atom  is 
to  alter  the  potential  of  the  other  atom.  Since  the  nearest  neighbor 


gallium  sites  are  along  a  <  110  >  direction,  which  is  not  a  symmetry 
rotation  axis  of  the  tetrahedral  point  group,  the  triply  degenerate  Si 
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band  will  be  split  into  three  bands  with  frequencies  nearly  independent 

of  the  lithium  isotope.  These  bands  will  be  located  near  the  Si„  fre- 

r  Ua 

~  1  (8) 

quency  of  384  cm  .  A  simple  perturbation  argument'  gives 


2 

(V 


l 

3 


3 

I 

i=l 


0). 


> 


where  u>  is  the  frequency  of  the  triply  degenerate  mode  and  the  uk  values 
are  the  perturbed  frequencies.  In  the  (^Si^-  Li^^1  defect  case  there 
will  also  be  three  bands  involving  primarily  lithium  motion.  The  fre¬ 
quencies  of  both  the  predominently  silicon  modes  and  those  due  to  lithium 
are  given  in  Table  I. 

The  ( Si^  Si.  ^  defect  is  more  complicated^’  ^  in  that 
V  Ga  As  y 

the  two  impurities  are  strongly  coupled,  and  both  are  light  compared  to 
the  atoms  they  replace.  The  pair  axis  is  in  a  <  111  >  direction  which  is 
a  three-fold  rotation  axis,  and  th<*  six  vibrational  modes  are  expected  to 
contain  two  sets  of  doubly  degenerate  modes.  Thus  there  can  be  four 
frequencies  of  which  three  are  identified  in  Table  I,  The  fourth  mode  is 

(^0 

either  too  small  in  absorption  to  be  observed  or,  as  calculation  indicates' 

falls  in  a  region  of  large  absorption  from  the  GaAs  host  lattice.  The 

location  of  all  bands  observed  in  lithium  or  copper  diffused,  silicon-doped 

GaAs  is  listed  in  Table  I.  For  completeness,  the  location  of  some  bands 

(10) 

related  to  lithium  paired  with  lattice  defects  in  undoped  GaAs  is  also 

given.  These  latter  bands  are  not  related  to  the  silicon  impurity  but  can 

(4) 

occur  in  al1  samples  which  are  lithium  diffused  at  sufficiently  high 
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Table  I 


Defect  Impurity 


Li-lattice 

defects 


Li 


SlGa‘  LlGa 


f 


Li 


Ga 


Si 


Ga 


SiGa"  °UGa 


Si 


Ga 


Si 


Ga 


Si 


Ga 


Si 


As 


Si_  -  Si  J>  < 
Ga  As  f  N 


Si 


As 


Si  -  SiA 
Ga  As 


Mode  Frequency  (cm  ) 


7  6 

Li  Li  Cu 


365 

379 

438 

447 

454 

374 

379 

405 


3  52 
389 
406 

470 

480 

487 

374 

379 

405 


374 

376 

399 


...  384  --- 


...  399  --- 

...  367  --- 

...  393  --- 

...464  --- 


The  frequencies  of  infrared  absorption  bands  in  lithium-  or  copper- 
diffused,  silicon  doped  GaAs. 
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temperature. 


There  have  been  a  numbei  of  experiments  which  have 
suggested  that  the  distribution  of  silicon  among  the  various  defects 
mentioned  will  be  influenced  by  the  presence  of  other  impurities.  Some 
electrical  and  radio-tracer  measurements  of  selenium  plus;  silicon-doped 
GaAs  indicated^  that  the  Se  donors  could  change  T)  =  [Si  ]/[Si  ] 

AS  AS  v_j£l 

from  <  1  to  >1.  Local  mode  measurements  of  a  sulfur  plus  silicon- 

(4)  -1 

doped  sample  showed  the  Si  band  at  3  99  cm  to  be  larger  than  the 

AS 

Si_  band  at  384  cm  In  similarly  silicon-doped  samples  but  without 
Ga 

sulfur the  384  cm  ^  band  is  ~  4-5  times  as  strong  as  the  3  99  cm 
band.  However,  the  bands  in  the  sulfur-doped  sample  were  small  and 
close  to  seme  very  large  Li-lattice  defect  bands,  thus  making  the  results 
uncertain.  Another  study  showed  that  silicon  site  transfer  could  take 
place^  as  a  result  of  lithium  diffusion  at  temperatures  >  800°C.  The 
sulfur  plus  silicon-doped  sample^  was  diffused  at  ~  800°C  again  making 
this  result  uncertain. 

Experimental  Results  and  Discussion 

Table  II  lists  several  ingots  of  GaAs  with  the  doping  levels 

anticipated  near  the  seed  end  on  the  basis  of  the  impurities  added  and 

published^^  distribution  coefficients.  Samples  were  taken  from  each 

ingot  and  diffused  with  lithium.  In  some  cases  samples  were  taken  from 

both  the  seed  end  and  near  the  back  end  of  the  ingot.  The  temperature  of 

diffusion  was  T  700°C  in  order  to  avoid  the  site  redistribution^  ^ 
d 
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Table  II 


1  Si  ~  1-2  x  10^®cm-^ 

2  Si  ~  2-3  x  lO^cm 

3  Si,Te  [Si  ]  ~  1  x  lO^cm  ^ ,  [Tel  ~  3-4  x  lO^cm  ^ 

4  Si,Te  [Si  1  ~  6-8  x  lO^cm  \  [Te]  ~  5-6  x  lO^cm  ^ 

18-3  18-3 

5  Si,  Zn  [  Si  ]  ~  6  x  10  cm  ,  [Zn]  ~  1  x  10  cm 

18-3  18-3 

6  Si,  Zn  [  Si  ]  ~  6  x  10  cm  ,  [  Zn]  ~  3-4  x  10  cm 

7  Si,  Zn  [  Si  1  ~  6  x  10*  \m  [  Znl  ~  3  x  10 *^cm  ^ 

GaAs  ingots  used  in  the  present  study  with  doping  concentrations 
near  the  seed  end  of  the  crystal 

- . - Q - - - 

observed  for  T^  >  800  C.  The  diffusion  was  done  from  a  surface  alloy 
phase  while  in  an  argon  atmosphere.  The  method  of  diffusion  has  been 
described  previously 

Samples  from  the  silicon -doped  ingots  1  and  2  were  lithium 
diffused  and  the  infrared  absorption  measured  at  liquid  nitrogen  tempera¬ 
ture  is  given  in  Fig.  1.  Also  shown  for  comparison  is  the  measured 
absorption  curve  for  pure  GaAs.  The  curves  show  the  spectral  region 
containing  most  of  the  local  mode  bands  which  are  labeled  with  the  related 
defect.  Although,  as  indicated  in  Table  I,  tome  bands  occur  at  uu  >  430 
cm  they  have  been  shown  in  previous  measurements^  ’  ^  and  will  not 
be  of  major  importance  for  the  present  work. 

The  (.  rosses  of  Fig.  2  show  the  peak  absorption  coefficient 

a  for  the  Si  _  band  at  384  cm  *  plotted  against  a  tc-r  the  Si.  band  at 
p  Ga  p  As 
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70  390  410  430 

cu(cm") 


nitrogen  temperature  absorption  of  samples  from 
1  and  2  of  Table  II.  The  samples  were  ^Li  diffused 
'C  for  24  hours. 


8- 


a  (  384  cm" 


CL  ( 399  cm' 1 ) 

Fig.  2.  A  plot  of  the  peak  absorption  coefficient  of  the  SiQa  band  at 
384  cm'*  vs  that  for  the  Si^g  band  at  399  cm'*.  The  x  are 
for  silicon  doped  samples,  the  o  tor  silicon  and  tellurium 
doped  samples,  and  u.  •  o  are  tor  silicon  and  zinc  doped 
Samples.  The  numbers  indicate  the  ingot  ot  Table  II.  All 
samples  are  lithium  diffused. 
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399  cm  ^ .  The  data  are  taken  from  a  number  of  silicon-doped  samples 

o 

which  are  all  diffused  with  lithium  at  <  800  C.  Most  of  these 

(6) 

samples  were  measured  as  part  of  an  earlier  etudy  .  The  points  are 
close  to  a  line  of  slope  ~  0.  9  indicating  that  the  Si ,  band  grows  slightly 

/i  8 

more  rapidly  than  the  Si  band  with  increasing  [Si  ]  .  It  is  to  be  empha 
sized  that  silicon  is  the  only  major  dopant  in  these  samples. 

The  absorption  of  samples  from  ingots  3  and  4,  which  are 

7 

doped  with  silicon  plus  tellurium  and  diffused  with  Li  is  shown  in  F\g.  3 

Comparison  with  Fig.  1  shows  several  new  features.  The  band*  near 

475  cm  *  and  391  cm  *  are  due  to  (  Li_  -  Te,  )  and  have  been  pre- 

viously  studied^^.  The  band  at  379  cm  *  i*  largely  a  Li-lattice  defect 

band  and  shows  a  large  lithium  isotope  shift,  i.e.  ,  379  cm  -  406  cm 

when  ^Li  “•  ^Li.  The  bands  at  464,  3  93,  and  367  cm  *  in  the  sample 

from  ingot  4  are  the  i  Si  ,  -  Si  ,  bands  mentioned  earlier  and  listed 

in  Table  I.  A  principal  difference  between  the  samples  of  Fig.  3  and 

Fig.  1  is  the  change  in  the  relative  strengths  of  the  Si  _  and  Si  bands 
®  Oa  As 

at  384  and  399  cm  .  The  absorption  peak  values  less  background  for 
these  bands  of  the  silicon  plus  tellurium  samples  of  Fig.  3  are  shown  as 
circles  on  Fig.  2  and  labeled  with  the  ingot  number.  The  effect  of  the 
tellurium  has  been  to  subs tanttally  enhance  the  fraction  of  the  (Si1 
which  resides  on  arsenic  sites.  More  extensive  measurements  should 
give  results  which  can  be  compared  with  thermodynamic  predictions. 

The  most  interesting  results  are  obtained  { r cm  samples 


.a. 
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taken  from  ingots  5,  6,  and  7  which  were  doped  with  silicon  plus  zinc. 

7 

Figure  4  shows  the  absorption  for  several  samples  diffused  with  Li. 

Figure  5  shows,  for  each  of  two  ingots,  a  comparison  of  two  samples 

T  & 

taken  from  adjacerf  wafers  and  diffused  with  Li  and  Li  respectively. 

/Q\ 

From  previous  studies v  ;  of  zinc -doped  and  lithium  compensated  samples 
it  is  known  that  there  are  four  lithium  bands,  presumably  associated  with 
(Li-Zn)  pairs.  These  bands  are  nearly  equal  for  a  given  [  Zn"!  and  are 
located  at  405,  378,  361 ,  and  340  cm  *  for  ^Li  and  at  433,  404,  385, 
and  361  cm  *  for  ^Li.  Many  of  these  bands  are  observed  in  Figs.  4  and 
5  and  are  taken  into  account  in  some  of  the  quantitative  discussion  which 
will  follow. 

The  Si  acceptor  band  at  399  cm  *  is  observed  only  in  the 

A  S 

most  lightly  zinc-doped  ingot  (5).  Within  experimental  accuracy  this 
band  ir  totally  absent  in  Fig.  4  from  the  samples  of  ingots  6  and  7.  We 
also  note  from  Figs.  4  and  5  that  the  {Si  -  Si  )  bands  are  totally  absent 

vjcL  A  S 

(at  464,  393,  and  367  cm  )  in  the  heavily  zinc-doped  ingots.  Moreover, 
several  of  the  samples  show  the  appearance  of  three  new  bands  of 

approximately  equal  strength.  These  bands  are  at  395,  ~382,  and  378 

-1 

cm 


The  378  cm  ^  band  requires  some  discussion.  In  Fig.  5, 

the  band  at  378  cm  *  could  involve  contributions  from  ?  absorption  bands 

7 

seen  in  previous  work.  In  CaAs  with  <  Li-Zn),  there  is  a  band  at  378 

-1  7  -1 

cm  ,  but  it  should  be  equal  in  strength  to  the  {  Li-Zn)  band  at  405  cm 

There  is  also  a  (^Li-lattice)  defect  band  at  379  cm  *  which  could  be 
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7  7 

contributing.  However,  both  the  (  Li-lattice)  and  the  (  Li-Zn)  bands 

“  1  6  7  6 

shift  to  near  405  cm  if  Li  is  used.  Comparison  of  the  Li  and  Li 

- 1 

samples  of  ingot  7  in  f’ig.  5  shows  little  change  in  the  378  cm  band. 

The  remaining  previously  known  contribution  to  the  absorption  near  378 

cm  *  comes  from  a  379  cm  *  silicon  mode  for  (Si _  -  Li_  )  defect  (see 

Ga  Ga 

(6) 

Table  I).  However,  this  contribution  should  only  be  as  large  as  the 

much  weaker  band  at  374  cm  *  (see  Fig.  5,  ingot  7).  Therefore,  we 

-1 

have  the  result  that  most  of  the  strength  of  the  378  cm  band  of  ingot  7, 
Fig.  5  cannot  be  accounted  for  in  terms  of  any  of  the  previously  known 
defects. 


The  enhanced  strength  of  the  378  cm  band  occurs  only  in 

those  cases  where  the  395  and  382  cm  *  bands  are  also  observed.  The 

three  new  bands  do  not  show  any  lithium  isotope  shift,  and  they  occur 

only  in  samples  with  large  [  Si  3  and  [  Zn  1.  It  is  reasonable  to 

assume  that  these  new  bands  arise  from  (Si_  -  Zn_  )  defects  where  the 

Ga  Ga 

silicon  and  zinc  are  on  nearest  neighbor  gallium  sites  (second  neighbor 

positions).  or  such  a  pair,  the  axis  is  <  110  >  and,  as  discussed  for 

(Sl_  -  Li_  ),  the  three  fold  degeneracy  of  the  tetrahedral  potential  will 
Ga  Ga 

be  completely  lifted.  Since  the  mass,  M(Zn)  =  65,  is  close  to  theM(Ga) 
70  which  it  replaces,  we  expect  to  see  only  three  silicon  modes  split  in 
frequency  about  the  Si  band.  Application  of  the  frequency  rule  given 

\-l  <* 

previously  yields 


U) 


=  {  y  [(395)2+  (382)Z+  (378)21  |  =  384.- 


cm 


-1 
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This  result  is  in  good  agreement  with  the  384  cm  ^  value  for  the  Si 


Ga 


band. 


The  concentration  of  (Si  -  Zn  )  may  be  crudely  estimated 

Oct  Gel 

as  follows:  It  is  assumed  that  the  total  absorption  cross  section  for  the 

{Si _  -  Zn  )  bands  is  the  same  as  that  for  the  Si  band,  i.  e.  ,  the  only 
ua  ua  Ga 

effect  of  the  zinc  has  been  to  lift  the  degeneracy  of  Si„  but  the  total 

Ga 

[SiGa-  ZnGa^ 

absorption  per  center  is  unchanged.  Then,  — - r- -  = 

[  Si  J 


S395  +  a382  *  a378 
a  384 


Ga 


.  The  prime  on  i®  1°  specify  that  the  other 


sources  of  absorption  at  that  frequency  have  been  already  subtacted. 
From  other  studies^"’  ^  it  is  known  that  the  absorption  cross  section  for 

■I  - 1 8  2 

the  3 t»4  cm  band  is  given  by  a,_  . /[  Si_  Is”  7.3  xlO-  cm  .  Thus 

384  Ga * 

19-3 

[Si^a~  Zn^  ]  for  the  ingot  7  samples  of  Fig.  5  is  ~  3  x  10  cm  and 


Ga 


r  .  t  19-3  19 

LSi^J  —  4-5  x  10  cm  .  The  total  silicon  concentration  of  7-8  x  10 

-3  19  -3 

cm  compares  favorably  with  the  doping  estimate  of  ~  6  x  10  cm 


given  in  Table  II.  The  total  zinc  concentration  of  ingot  7  is  ~  4-5  x  10 


19 


-3  19  -3 

cm  .  This  concentration  includes  both  [  Zn_  -  Si_  ]  ~  3  x  10  cm 

Ga  Ga 

and  [  Zn_  -  Li_  ]~  102  x  lO^cm  ^  as  estimated  from  the  405  cm  * 

Ga  Ga 

■  1  (8 ) 

and  361  cm  bands  and  previous  published'  7  a  vs.  [  Zn_  -  Li  ]  data. 

Ga  Ga 

Again  the  total  [  Zn]  is  regarded  as  in  reasonable  agreement  with  the 
19  -1 

estimate  of  3  x  10  cm  given  in  Table  II.  It  may  also  be  noted  that  if 


the  Zn„  and  Si_  were  randomly  distributed  on  the  gallium  sublattice 
Ga  Ga 

17  18  -3 

then  the  [  Zn_  -  Si  ]  is  10  to  10  cm  for  these  samples.  This 
Ga  Ga 

concentration  is  substantially  smaller  than  the  estimate  given  above. 
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The  measurements  of  samples  taken  from  ingot  6  and  7  have 
been  used  to  obtain  the  points  given  by  solid  'lots  on  Fig,  2.  Recall  that 
in  these  samples  no  detectable  band  was  observed  at  399  cm  * ,  and  thus 
the  arrows  indicate  that  the  points  have  been  placed  at  the  maximum 
possible  value  for  a  (399  cm  *)•  The  true  value  is  probably  substantially 
smaller.  The  points  indicate  that  for  ingots  6  and  7,  the  [  Si  ]  has 

A  0 

been  decreased  by  at  least  one  order  of  magnitude,  and  a  more  realistic 

estimate  gives  ~  two  orders  of  magnitude  for  a  given  [Si  ].  Although 

Oa 

not  plotted  in  Fig.  2,  a  comparison  of  Fig.  1  with  Figs.  4  and  5  gives  a 

similar  conclus ion  for  [  Si,„  -  Si.  ]. 

Ga  As 

It  is  reasonably  clear  that  studies  of  silicon  site  distribution 
as  functions  of  [Si]  ,  [  Te],  and  [  Zn]  should  lead  to  quantitative  esti¬ 
mates  of  distribution  coefficients  and  the  dependence  of  specific  defects 
on  impurity  concentrations.  Moreover,  correlations  with  electrical 
measurements  should  yield  information  on  the  electrical  nature  of  each 
of  the  defects.  This  work  is  now  being  carried  out. 
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1.1.3  Thin  Film  and  Interface  Phenomena 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
C.  R  Crowell 

A  program  of  device-oriented  studies  related  to  current 
transport  and  charge  storage  in  metal-semiconductor  {Schottky)  contacts. 


1.  1.3.  1  S  rface  State  and  Interface  Effects  on  the  Capacitance- 
Voltage  Relationship  in  Schottky  Barriers 
C.  R.  Crowell  and  G.  I.  Roberts 


The  analysis  given  in  the  previous  report  has  been  completed 
and  submitted  for  publication.  A  further  modification  gives  the  relation¬ 
ship  between  the  surface  state  charge  density,  o,  in  zero  temperature 
approximation  and  for  a  finite  temperature  T.  The  precise  result  is 
shown  to  be 


2  , 

d2n 

n  IE  A  + 
as  t 

Tz  |kT> 

88 

dE2 

+  8T4  +  •  • . 

Ef 

where  n  (Ef)  is  the  surface  state  density  per  unit  area  per  unit  energy 

8  S  I 

at  the  Fermi  energy  E  .  Thus  as  long  as  the  curvature  of  the  n  versus 

I  8  8 

i  *2 

E  relationship  at  the  Fermi  energy  is  small  compared  with  n  (E  )  |  kT  |  , 

8  8  I 

the  zero  temperature  approximation  will  be  satisfactory. 


The  results  of  this  analysis  will  be  used  as  a  tool  for  further 
experimental  studies  of  surface  states  in  metal  semiconductor  systems. 
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1.  1.  3.2  Thermionic-Field  Resistance  Maxima  in  Metal- 
Semiconductor  (Schottky)  Barriers 
C.  R.  Crowell  and  V.  L.  Rideout 


A  maximum  in  the  differential  resistance  versus  applied  bias 
relationship  of  metal-semiconductoi  contacts  is  predicted  to  occur  when 
current  flow  is  predominantly  by  thermionic -field  (thermally  excited 
tunnel)  emission.  The  predicted  resistance  peaks  are  generally  asym¬ 
metrical  with  respect  to  voltage  and  may  occur  on  either  side  of  zero 
bias.  The  peak  location  has  only  an  indirect  correlation  with  the  Fermi 
kinetic  energy  in  the  semiconductor.  The  theoretical  approach  is 
generally  applicable  to  any  metal-semiconductor  system  when  the 
dominant  carrier  flux  is  associated  with  the  tail  of  a  Fermi-Dirac 
distribution.  The  theory  is  in  reasonable  agreement  with  recent  experi¬ 
mental  resistance  measurements  on  Cr-Si  Schottky  barrier  diodes  at 
77°K. 


Anomalies  in  the  differential  resistance  versus  applied  bias 

.  1-5 

relationship  near  zero  bias  have  been  observed  in  p-n  junctions  , 

6-8  9-12 

metal-insulator-metal  or  semiconductor  structures,  and  metal- 

13- 19  9 

semiconductor  (Schottky)  barriers  .  In  addition  Gray  ,  Esaki  and 

Stiles^,  and  Chang^°  have  used  the  WKB  tunneling  approximation  for 

the  insulator  barrier  to  calculate  the  incremental  resistance  of  metal- 

2  1 

insulator-semiconductor  structures.  Stratton  and  Padivani  have 
similarly  applied  the  WKB  expression  to  Schottky  barriers.  The  present 
report  is  concerned  with  the  differential  resistance  associated  with 
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thermionic -field  (thermally  excited  tunnel)  emission  currents  in  metal- 


semiconductor  contacts. 

22 

Recently  Crowell  ana  Rideout  (C-R)  have  analyzed  the 
thermionic -field  current- voltage  characteristics  of  Schottky  barriers  in 
a  normalized  formulation  which  uses  the  WKB  tunneling  approximation. 
From  consideration  of  the  principles  of  detailed  balance  and  reciprocity 
of  tunneling  coefficients,  C-R  showed  that  the  current  density,  J,  versus 
applied  bias,  V^,  relationship  for  thermionic -field  emission  (i.e.  , 
current  flow  dominated  by  the  Maxwell- Boltzmann  tail  of  Fermi-Dirac 
distributions  in  both  semiconductor  and  metal)  is  of  the  form 

J  =  J  exp  (qV  / nkT ) [  1  -  exp  (-qV  /kT)j  (1) 

8  1  1 

(cf.  ,  equation  (23)  of  C-R).  J  is  the  saturation  current  density  and 

8 

n  =  (q/kT)  dV^/dln  J^,  where  J  is  the  forward  component  of  the  carrier 
flux  (J^  J  for  »  kT/q).  Equation  (1)  will  not,  however,  be  valid 
when  pure  field  emission  becomes  dominant  at  high  semiconductor 

carrier  density  and  low  temperatures.  Band-tailing  density  of  states 

13-15  .  .  ,  ,  .16  J  t 

effects  ,  many  body  electron-po.aron  interactions  and  phonon 

17,18  L  .  ...  , 

excitations  ha  e  been  proposed  as  responsible  for  resistance 

anomalies  observed  in  metal-degenerate  semiconductor  contacts. 

Although  such  effects  are  not  considered  per  se  in  this  report,  equation 

(1)  is  generally  applicable  to  any  me  tal  -  s  emic  onduc  tor  contact  when  the 

dominant  carrier  flux  is  associated  with  th^.  .ail  of  a  Fermi-Dirac 

die  t.-ibution. 
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From  equation  (1)  the  incremental  resistance,  R  =  dV^./dJ, 

normalized  by  the  incremental  resistance  at  zero  bias,  R  =  kT/qJ  ,  is 
1  0  8 

R/R  =  n  exp  (-qV,/nkT)[  1  +  (n-l)exp  (-qV./kT)]'1  .  (2) 

o  i  * 

The  applied  voltage,  ,  at  which  a  resistance  maximum  (cond  ictance 

2  2 

minimum)  occurs  (i.e.  ,  when  d  J/dV^  =  0)  is 

V,  -  (2kT/q)  In  (n-1)  .  (3) 

im 

In  equations  (2)  and  (3)  we  have  assumed  that  there  is  only  a  slow  varia¬ 
tion  of  the  parameters  J  and  n  with  V,. 

r  8  f 

Figure  1  shows  the  normalized  incremental  resistance, 

R/R  ,  as  a  function  of  the  normalized  applied  voltage,  qV  /kT,  for 
o  t 

selected  values  of  n  (cf.  equation  (2)).  The  pip  mark  on  each  curve 

indicates  the  value  of  qV,  /kT  (cf,  equation  (3)).  In  general  the 

fm 

resistance  peaks  do  not  occur  exactly  at  zero  bias  nor  are  they  symme¬ 
trical  with  respect  toqV^/kT.  Equation  (2)  does,  however,  have  the 
following  symmetry  characteristics:  a  curve  with  a  maximum  at  some 
negative  value  of  and  defined  by  n  -  n  has  a  mirror  image  reflected 
in  the  zero  bias  line  it  n  :  n+  where  1/n  =  i-i/n+  (cf.  Figure  1).  For 

values  of  n  >  2  the  metal- aemic  onduc  tor  contact  behaves  like  a  backward 
diode.  When  n  2  the  resistance  peak  is  symmetrical  about  zero  bias: 

R/R  |  =  L  cosh  (qV ,/2kTp  "l  .  (4) 

o  1  n  2  f 
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Normalized  incremental  resistance  R/'R^  versus  normalized  applied 
bias  qVj/kT  for  selected  diode  n  values.  The  pip  mark  on  each  curve 
indicates  the  position  of  the  resistance  maximum  as  defined  by- 
equation  (3).  Note  that  each  curve  has  an  image  reflected  across 
the  zero  bias  line. 


Once  any  experimental  R-V^  characteristic  has  been  fitted 

to  equation  (2).  the  subsequent  degree  of  agreement  w i th  any  detailed 

theorv  of  thermionic-field  emission  depends  upon  the  accuracy  with 

which  the  parameters  n  and  R  {or  J  )  can  be  correlated  with  other 

o  & ' 

measurements.  C - R  have  calculated  n  and  j  as  a  function  of  RT/'E 

s  oo 

with  normalized  bandbending  E^/kT  as  parameter  (cf.  Figures  6  and  7 

16  21 

of  C-R).  The  term  E  ,  a  constant  of  the  material  ,  is  associated 

oo 

with  the  one-dimensional  single-particle  WKB  expression  for  the  trans¬ 
mission  of  the  barrier  for  carriers  at  the  bottom  of  the  conduction  band 
(cf.  equation  (8)  of  C-R): 

E  =  (qh/4n)  (N/m*c)*  =  18.  5  x  10* 12 (N/m  £  )*  eV  ,  (5) 

co  r  r 

where  N  is  the  semiconductor  doping  concentration  in  units  of  cm  \  m 

r 

the  tunneling  effective  mass  measured  in  units  of  the  free  electron  mass, 
and  c  the  semiconductor  static  dielectric  constant.  The  bandbending  in 
the  semiconductor  depletion  region  in  units  of  kl  is 

Eb/kT  =  (q/kT»(cpb-epg.  Vf)  ,  (6) 

where  cpb  is  the  metal-semiconductor  barrier  height,  and  cpg  is  the 
difference  between  the  semiconductor  Fermi  Level  and  the  bottom  of  the 
conduction  band  (<pg  is  negative  for  degenerate  semiconductors).  The 
quantity  kT/E  determines  the  magnitude  of  the  kinetic  energy  relative 
to  the  bandbending  Eb  at  which  the  maximum  injec  ;on  of  carriers  through 
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the  barrier  occurs  (cf.  equation  (14)  and  Figures  3  and  4  of  C-R).  Thus 
kT/E^  indicates  whether  the  predominant  character  of  the  current  flow 
is  thermionic,  iherrnionic-field,  or  field  emission. 


Figure  2,  deduced  from  Figures  7  and  11  of  C-R,  shows  the 

close  correlation  of  the  diode  n  value  with  kT/'E  .  Larger  values  of  n 

oo 

and  would  be  predicted  from  theories  which  added  effects  of  image 
24 

force  lowering  ,  a  two  band  model  for  the  complex  k  portion  of  the 
15,25 

E-k  relationship  ’  ,  an  interfacial  layer  between  the  metal  and  the 

26,27 

semiconductor  ,  or  an  increased  semiconductor  surface  doping. 

12 

For  resonant  tunneling  phenomena  ,  n  and  are  sufficiently  strong 
functions  of  V^.  that  equation  (1)  s'.ould  not  be  expected  to  be  an  appro¬ 


priate  form  for  analysis  of  the  experimental  data.  It  would  also  not  be 

desirable  to  apply  the  C-R  analysis  to  cases  where  (n-1)  <  0.  02  for 

24 

which  image  force  lowering  is  important  ,  or  where  n  »  1  for  which 

field  emission  is  dominant.  C-R  have  derived  criteria  for  the  bias 

range  for  the  transition  from  thermionic -field  emission  to  field  emission. 

For  reverse  bias  voltages  less  than  V  ,  the  maximum  of  the  energy 

r-ma.4  07 

distribution  of  emitted  carrier  flux  occurs  above  the  Fermi  level  in  the 
metal: 


V 

r-max 


V’b 


sinh  (E  /kT)  . 
'  oo 


(7) 


i  or  degenerate  semiconductors  a  corresponding  maximum  forward  bias 
voltage  is  given  by 


f-max 


V  sinh2(E  /kT)  . 

r-max  oo 


(8) 
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kT/Eoo 


Fig.  2.  Barrier  constant  (n-1)  versus  normalized  materials  parameter 

kT/E00  for  selected  values  of  the  normalized  bandbending  E^/kT. 
The  right  hand  ordinate  indicates  the  bias  in  units  of  kT/q  for 
occurance  of  the  resistance  maximum. 
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Saxena  has  recently  measured  the  incremental  resistance 

in  the  neighborhood  of  zero  bias  for  Cr  n-type  Si  Schottky  diodes  under 

conditions  where  appreciable  contributions  from  T-F  emission  could  be 

anticipated.  His  experimental  data  at  ?7°K  are  shown  in  Figure  3.  From 

equation  (2)  n  «  (-q/kT)dVf/d  In  R  for  qVf/kT  »  1.  From  Saxena's 

data  we  obtain  a  value  for  n  of  5.  7.  For  E^kT  90  and  n  =  5.  7, 

Figure  2  gives  kT/E  =  0.18  which  corresponds  to  N  =  1.2  x  1019cm"3 

z\:  :  (  K  {cl.  equation  (5)},  The  hulk  doping  concentration  of  the  Si  used 

by  Saxena  as  determined  from  Hall  coefficient  measurements  was 
18  3 

NHall  =  7  X  10  cm  .  The  doping  concentration  throughout  the  wafer 
may  of  course  differ  from  the  average  measured  by  the  Hall  effect.  In 
addition,  the  local  surface  doping  concentration  may  be  different  from 
^Hall*  ^  possible  that  the  Cr-Si  diodes  contain  an  interfacial 

layer  which  has  resulted  in  an  increased  diode  n  value. 

When  Rq  and  n  of  equation  (2)  are  chosen  to  match  Saxena  s 
data  for  V{  »  kT  as  shown  by  the  dashed  line  in  Figure  3,  the  experi¬ 
mental  data  are  matched  at  the  peak  of  the  T  F  curve.  This  occurs  at 

9V  /kT  =  3.  1.  The  barrier  height  associated  with  the  value  of  R  is 

o 

0.61  eV,  in  agreement  with  that  typically  observed  by  Saxena  for  diodes 
of  low  doping  concentration.  The  value  of  the  barrier  height  used  by 
Padovani  and  Stratton  is  not  known.  »  he  disagreement  at  low  bias 
between  Saxena's  experimental  data  and  the  prediction  of  the  T-F  model 
of  C-R  is  expected  since  the  criteria  for  T-F  emission  implied  by  equa¬ 
tions  (7)  and  (8)  can  only  be  satisfied  by  non-zero  V  and  Vr 

r-max  f-max 
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is  the  appropriate 


o  1  o 

above  120  K  (assuming  that  N  T  ,=7x10  cm 

Hall 

doping).  A  sufficiently  large  percentage  of  the  tunnel  current  flux  arises 
from  carriers  above  the  Fermi  energy  that  the  form  of  equation  (3)  is  a 
good  approximation  for  moderate  forward  bias.  The  experimental  curve 
near  zero  bias  is  unquestionably  characteristic  of  Fermi-Dirac  distri¬ 
butions.  Note  that  the  T-F  emission  theory  fits  as  well  at  77°K  as  the 

23 

zero  temperature  calculation  of  Padovani  and  Stratton  .  The  shapes  of 
the  R-Vj  characteristics  for  the  two  models  are  even  remarkably  similar 
considering  that  the  Maxwell  -  Boltzmann  approximation  in  the  model  of 
C-R  does  not  contain  the  inflection  point  in  the  occupation  probability  at 
the  Fermi  level.  It  is  apparent  that  measurements  similar  to  those  of 
Saxena  over  a  wide  range  of  temperature  and  semiconductor  doping  would 
be  of  interest. 

We  wish  to  thank  Dr.  A.  N.  Saxena  for  permission  to  quote 
his  experimental  results  prior  to  publication. 
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1.1.4  Radiative  Recombination  in  Semiconductors 


NCR -05-018-044,  National  Aeronautics  and  Space 
Administration 

GKF  1764,  National  Science  Foundation 
AF-AFOSR-496,  Joint  Services  Electronics  Program 
DA -A RO-D-3 1  -  124-G1 054 ,  U.S.  Army  Research  Office, 
Durham 

M.  Gershenzon 

General  Objectives 

1.  To  determine  the  mechanisms  of  radiative  and  of  non-radiative 
recombination,  particularly  in  wide  band-gap  semiconductors. 

2.  To  correlate  extrinsic  recombination  processes  with  the  defects 
responsible  for  them. 

3.  To  survey  a  number  of  relatively  unexplored  wide  band-gap 
materials,  seeking  efficient  luminescence  mainly  in  the  visible 
spectrum. 

4.  To  study  competing  recombination  processes  as  a  function  of 
temperature. 

5.  To  investigate  luminescent  devices  where  the  minority  carriers  are 
created  by  electrical  injection. 

Progress  and  Plans 

1.  GaN 

In  the  last  report  we  described  various  vapor  transport 


methods  which  produced  small  crystallites  of  GaN  and  indicated  the 
feasibility  of  epitax'al  deposition.  In  order  to  grow  larger  platelets 
required  for  the  optical  work  we  attempted  to  reduce  the  growth  varia¬ 
bles  by  using  the  material  grown  previously  as  the  source  in  a  closed- 
tube  halide  (chloride)dispropor donation  vapor  transport  technique.  After 
optimizing  temperatures,  temperature  gradients  and  halide  concentration 
to  yield  re  atively  large  crystals,  it  was  discovered  that  these  were  in 
reali  •  B  Ga^O^  rather  than  GaN  which  was  introduced  as  a  contaminant 
in  the  source  material.  Thus  we  have  recently  prepared  uncontaminated 
source  nia  arial  by  a  direct  open  tube  liquid  Ga-NH^  reaction  and  are 
repeating  the  optimization  process  for  single  crystal  preparation  via  the 
closed  tube  ha  de  transport  process. 

Our  future  plans  remain  the  same  as  before,  namely  the 
study  f  the  optical,  electrical  and  luminescent  properties  of  the  "pure" 
material  as  wall  as  the  effects  of  doping  with  elements  which  should  act 
as  shallow  donori  or  acceptors  or  as  luminescent  centers. 

2.  A1P 

In  the  last  report  we  described  r  closed  tube  iodide  transport 
method  which  transported  A1P  source  material  obtained  by  several  methods 

v  my  rapidly.  Now  this  process  has  been  optimized  to  some  extent 

c  o 

t  T  ~  1100  ,  AT  -  40  )  to  produce  fiat  ribbons  and  wh  s kers  up  to  1  cm 
in  'engti  These  are  large  enough  for  luminescence  measurements  but 
too  narrow  for  absorption  work.  The  student  working  on  this  project 
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has  left  the  laboratory.  Hence  this  project  will  be  delayed  until  a 
replacement  is  found  for  him. 

3.  Si 

Our  plans  to  look  for  isolated  donor-acceptor  pair  recombi¬ 
nation  emission  lines  are  proceeding.  Double  doped  (B,  P)  Si  ingots 
were  grown  as  described  in  the  last  report.  A  high  resolution  grating 
spectrometer  (Perkin-Elmer  El)  is  being  calibrated  and  equipped  for 

photoluminescence  measurements  aimed  at  measuring  and  optimizing 

o 

the  broad,  distant  donor-acceptor  pair  band  at  4.2  K  prior  to  exciting 
a  high  power  CW  GaAs  laser  in  the  attempt  to  saturate  the  distant  pairs 
and  allow  the  closer,  isolated  pair  lines  to  be  obse  -'ed. 

4.  ZnS:  Cu 

We  hope  to  elucidate  the  mechanism  of  the  very  efficient 
green  emission  of  Cu  activated  ZnS  at  room  temperature  by  exciting  this 
emission  at  4.2°K  with  very  high  intensity  radiation  in  its  characteristic 
excitation  absorption  band  and  looking  for  any  changes  in  the  emission 
spectrum,  notably  the  possible  appearance  of  lines.  A  pulsed  molecular 
laser  (average  cw  power  100  mw)  whose  output  wavelength  matches  the 
excitation  band  very  well  has  been  set  up.  ZnS  crystals  obtained  from 
various  sources  are  being  doped  with  Cu  and  as  co-activator  A1  or  a  halide 
A  quartz  prism  monochromator  (Perkin-Elmer  Model  was  calibrated 
and  is  being  used  for  photolummesc ence  measurements  to  establish  the 
presence  and  efficiency  of  the  Cu-green  emission.  It  is  planned  to 
excite  the  efficient  samples  with  the  focussed  laser  beam  looking  for 

l 
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sharp  lines  in  the  time  resolved  spectra  using  a  high  resolution  grating 
instrument  (Bausch  and  Lomb,  2  Meter)  and  integrating  the  signal  either 
via  photographic  recording  or  via  a  photomultiplier  and  a  boxcar  integrator. 

5,  Cubic  ZnS 

We  have  obtained  some  platelets  of  ZnS  having  large  areas  of 
the  zincblende  structure  which  are  relatively  free  of  stacking  faults.  We 
hope  to  study  the  absorption  edge  and  luminescence  due  to  free  exciton, 
bound  exciton  and  "edge  emission"  in  such  cubic  crystals  utilizing  the 
high  resolution  capabilities  of  the  Bausch  and  Lomb  spectrograph  and 
invoking  uniaxial  stress  and  Zeeman  splittings  where  possible.  For  this 
work  a  quartz  dewar  is  being  constructed  with  capabilities  for  absorption, 
photoluminescence  and  photoexcitation  work  in  conjunction  with  the  high 
resolution  s  pec  t r om e te  r . 
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1.  1.  5  Electron  Tunneling  in  Soiids 


R.  T.  Payne 

Our  technology  has  produced  a  large  number  of  electronic 
solid  state  devices  which  rely  on  the  quantum  mechanical  phenomena  of 
tunneling  (the  ability  of  a  particle  to  go  through  a  barrier  rather  than 
over  it).  Typical  examples  of  these  devices  are  p-n  diodes,  Esaki  or 
Tunnel  diodes,  Schottky  diodes,  Gunn  Effect  oscillators,  MOS  or  Cold 
Cathode  structures,  superconducting  junctions,  and  Josephson  junctions. 

Previous  work  by  the  author  indicated  that  by  using  low 
temperatures,  high  magnetic  fields,  and  high  pressures,  these  electron 
tunneling  devices  could  be  uP**d  as  probes  to  study  the  properties  of  the 
solid  lattice  and  properties  of  the  electron  gas  in  the  solid.  Other 
investigators  have  shown  potential  uses  for  these  devices  in  a  variety 
of  fields,  such  as  microwave  spectroscopy,  ultrasonics,  infrared 
spectroscopy,  and  low  energy  electron  spectroscopy. 

In  October  1966,  the  Electron  Tunneling  Laboratory  was 
established  to  use  these  electron  tunneling  devices  in  the  study  of  the 
fundamental  properties  of  solids,  solid  surfaces,  and  general  physics. 
The  laboratory  will  eventually  have  facilities  lor  high  magnetic  fields 
(100  kilo=oersted),  high  pressure  (60,000  psi),  and  low  temperatures 
(1 . 0°Kelvin).  It  will  require  a  variety  of  measuring  devices,  a  multiple- 
derivative  plotter,  a  susceptibility  bridge,  microwave  sources, 
receivers  and  frequency  standards,  a  nano-volt  potentiometer  and 
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voltaee  references. 

To  date,  the  mulci-derivative  plotter,  high  pressure  system 
(45,000  psi),  low  temperature  system  to  '  3°K,  medium  duty  evaporator, 
and  cold  cathode  sputtering  module  have  been  put  into  operation.  These 
systems  are  being  used  in  the  research  to  be  discussed  in  the  following 
sec  tions. 

1 .  i .  5,  1  Study  for  a  Measurement  of  Galactic  and  Stellar 

Emission  Spectra  in  the  500  Micron  Region.  Part  II 
NGL-05-018-044  Suppl.  #2,  National  Aeronautics 
and  Space  Adminstration 
R.  T.  Payne  and  F.  Schrey 

This  is  a  contract  to  study  wide  band  detectors  for  use  in  the 
500  micron  region.  One  such  detector  is  a  Josephson  Junction.  Since 
transition  metal  alloys  provide  the  highest  transition  temperature  for  the 
Josephson  effect  these  materials  are  under  study.  Before  the  contract 
per’"1  '!  J'u1'  1068-'  ' t  w^>«  ’  t*>  ferm’nrfr  o. 

the  preparation  of  the  alloy  sputtering  source  was  unsuitable.  Apparently, 
the  composition  of  the  alloy  was  not  close  enough  to  the  Nb^Sn  crystal 
phase  to  produce  superconducting  films.  Further,  the  shape  of  the 
source  made  it  difficult  to  use.  Examination  of  the  material  is  in 
progress. 

At  the  beginning  of  the  contract  period  a  tantalui.  metal 
source  was  obtained  for  sputtering  experiments  at  R.  D.  Mathis  Co. 

\ 
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Mr.  Schrey  commuted  (and  still  Hops)  to  the  company  to  sputter  tantalum 
films  on  glass  substrates  for  study.  The  resulting  films  look  satisfactory 
and  have  resistivities  close  to  that  of  the  bulk  material.  Cryogenic  tests 
of  these  films  are  in  progress.  A  V^Ga  metal  source  which  has  been 
on  order  since  April  should  be  arriving  soon.  This  source  will  be  used 
to  study  the  sputtering  of  alloy  films  in  place  of  the  Nb^Sn. 

1.  1.  5.2  Excitonic  Insulators  in  Semiconductor  Junctions 
AF-AFOSR-496- 67,  Joint  Services  Electronics 
Program 
R.  T.  Payne 

The  excitonic  insulator  is  a  many  body  state  similar  to  that 
of  a  superconductor.  The  author  discovered  the  state  while  investigating 
the  properties  of  semiconductor  junctions.  Experiments  are  still  in 
progress  to  discover  the  transition  temperature  dependence  on  pressure, 
carrier  concentration,  and  impurity  atoms. 

Another  possible  material  for  study  is  chromium.  During 
the  summer  preliminary  tests  on  chromium  junctions  were  made.  No 
results  have  been  obtained  to  date.  Further  experiments  are  being 
encouraged  at  the  Naval  Ordinance  Laboratory  at  China  Lake,  California. 
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1 .  1 .  5.  3  A  100  Kilo-oersted  Supercouduct  ig  Magnet  Facility 


AF-AFOSR -4%-67,  Joint  Serve  es  Electro. lies 
Program 

M.  H.  Halloran  and  R.  T.  Payne 

The  superconducting  magnet  facility  has  been  started  jointly 

with  M.  H.  Halloran.  A  magnet  has  been  ordered  which  is  capable  of 

4 

producing  a  100  Kilo-oersteds  with  a  homogeneity  of  4  parts  in  10  in  a 

5 

1  cm.  sphere,  a  stability  of  1  part  in  10  and  a  total  working  volume  of 
3.810  cm.  in  diameter.  Using  this  facility,  experimental  research  will 
be  pursued  in  the  general  area  of  electronic  properties  of  metals  and 
degenerate  semiconductors  at  low  tempe  ture. 

The  specific  scientific  goals  of  this  research  are  five-fold: 
to  obtain  complete  Fermi  surface  information  for  the  transition  metal 
elements  and  certain  intermetallic  compounds  (Halloran),  to  determine 
whether  the  theoretical  explanation  of  the  large  B-H  effects  observed  in 
the  noble  metals  holds  true  (Halloran),  to  determine  whether  the  effec¬ 
tive  mass  of  electrons  in  a  metal  is  magnetic  field  dependent  as  suggested 
recently  on  the  basis  of  field  theory  consideration  (Halloran),  to  measure 
pressure  effects  of  the  Fermi  surface  of  transition  metals  and  compare 
these  effects  with  theories  of  critical  spin  fluctuations  (Payne),  and  to 
study  the  properties  of  electrons  tunneling  through  semiconductor  inter¬ 
faces  in  the  presence  of  intense  magnetic  fields  and  high  pressures  (Payne). 

After  repeated  attempts  to  negotiate  with  the  suppliers  of  the 
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magnet  to  obtain  a  properly  working  magnet  the  prospective  users 
decided  to  change  vendors.  Negotiations  are  now  being  made  with 
different  suppliers  for  the  same  type  of  magnet.  All  research  requiring 
the  high  magnetic  field  has  been  suspended  until  a  new  supplier  is  found. 

1.  1.  5.4  Magneto-Tunneling  in  Semiconductors  under  High 
Pressure 

AF-AFCSR-4%-67,  Joint  Services  Electronics 
Program 

R.  T.  Payne  and  J.  Downs 

It  is  planned  to  measure  magneto-tunneling,  the  magnetic 
field  dependence  of  the  tunneling  current  of  semiconductor  p-n  junctions, 
in  the  presence  of  hydrostatic  pressure  and  at  low  temperatures.  The 
change  in  effective  mass  with  pressure,  the  change  in  quasi-particle 
"g"  factor  with  pr  ssure  and  other  properties  of  the  semiconductor  may 
be  derived  from  this  experiment  and  compared  with  theory.  Preliminary 
measurerm  .its  of  Ge  back  diodes  have  been  made  to  32,000  psi  at  room 
temperature.  Good  comparison  has  been  made  with  the  work  of  M.  E. 

Sikorski  and  P.  Andreath. 

The  summer  was  spent  revising  the  high  pressure  equipment. 

Present  operation  is  now  extended  to  45,000  psi. 
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1.  1.  5.  d  Pressure  Dependence  of  the  Magnetic  Susceptibility, 
the  Fermi  Surface,  and  the  Effective  Mass  of 
Palladium 

AF-AFSR-496-67,  Joint  Services  Electronics 

Program 

R.  T.  Payne 

It  is  planned  to  measuie  the  magnetic  susceptibility,  the 
effective  mass,  and  several  extremal  cross  sections  of  the  Fermi  surface 
of  palladium  as  a  function  of  pressure.  These  measurements  will  be 
compared  with  similar  magnetic  measurements  on  alloys  of  palladium 
and  with  the  theory  of  Berk  and  Schrieffer  on  critic^1  spin  fluctuations, 

A  better  understanding  of  the  properties  of  palladium  and  of  ferro¬ 
magnetic  exchange  may  be  obtained. 

At  present  a  palladium  single  crystal  has  been  bought. 
Preliminary  residual  resistance  measurements  indicate  the  sample  is 
not  pure  enough  for  the  primary  objective.  An  electron  beam  zone 
refiner  now  on  order  will  be  used  to  further  refine  the  palladium.  In 
the  meantime  the  possibility  of  using  pieces  of  the  sample  for  tunneling 
experiments  are  being  investigated. 
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1.1.6  Growth  of  GaAs  and  Ga  In.  As  Thin  Films 

x  i-x 

GK-698,  National  Science  Foundation 

DA-3  1  -  124-ARO-D-4 50 .  U.S.  Army  Research  Office 

L.  R.  Dawson,  C.  T.  Li  and  J.  M.  Whelan 

The  technique  has  been  extended  for  growing  thin  films  (l-10|j) 
of  GaAs  from  Ga  solutions  saturated  with  arsenic.  Previous  difficulties 
with  variable  surface  textures  have  been  substantially  reduced  by  protect¬ 
ing  the  substrate  during  the  reduction  of  the  residual  oxide  film  on  the 
liquid  Ga  surface.  Surface  defects  are  presently  limited  by  dislocation 
densities  in  the  substrates.  These  are  particularly  pronounced  for  the 
substrate  orientations  ,  (2 1  l )  and  (3  1 1 ).  Except  for  dislocations  ,  these 

orientations  are  favorable  for  thin  film  growth.  The  relative  advantages 

-2 

of  low  dislocation  substrates  (less  than  500  cm  )  oriented  in  the  [  1 1 1 J , 
[Ttt],  [110],  [  3 1 1  ]  ,  L311]and[21l]  directions  are  being  presently 
evaluated. 

Film  structures  containing  two  films  have  been  grown. 

Suitable  structures  for  bipolar  transistors  and  Gunn  oscillators  have 
been  grown.  Bipolar  transistors  wre  made  using  a  |  u  n  film  on  a 
semi-insulating  substrate  which  exhibited  pinch  off  at  0.6  V  bias.  An 
nt  10pn  film  structure  on  an  n+  substrate  was  operated  as  a  Gunn 
oac  'lator. 

Techniques  have  been  developed  for  doping  the  above  GaAs 
films  with  Te,  Sn  and  Ge.  Germanium  presently  seems  to  be  a  superior 
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acceptor  dopant  than  either  An  or  Cd.  Attempts  to  make  npn  o-  pnp 
transistor  structures  using  Zn  were  unsuccessful  because  of  the  Zn 
diffusion. 

Additional  facilities  for  growing  III  -  V  alloys  are  completed 
and  are  being  tested  with  GaAs. 

Evaluation  of  electrical  contacts  are  in  progress.  Electro¬ 
plated  Au-Sn  followed  by  alloying  looks  promising  as  a  reliable  n  type 
contact.  The  alloy  depth  is  approximately  1  pi . 


1.2  QUANTUM  ELECTRONICS  AND  LASERS 


1.2.1  Quantum  Theory  of  Noise 

DA-AR O-D-3  1  -  124- G990 ,  U.S.  Army  Research  Office 
W.  H.  Louisell,  J.  H.  Marburger,  D.  White 

1.  Lax  has  recently  given  a  "classical"  correspondence  for 

a  quantum  stochastic  process.  We  have  been  able  to  give  a  reformulation 

in  such  a  way  that  we  can  show  its  connection  with  a  classical  corres- 

2 

pondence  ntroduced  by  Gordon  .  In  addition  we  have  shown  why  Gordon's 
distrib  kion  function  obeys  a  "classical"  Fokker- Plane  k  equation  in 
which  the  diffusion  matrix  is  non-singular  and  why  his  P^  function  has  a 
diffusion  matrix  with  negative  eigenvalues. 

In  addition  we  have  obtained  a  very  simple  derivation  of  the 
quantum  regression  theorem  for  general  quantum  systems  in  which  multi¬ 
time  averages  can  be  computed  from  single  time  averages  for  Markoffian 
systems.  This  work  is  about  to  be  written  up  for  publication  with  J.  H. 
Ma  rburger . 

2.  Further  work  on  the  quantum  theory  of  the  laser  has  been 
done  with  Lax  and  is  about  ready  to  be  submitted  for  publication  as  QX!' 
in  the  series  of  quantum  stochastic  processes.  In  this  report,  after  we 
eliminate  the  dipole  moment  operator  adia!>atically  from  the  Langevin 
equations,  we  use  the  classical  correspondence  to  convert  to  a  classical 
problem  in  the  population  and  field  associated  viriables.  The  populations 
are  eliminated  classically  which  gives  a  further  check  on  prior  quantum 


eliminations  where  there  were  ordering  ambiguities.  We  express  the 
results  in  the  photon  number,  intensity  and  phase,  and  phase-number 
representations.  We  also  obtain  the  linexwdth  due  to  phase  diffusion  and 
steady  state  solutions  in  the  various  representations 

3.  We  a  e  also  studying  the  optical  parametric  oscillator 
in  which  the  atoms  are  treated  explicitly  rather  than  phenomenologically 
using  the  model  of  Graham  and  Haken  .  In  particular  we  are  studying 
critically  the  various  possible  methods  of  eliminating  the  atoms.  This 
work  is  being  done  with  D.  White  and  results  are  imminent. 
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1.2.2  Interaction  of  High  Intensity  Light  Beams  with  Matter 


A^-AFOSR-496-67,  Joint  Services  Electronics  Program 

W.  G.  Wagner 

Our  investigations  center  on  the  interactions  of  very  intense 
laser  beams  with  materials,  which  are  governed  by  noniinearities  pro¬ 
ducing  a  number  of  qualitatively  striking  phenomena.  Our  immediate 
objective  is  to  study  the  propagation  properties  of  intense  light  beams, 
since  only  when  these  are  known  in  some  detail  can  reliable  extraction 
of  material  nonlinear  coefficients  be  made  from  various  experiments. 
During  the  past  period  we  have  completed  our  study  of  the  initial  growth 
of  the  instabilities  in  the  propagation  of  electromagnetic  waves  of  very 
broad  extent  through  low  density  fluids,  and  we  have  begun  an  attempt  to 
calculate  the  distortions  in  beam  profiles  resulting  from  density  and 
thermal  variations  induced  in  the  medium  by  the  passage  of  the  optical 
beam. 

Our  previous  work  had  resulted  in  a  very  complicated  disper¬ 
sion  relation  which  enabled  one  to  calculate  the  temporal  or  spatial  growth 
rate  for  small  amplitude  disturbances  to  either  the  electromagnetic  beam, 
density,  or  temperature  of  th**  fluid.  The  theory  was  sufficiently  general 
co  include  the  couplings  of  thermal  deposition  by  absorption  of  radiation, 
electrostriction,  and  a  nonlinear  dependence  of  the  dielectric  response 
upon  the  electric  field  (Kerr  effect).  This  dispersion  relation  was  studied 
numerically  for  cases  of  interest,  and  it  was  found  that  the  correct 
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answers  for  the  maximum  instability  growth  rate  and  for  the  velocity  of 
propagation  of  the  disturbance  was  in  very  close  agreement  with  our 
previously  derived  analytical  estimates.  The  accuracy  of  the  agreement 
(about  1%)  is  somewhat  accidental  as  the  analytical  estimates  were  not 
expected  to  be  that  good. 

We  have  also  clarified  the  relation  of  our  work  to  that  of 
Brueckner  and  Jorna,  and  that  of  Livingston.  Our  results  are  in  sub¬ 
stantial  agreement  with  Brueckner  and  Jorna,  and  in  substantial  disagree¬ 
ment  with  Livingston. 

Further  work  on  this  problem  is  required  in  order  to  under¬ 
stand  the  nature  of  the  beam  distortion  and  possibly  the  saturated  state 
of  the  instability  development.  Our  only  hope  tor  the  answers  to  such 
questions  lies  in  an  advanced  digital  computer  study  of  the  problem. 
Toward  this  end,  a  program  for  the  solution  of  the  system  of  seven 
partial  differential  equations  in  three  independent  variables  is  presently 
under  development.  Substantial  effort  has  been  made  to  invent  a  rela¬ 
tively  sophisticated  iteration  technique  so  that  the  mesh  size  can  be 
increased,  but  we  will  have  to  prove  that  this  n^w  t,  ^hnique  is  a  stable 
one. 

During  this  past  period  we  have  also  completed  work  on  a 
few  projects  that  have  been  under  less  intensive  investigation  over  a 
period  of  eighteen  months.  The  first  of  these  is  the  investigation  of  the 
effect  of  saturation  on  self-trapped  optical  beams  in  liquids.  The  condi¬ 
tions  under  which  an  electromagnetic  beam  can  produce  its  own  dielectric 
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waveguide  and  propagate  without  spreading  have  been  discussed  for  a 
model  cf  the  nonlinear  dielectric  response  which  includes  the  effects  of 
saturation.  The  model  assumes  that  the  principal  contribution  to  the 
nonlinear  response  c-’ses  from  the  reorientation  of  anisotropic  molecules. 
For  the  lowest  order  'modes',  the  variation  of  the  f ;eid  amplitude  over 
the  cross  section  of  the  beam  has  been  derived  and  presented  graphically 
in  a  paper  published  in  the  IEEE  Journal  of  Quantum  Electronics,  It  has 
also  been  shown  that,  for  'cigar'-type  axially  symmetric  molecules, 
there  are  two  trapped  beams  possible  which  have  no  nodes,  but  which 
have  the  same  total  power. 

In  another  paper  published  during  this  period  in  the  IEEE 
Journal  of  Quantum  Electronics,  we  have  examined  the  equations  of  laser 
Q  switching  for  a  homogeneously  broadened  two-level  material.  It  was 
shown  that  inertial  effects  in  the  material  polarizability  may  play  an 
important  role  in  Q  switching  of  lasers  of  narrow  material  linewidth.  It 
was  further  shown  that  it  is  possible  to  amplify  pulses  of  durations  shorter 
than  the  inverse  linewidth  of  the  active  material,  provided  the  pulses  are 
energetic  enough  to  produce  appreciable  nonlinear  behavior  of  the  material. 
In  this  way  it  should  be  possible  to  produce,  from  a  relatively  long  pulse, 
pulses  of  short  duration  as  compared  to  the  inverse  linewidth  of  the  ampli¬ 
fying  material.  Numerical  examples  were  given  based  on  parameters 
appropriate  to  the  CO^  laser.  Even  though  the  model  of  a  homogeneously 
broadened  two-level  system  is  an  oversimplification  of  CO^  laser  opera¬ 
tion,  it  is  believed  that  guidelines  for  future  work  on  CO^  systems  can  be 
gleaned  from  these  results. 
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1.2.3  Optical  Experiments  with  Laser  Sources 

AF-AFOSR -496-67,  Joint  Services  Electronics  Program 
DAHC-04-69-C-0003,  U.  S.  Army 
W.  L.  Faust 

During  the  past  period  we  have  successfully  operated  the 

large  CO^  laser  contemplated  in  the  previous  report.  The  unit  delivers 

Q  switched  pulses  at  a  rate  of  400  sec  ,  on  various  vibrational- 

rotational  lines  of  CO^  selected  by  a  grating.  We  have  performed  a 

simple  (but  very  direct,  and  inexpensive)  calorimetric  measurement  of 

the  time-average  power,  making  use  of  the  known  heat  of  vaporization 

of  liquid  nitrogen.  This  average  power  is  about  one  watt.  We  estimate 

trie  peak  power  during  a  pulse  to  be  of  the  order  of  25  KW ,  but  we  have 

not  yet  equipped  ourselves  with  an  adequately  fast  detector  to  make  a 

reliable  measurement.  Unfortunately,  the  laser  tends  to  run  on  two 

rotational  lines  (or  more)  at  once;  the  grating  and  cavity  configuration  at 

present  do  not  afford  sufficient  discrimination,  even  with  the  use  of  an 

iris.  We  are  presently  contriving  to  double  the  effective  optical  length  of 

the  cavity,  which  should  solve  this  problem.  We  are  also  under  way  in 

the  construction  of  a  second  such  laser,  since  the  physical  experiments  we 

anticipate  (two  of  them,  quite  distinct)  involve  interaction  between  beams 

at  two  different  frequencies.  These  experiments  are  i)  high  intensity  (for 

double-quantum  effects)  two-beam  studies  of  the  local  mo  ’»  of  Ca  F  :  H 

2 

and  ii)  rotational  relaxation  studies  on  CCL,  itself.  On  another  front,  we 
have  inherited  a  Korad  ruby  laser  capable  of  about  25  KW  pulses  at  a  1 
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sec  rate.  We  hope  to  be  able  to  perform  Stokes  Raman,  scattering 
experiments  with  this  unit  as  a  source.  The  principal  problem  is  expected 
to  be  discrimination  against  intense  scattered  light  from  the  primary 
laser  beam,  which  (apart  from  likely  destruction  of  a  photo  mutliplier 
tube!)  could  easily  overwhelm  Raman  scattering  having  an  absolute 
intensity  which  would  otherwise  be  readily  detectable.  We  are  consi¬ 
dering  as  filters  GaA  P  alloy  films,  which  should  be  able  to  be 

s  x  ' 

1  -x 

fabricated  if  sufficient  interest  can  be  aroused.  If  upon  a  pure  GaP 

substrate  slab,  there  were  built  up  an  alloy  layer  >  20  |i  thick  on  each 

side,  with  0.30  <_  x  0.32,  the  result  should  be  ideal.  The  filter  should 

o 

be  essentially  transparent  for  \  >  7000  A  Stokes  Randan  scattering,  but 
should  exhibit  (very  strong)  direct  interband  electronic  absorption  at  the 
ruby  laser  wavelength.  Consultation  with  an  autority  (J.  Whelan,  of  the 
Materials  Science  Dept,  of  this  university)  leads  us  to  believe  that  fabri¬ 
cation  should  be  possible,  giving  collaboration  with  workers  having  the 
necessary  facilities.  Our  particular  scattering  problem  of  interest  is  in 
local  modes  of  impurities  in  GaP  (we  are  in  collaboration  with  W.  G. 
Spitzer  at  this  point),  but  it  is  clear  that  such  a  filter  would  be  a  boon  for 
any  sort  of  Stokes  Raman  scattering  study  with  a  ruby  laser  source. 

We  continue  expansion  of  other  general  facilities,  such  as  a 
vacuum  evaporation  deposition  system  for  making  our  own  metal /dielec trie 
mirrors,  special  gratings,  etc. 


1.2.4  F i r3_t  Order  Raman  Effect  in  Wurtzite  Crystals 
GP  7804,  National  Science  Foundation 
C.  A.  Arguello,  D.  L.  Rousseau,  and  S.  P.  S.  Porto 

First  order  Raman  scattering  from  four  crystals  (BeO,  ZnO, 

ZnS,  and  CdS)  with  the  Wurtzite  (C,,.)  structure  has  been  investigated  to 

o  V 

better  understand  the  effects  of  the  competition  between  the  long  range 
electrostatic  and  the  short  range  interatomic  forces,  and  to  determine 
the  contribution  of  the  elec tro- optic  coefficient  to  the  polar  phonon 
scattering  intensity. 

The  Raman  scattering  data  from  ZnO  and  CdS  are  in  substan- 

1  2 

tial  agreement  with  prior  investigations  on  these  crystals.  ’  In  the 

analysis  of  BeO,  both  the  transverse  and  longitudinal  components  of  the 

Aj  and  the  species  were  observed  and  both  E^  modes  were  found. 

The  BeO  frequencies  are:  AjLO-1081  cm  *;AjTO-678  cm  E^LO- 

1097  cm  ;  EjTQ-722  cm  ;  E^'684  cm  *  and  338  cm  \  For  all  of 

these  crystals  the  angular  dependence  of  the  transverse  and  longitudinal 

3 

modes  is  in  agreement  with  theory.  The  analysis  of  the  Raman  data 

4  5 

from  ZnS  is  in  disagreement  with  the  previous  Raman  investigations.  ’ 

In  the  present  study  it  is  shown  that  the  crystal  anisotropy  does  not  split 
the  transverse  or  longitudinal  phonons.  They  were  located  at  274  and 
352  cm  respectively.  The  E,,  phonons  were  demonstrated  to  occur  at 

2  74  cm  *  (accidentally  degenerate  with  the  transverse  phonon)  and  at 

.  _  -1  4  -  ] 

^5  cm  .  A  line  previously  assigned  as  an  E.,  phonon  at  68  cm  was 

shown  to  have  the  wrong  depolarization  characteristics  to  be  of  E, 
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symmetry. 


3  6 

Loudon  ’  has  shown  that  the  scattering  cross  section  of  a 
longitudinal  phonon  depends  on  a  delormaticn  potential  term  and  an  electro¬ 
optic  term  while  the  intensity  of  a  transveise  phonon  depends  only  on  a 
deformation  potential  term.  The  electro-optic  coefficient  may  be 
written  in  the  form 


Z  =  ±  (AS.  ±  BS  ) 
ij  l  f 

where  and  are  the  scattering  efficiencies  for  the  longitudinal  and 
transverse  phonons  respectively,  and  A  and  B  a,e  known  constants.  By 
measuring  intensities  of  the  transverse  and  longitudinal  phonons  the  three 
independent  electro-optic  coefficients  were  determined  for  each  of  the 
crystals  BeO,  ZnO,  and  CdS.  Comparison  of  these  values  with  those 
measured  by  other  methods  at  constant  strain^  demonstrates  that  very 
accurate  measurements  of  electro-optic  coefficients  may  be  made  by  tb 
Raman  effect. 
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1.2.5  Quantum  Electronic  Investigation  of  Cross-Relaxation 
in  Rare-Earth  Crystals 

NGR-05-018-044,  National  Aeronautics  and  Space 
Adminis  tration 

L.  G.  DeShazer,  M.  M.  Mann,  E.  A.  Maunders  and 
D.  K.  Rice 

Determination  of  relaxation  processes  in  laser  media  is 
currently  of  great  importance  in  understanding  interaction  of  radiation 
with  optically  pumped  media  We  are  exploring  correlations  between 
relaxation  effects  in  gain  saturation  of  laser  amplifiers,  time-resolved 
emission  spectra  of  laser  oscillators,  and  mode-locked  operation  of 
laser  oscillators.  In  particular,  we  are  examining  neodymium  crystals 
and  neodymium  glass  as  the  laser  media.  This  work  is  directed  towards 
development  of  a  model  for  the  relaxation  mechanism,  and  determination 
of  rate  equations  for  the  spectral  flux  and  spectral  inversion  density. 

The  model  will  be  used  to  derive  the  relaxation  term  for  the  rate  equation. 

Details  of  this  investigation  are  described  in  six  publications 
from  this  laboratory  and  are  briefly  summarized  below: 

1.  Employing  relaxation  effects  for  controlling  the  frequency  of 

1 , 2 

a  laser  oscillator. 
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Frequency  and  polarization  selection  of  neodymium  glass 
laser  oscillators  was  achieved  by  utilizing  light  from  a  controlled 
secondary  source.  When  light  from  the  secondary  source  was  injected 
into  the  laser  medium,  the  laser  emission  duplicated  the  spectral 
properties  of  this  light.  Thus,  a  narrowband  emission  was  obtained 
from  the  normally  broadband  laser  by  injecting  narrowband  light  at  a 
matching  wavelength.  Similarly,  a  polarized  emission  was  obtained 
from  the  normally  unpolarized  laser  by  using  polarized  light  from  the 
secondary  source.  The  rapid  energy  transfer  between  the  neodymium 
ions  allowed  the  laser  emission  to  duplicate  the  spectral  properties  of 
this  secondary  radiation. 

3,4,5 

2.  Spectral  broadening  of  rare-earth  ions  in  laser  media.  ’ 

This  program  is  studying  classical  experimental  procedures 
of  determining  line  broadening  parameters  of  rare-earth  ions  in  solid 
hosts.  Systematic  spectroscopic  studies  were  made  on  trivalent  neo¬ 
dymium  ions  in  various  laser  glasses.  Analysis  of  the  absorption  and 

o 

fluorescence  spectra  and  their  temperature  dependence  (from  4  K  to 
o 

300  K)  determined  the  energy  levels  of  neodymium.  The  complete  Stark 
splitting  of  the  low-lying  multiplets  showed  that  the  rare-earth  ions 
occupied  sites  uf  low  symmetry  in  the  glass.  The  magnitude  of  the 
splittings  correlate  closely  with  those  determined  from  oxide  crystals. 
The  inhomogeneity  of  the  glass  structure  is  manifested  in  the  large 
broadening  of  the  spectral  lines.  The  profile  of  the  component  lines  is 
Gaussian  with  halfwidths  from  40  to  85  cm  There  is  no  thermal 
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variation  of  lmewidths,  indicative  of  inhomogeneous  broadening. 

Previously  reported  fluorescence  linewidth  variation  with  temperature 
is  explained  by  the  thermal  variation  in  the  population  distribution  over 
the  Stark  components. 

A  study  of  europium  glass  was  made  to  supplement  and 
support  these  results  on  neodymium  glass.  The  spectrum  of  europium 
contained  many  simple  features  greatly  assisting  this  analysis.  Absorp¬ 
tion  and  fluorescence  from  singlet  levels  of  europium  made  it  possible  to 
study  the  linewiaths  without  the  complication  of  Stark  structure.  Observa¬ 
tion  of  J  -  0  -  -  J  -  0  transitions  in  europium  glass  limits  the  ion  site 
symmetry  to  a  low  symmetry  point  group,  regardless  of  the  mechanism 
invoked  for  allowing  the  transition. 

3.  Measurements  on  c  ross  -  relaxation  processes  between  ions 

6 

having  an  inhomogeneously  broadened  fluorescence  line. 

The  laser  injection  oscillator  experiments  demonstrated  their 
feasibility  in  studying  ionic  relaxation.  By  this  technique,  laser  oscilla¬ 
tions  were  forced  to  develop  from  the  narrowband  radiation  injected  into 
the  laser  cavity,  instead  of  broadband  optical  noise.  The  interionic 
transfer  allowed  the  feeding  of  this  narrowband  radiation  by  ions  having 
frequencies  outside  this  interval.  The  efficiency  of  energy  extraction 
within  the  narrow  bandwidth  measured  the  c  ross  -  relaxation  rate.  A  new 
experiment  is  being  constructed  utilizing  amplification  of  saturating 
giant-pulse  laser  radiation.  If  successful,  this  experiment  will  give. 
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for  the  first  time, all  the  line  broadening  parameters  of  a  rare-earth  ion 
in  a  solid  host. 
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1.2.6  Intensities  of  Crystal  Spectra  of  Rare-Earth  Ions 

AF-AFOSR-496-67 ,  Joint  Services  Electronics  Program 
M.  M.  Mann,  T.  S.  LaF ranee,  and  L.  G.  DeShazer 

In  this  program  the  mechanisms  which  give  rise  to  parity- 
ferbidden  transitions  in  the  spectra  of  rare-earth  ions  in  c  rystals  are 
being  investigated.  The  parity-forbidden  transitions  are  the  most 
commonly  occuring  transitions  in  rare-earth  ions,  and  all  presently 
observed  laser  action  in  rare-earth  crystals  is  due  to  these  transitions. 
The  present  theories  of  laser  ac  tion  in  rare-earth  c  .ystals  treat  transi¬ 
tion  probabilities  in  a  phenomenology':?!  manner.  In  order  to  develop  a 
more  fundamental  description,  it  is  necessary  to  have  a  suitable  mode! 
for  the  transition  probabilities.  Thus,  it  is  hoped  that  this  study  will 
lead  to  a  more  fundamental  understanding  of  laser  processes  in  rare-earth 
crystals.  Also,  this  understanding  will  hopefully  help  choose  the  "best" 
laser  host  for  rare-earth  ions,  a  very  illusive  host  at  present. 

1.  Enforced  Electric  Dipole  T  ransttions  of  Rarc-Furth  Ions. 

Doth  the  initial  and  final  states  associated  with  the  emission 
transitions  of  interest  belong  to  the  same  electronic  ontiguration.  Then, 
two  types  ot  radiation  m  crystal  spectra  must  be  considered:  a)  magnetic 
dipole  and  elect  ru  qiudrupule  which  also  occur  m  the  free  ion  spectra, 
and  b)  electric  dipole  which  does  not  occur  m  the  tree  ion  spin  tia  hut 
does  in  the  crystal  spectral.  Tins  parity -forbidden  electric  dipole 
radiation  is  called  enforced  electric  dipole  radiation  and  occurs  due  to 


-  (  . 


.  .c  interaction  of  the  rare-earth  ion  with  the  electric  field  of  the  crystal. 

1  2 

Judd  and  Ofelt  have  presented  theoretical  analyses  of  the  intensities 
associated  with  the  enforced  electric  dipole  transitions.  They  have  shown 
that  a  substantial  admixture  of  4f  and  5d  wavefunctions  occurs  when  the 
rare-earth  ion  is  not  located  at  an  inversion  center.  This  removes  the 
parity  restriction  and  allows  electric  ditole  transitions. 

We  have  studied  the  fluorescence  and  absorption  oscillator 
strengths  of  the  transitions  of  europium  in  LaCi^.  A  preliminary  fit  of 
the  Judd-Ofelt  model  to  this  data  was  made.  It  was  observed  that  in 
addition  to  the  intensity  measurements,  the  anisotropy  of  LaCl^  must  be 
known  to  permit  comparison  of  the  experimental  results  with  theoretical 
predictions.  Therefore,  a  theoretical  analysis  of  atomic  radiation  in 

3 

anisotropic  media  was  made  for  non- resonant,  linear  and  lossless  media. 
The  dei  ;ved  expressions  provide  an  analytical  rr  eans  of  relating  oscillator 
strengths  to  measured  intensity  data.  It  has  also  been  shown  that  the 
multiple  intensity  distribution  sensitively  reflects  both  the  macroscopic  and 
microscopic  dielectric  properties  of  the  host. 

No  general  expression  cai  be  given  for  the  local  field  param¬ 
eters.  A  suitable  analytical  model  must  be  chosen  for  each  case.  The  local 
field  parameters  may  be  more  profitably  viewed  as  experimentally  deter¬ 
mined  constants.  Thus  measurements  of  spatial  intensity  distribution 
can  be  used  to  obtain  information  on  the  local  environment  of  the  radiating 
ion. 
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2.  Ps cudo-Quadr upole  Transitions  of  Rare-Earth  Ions. 


It  was  observed  that  certain  fluorescence  lines  of  Eu:  LaCl^ 
were  unduly  intense  and  caused  a  peculiar  parameter  fit  to  the  Judd- 
Ofelt  model.  When  these  intense  transitions  were  neglected,  the  Judd- 
Ofelt  model  predicted  tne  remaining  intensity  data  in  a  reasonable  way. 
These  intense  transitions  were  named  "seudo-quadrupole"  because  they 
mimic  electric  quadrupole  selection  rules.  After  recognition  of  these 
transitions,  it  was  observed  that  almost  all  of  the  strong  laser  transitions 
from  rare-earths  appear  to  be  pseudo-quadrupole. 

4 

For  rare-earths  in  various  solutions,  Jorgensen  and  Judd 
noted  that  the  intensities  of  a  few  transitions  were  hypersensitive  to  the 
solvent,  while  the  intensities  of  the  other  transitions  were  insensitive. 

The  hypersensitive  transitions  mimicked  electric  quadrupole  selection 
rules.  Follov/ing  this  hint,  we  are  preparing  an  experiment  on  the 
europium  ion  in  europium  ethyisulfate  and  europium  magnesium  nitrate 
crystals  under  conditions  of  varying  pressure.  Large  single  crystals  wf 
europium  ethyisulfate  have  been  grown  and  are  being  prepared  for  this 
study. 
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1 . 2 u  7  Experimental  Studies  in  Nonlinear  Optics 

AF-AFQSR-496-67,  Joint  Services  Electronics  Program 
G.  L.  McAllister,  L,  Huff,  L.  G.  DeShazer 

This  experimental  program  is  investigating  the  nonlinear 
effects  in  the  propagation  of  light  in  organic  liquids.  The  main  emphasis 
in  this  program  is  on  self-focussing  effects  and  saturated  absorption 
effects.  Considerable  attention  was  given  to  the  development  of  a  high- 
power  narrow-pulse  laser.  Recent  studies  on  mode  classification  of  a 
giant-pulse  laser  are  being  submitted  for  publication.  The  investigation 
of  pulse  shaping  due  to  the  time-dependent  behavior  of  the  electric  field 
in  self-focussing  situations  has  been  started.  For  intensities  below  the 
threshold  of  seif-focusing,  such  as  occurs  at  the  beginning  and  end  of  the 
laser  pulse,  the  on-axis  intensity  of  a  beam  decreases  due  to  the  diffrac¬ 
tive  spreading.  If  this  beam  passes  through  a  nonlinear  medium,  such  as 
nitrobenzene,  the  spreading  effect  can  be  counteracted  by  self-focusing. 
For  intensities  above  the  self-focusing  threshold,  the  beam  is  squeezed 
and  the  net  result  is  suppression  of  the  leading  and  lagging  edges  and 
enhancement  of  the  pulse  center.  This  appears  as  a  time  compression 
of  the  pulse.  By  this  experiment,  much  information  can  be  obtained  on 
threshold  power  and  self-focusing  length. 

Experiments  on  saturated  transmission  of  ruby  laser 
radiation  by  organic  dye  solutions  were  continued.  Saturated  transmission 
is  a  reversal  ^f  the  well-known  observation  that  the  absorption  of  high 
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power  laser  radiation  by  many  organic  materials  reduces  with  increased 

laser  power.  The  sulfonatec1  forms  of  the  indanthrone  dye  become  opaque 

to  intense  radiation.  At  present,  four  sulfonated  forms  have  been  care- 

2 

fully  studied  for  powers  from  0.  1  to  100  megawatts  per  cm  . 

1.2.8  Computer  Studies  of  Self-Focusing  Optical  Beams 

DA-ARO-D-3  1  -  124-6920,  U.S.  Army  Research  Office 
J.  Marburger,  E.  Dawes,  R.  Reynolds,  L.  Huff, 

J.  Wagner,  J.  Reichert 

The  numerical  work  reported  in  the  previous  semiannual 

progress  report  has  been  completed  and  either  published  or  submitted 
1,2,3 

for  piolication  ’  ’  .  During  this  period  we  began  a  theoretical  investi¬ 
gation  of  the  effects  of  relaxation  of  the  nonlinear  refractive  index  on 
the  dynamical  self-focusing  of  an  intense  laser  beam.  We  have  found  on 
the  basis  of  simple  arguments  that  relaxation  effects  (which  cause  a 
delay  between  the  nonlinear  response  and  the  applied  field)  should  have 
the  same  qualitative  effect  on  the  focusing  beam  as  the  saturation  effects 
studied  during  the  previous  reporting  period.  Work  is  now  underway  to 
establish  this  connection  more  rigorously  and  use  it  to  estimate  the 
properties  of  small  scale  trapped  filaments  of  light. 

Additional  numerical  solutions  of  the  "steady  state"  self 

trapping  equations  have  been  obtained.  This  work,  which  augments  that 

4  5 

of  Reichert  and  Wagner  and  Gustafson  et  al  ,  was  motivated  by  the 
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discovery  of  a  technique  for  treating  the  nonlinear  local  field  theory  of 
Lorentz  correctly  to  all  orders  in  the  applied  field.  In  references  4  and 
5  the  local  fields  were  treated  only  in  a  rough  way.  We  have  found  that 
the  saturation  properties  of  the  nonlinear  refractive  index,  which  accord¬ 
ing  to  some  theories  determine  the  size  of  small  scale  trapped  light  fila¬ 
ments,  depend  sensitively  upon  the  local  field  corrections.  The  results 
of  this  investigation  are  reported  in  reference  6. 
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1.3  MAGNETISM 


1.3.1  Nuclear  Magnetic  Resonance  in  Ni  Rich  Ni-Cu  Alloys 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
S.  Ogawa 

The  magnetic  moment  of  the  Ni  atoms  originates  from  0.6 
holes  per  atom  in  the  3d  band.  Since  the  magnetization  of  Ni  decreases 
by  1  Bohr  magneton  when  one  Ni  atom  is  replaced  by  one  Cu  atom,  it 
becomes  zero  at  about  60  atomic  percent  of  Cu.  However,  the  screening 
effect  of  the  Cu  nuclear  charge  attacts  the  excess  electrons  to  fill  Ni  holes 
near  the  Cu  atom. 

The  NMR  method  was  used  to  investigate  this  localization  of 
charge  by  observing  hyperfine  field  spectra  of  both  Ni  and  Cu  nuclei  in 
alloys  up  to  10%  Cu  in  Ni  metal.  The  hyperfine  field  is  proportional  to 
the  spin  density  of  electrons.  By  means  of  a  computer  calculation  the 
spectra  were  analyzed  assuming  a  statistical  distribution  of  the  Cu  atom. 

One  Cu  atom  replacing  one  Ni  atom  as  a  first  or  second  neighbor  to  a  Ni 
atom  reduces  the  field  by  4.3  and  2.  5%  of  that  of  pure  Ni  metal,  respec¬ 
tively.  The  effect  of  the  same  replacement  as  a  first,  second  or  *hird  j 

j 

neighbor  of  a  Cu  atom  causes  the  reductions  of  8.  2,  5.  4  and  1 . 8%  of  the 
Cu  hyperfine  field,  respectively.  From  these  data  the  following  conclu¬ 
sions  can  be  drawn: 

1.  Ni  hyperfine  field  is  mainly  produced  by  the  core  polari¬ 
zation  and  the  contribution  of  neighboring  Ni  atoms  for  it 
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is  about  i/3.  The  Cu  hyperfine  field  is  produced  by  the 
conduction  electron  polarization  of  the  neighboring  Ni  atoms. 

2.  The  effect  of  Cu  atoms  in  Ni  metal  is  screened  out  by  the 
several  neighbor  shells. 

3.  The  average  Cu  hyperfine  field. calculated  with  the  effective 
factors  and  probability  numbers  of  Cu  atoms  of  the  successive 
shells,  decreases  with  Cu  concentrations  as  does  the  change 
of  the  saturation  magnetization. 

1.3.2  Ion  Configuration  in  Spinels 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
J.  Smit 

The  saturation  magnetic  momeu  and  other  technical 
properties  of  ferrites  are  largely  determined  by  the  distribution  of  the 
various  metal  ions  over  the  available  crystallographic  sites.  In  spinels 
and  garnets,  which  are  the  crystal  structures  of  the  most  important 
ferrites,  the  iron  group  metal  ions  occupy  either  octahedrally  or  tetra- 
hedrally  coordinated  sites,  the  nearest  neighbors  being  anions  like 
oxygen  or  sulfur.  The  magnetic  moments  of  the  metal  ions  on  these 
sites  are  anti-parallel. 

The  anions  are  asymmetrically  surrounded  by  the  metal  ions, 
so  that  an  electrical  field  at  the  anions  results,  which  gives  rise  to  a 
polarization  which  lowers  the  energy.  It  is  shown  chat  in  spinels  this 
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energy  should  always  favor  the  so-called  normal  configuration  in  which 

divalent  ions  occupy  tetrahedral  sites.  Exceptions  occur,  but  then  the 

5  10 

ions  or,  tetrahedral  sites  always  have  ad  or  d  electronic  configuration. 
This  condition  is  favorable  for  a  partly  covalent  bonding,  which  is  parti¬ 
cularly  stable  in  these  configurations. 

Recently  ferromagnetic  spinels  of  the  type  CuCr^S^  have  been 
reported  in  literature.  There  seems  to  be  strong  evidence  that  the  Cu 
ion  in  these  compounds  is  monovalent,  and  one  of  the  Cr  ions  quadruvalent. 
It  is  shown  that  these  configurations  are  quantitatively  consistent  with  the 
above  mechanisms,  as  they  are  stabilized  by  the  strong  polarization  of 
the  sulfur  ions. 

A  report  of  these  considerations  is  published  in  "Solid  State 
Communications  ". 

1.3.3  Mosabauer  Effect  Study  of  Mixed  Lithium-Zinc  Ferrites 
AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
J.  Young  and  J.  Sin  it 


Abstract 


Mossbauer  spectra  were  measured  at  room  temperature  for 

2  +  3+  |  1  + 

six  samples  of  the  lithium-zinc  ferrite  system,  (Z  Fe,  |  L 

x  0. 5-x/2 

3+  )0  .  Hyperfine  magnetic  fields  decrease  and  eventually  vanish 

F'l.5,*/2  4 

as  the  zinc  content  is  increased.  For  pure  lithium  ferrite,  the  fields 
are  515  Koe  for  the  tetrahedral  sites  and  520  Koe  for  th  octahedral  sites. 
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Pure  zinc  ferrite  is  paramagnetic  with  a  quadrapole  splitting  of  0.40 
mm/sec.  Intermediate  samples  show  a  gradual  transition  from  ferri- 
magnetic  to  paramagnetic  behavior. 

The  main  features  of  the  observed  spectra  can  be  accounted 
for  by  a  model  assuming  only  nearest  neighbor  inter-sublattice  exchange 
interactions  are  important. 

Experimental  Procedure 

The  mixed  lithium- zinc  ferrites  have  the  molecular  formula 

(Z  ^  Fe^  1  L.  i  1.  .  ,_)0 . .  Samples  having  x  =  0,  0.2,  0.4,  0.6,  0.8, 

n  1 -x  0.  5-x/2  4 

x 

and  1. 0  were  made  by  wet  mixing  ZnO,  Li^O,  and  Fe^O^  powders  in  the 

appropriate  stoichiometric  ratios,  drying,  pressing  into  pellets,  and 

firing  overnight  at  1050°C.  The  resulting  reddish-brown  pellets  were 

ground  to  a  fine  powder  and  set  in  sample  holders  with  paraffin.  Sample 

2 

thickness  was  adjusted  to  provided  about  30  mg  of  iron  per  cm  . 

Mossbauer  spectra  were  measured  using  a  system  manu¬ 
factured  by  the  Technical  Measurement  Corp.  ,  operating  in  its  analog 
sample  mode  and  later,  with  a  special  modification,  in  a  synchronized 
multiscaling  mode.  The  measured  linewidths  of  a  thin  natural  iron  foil 
are  about  0.  35  mm/sec.  Measurement  times  varied  from  sample  to 

sample  in  order  to  achieve  suitable  3  ignal- to- noise  ratios  with  spectra 

.  -  57 

of  different  complexities.  The  Mossbaucr  source  was  5  me  Co  in  Pd. 

Discussion  of  Results 

The  Mossbauer  spectra  of  all  the  samples  are  shown  in 
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Fig.  1.  Six-line  spectra  characteristic  of  magnetic  ordering  are 
observed  for  the  samples  having  x  -  0  through  0.6.  In  pure  lithium 
ferrite  (x  -  0),  the  lines  corresponding  to  tetrahedral  (A)  and  octahedral 
(13)  sites  are  not  clearly  resolved.  However,  line  widths  and  line  shapes 
vary,  and  the  spectrum  can  be  decomposed  into  two  sets  of  lines  having 
slightly  different  magnetic  fields  and  isomer  shifts.  Neither  set  shows 
any  appreciable  quadrupole  splitting. 

The  x  -  0.  2  and  0.4  samples  show  a  set  of  small  narrow  lines 
clearly  separated  from  a  set  of  large  broad  lines.  The  large  lines  have 
asymmetrical  line  shapes  with  steeper  slopes  on  their  high  velocity  sides. 
Approximate  measurements  of  the  ratios  of  the  areas  under  the  large  and 
small  lines  give  results  close  to  the  ratios  of  iron  atoms  on  octahedral 
and  tetrahedral  sites.  This  suggests  assigning  the  small  lines  to  the 
tetrahedral  sites.  However,  these  small  lines  also  show,  a  considerable 
quadrupole  splitting  and  in  general  are  similar,  though  not  identical,  to 
the  spectrum  of  a-Fe^O^.  Hence,  they  may  be  the  result  of  a  separate 
phase  of  Fe^O^  present  in  these  samples.  It  is  expected  that  x-ray 
diffraction  studies  will  resolve  this  question. 

In  the  x  -0.8  sample,  the  small  lines  have  disappeared  and 
the  large  ltnes  have  become  very  broad  and  shallow.  Also,  a  quadrupole 
split  central  doublet  hafe  appeared  indicating  that  some  iron  on  octahedral 
sites  has  become  paramagnetic.  The  x  0.8  and  1.0  samples  show  only 
the  central  doublet.  Apparently  all  the  iron  in  these  materials  has 
become  paramagnetic. 


A  summary  of  the  hyperfine  magnetic  fields,  quadrupole 
splittings  and  isomer  shifts  of  all  samples  is  given  m  Table  I.  Isomer 
shifts  are  measured  relative  to  the  center  ot  symmetry  of  natural  iron. 


Table  I 
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X 

0 

0.  2 
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0.  6 
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1 
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1  i 

.  20 
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Theoretical  Calculations  of  the  Spectra 

We  have  attempted  to  explain  tv.  o  main  features  of  the  observed 

spectra,  the  decrease  in  the  hyperfine  fielu  experienced  by  B  sites  and 

the  broadening  and  increasingly  asymmetrical  shape  of  these  lines  as 

zinc  content  is  increased.  The  model  used  assumes  i..at  the  hyperfine 
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field  at  the  nucleus  ot  an  Fe  atom  is  proportional  m  its  time  average 
spin,  and  that  the  spin  is  determined  by  inter- sublattic  e  exchange  inter¬ 
actions  bc'ween  an  iron  atom  and  its  nearest  neighbor  iron  atoms  on  the 

opposite  sublattice.  In  this  approximation,  the  value  of  the  inter  -  subla  ttu  e 

A  [i 

exchange  integral.  J  .  is  -  kT  i  42.8  for  pure  lith'um  ter  rite  which  has 

c 

T  -  Q70  K.  We  assume  that  the  value  ot  J  is  unchanged  bv  the  addition 

c 

of  zinc.  Thus,  the  effort  of  the  zinc  is  assumed  to  be  due  only  to  a  >  hange 
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in  the  number  of  nearest  neighbor  Fe's  seen  by  atoms  on  the  two 
sublattices. 

In  the  spinel  lattice,  each  A  site  "sees"  12nn  B  sites, and 

each  B  site  "sees"  6  nn  A  sites.  On  the  average  an  A  site  "sees" 

n„(x)  =  9  +  3x  nn  Fe's  on  B  sites,  and  a  B  site  "sees"  n  (x)  =  6-6x  nn 
B  A 

Fe's  on  A  sites. 


The  time  average  spins  of  Fe  atoms  on  A  and  B  sites  are  then 


given  by: 


and 


2  B 5/ 2  ^  5n 


£^-sB 

B  kT  z 


) 


,B 


2  B  5/ 2 


C  5nA  T? 


AB 


,A  ^ 
kT  "*z  J 


o  JAB 

At  T  =  300  K,  -  =  -0.0755.  If,  as  a  first  approximation,  we  take 
SA  =  ~  and  SB  =  -  -j  within  the  arguments  of  the  Brillouin  functions,  we 

2  u  Z  far 

find  SA  -  ~  B  ,  (0.  944  n  )  and  SB  =  -  ~  B  .  (0.  944  n  ).  It  appears 

that  the  proportionality  constant  relating  the  spin  may  be  different  for 

A  i  A  13  B 

the  two  sites.  Thus,  H.  ,  =  a.  I  S  and  H.  ,  a  I  S  I.  These  constants 

hf  A  z  hf  B  z 

can  be  determined  from  the  measured  fields  of  pure  lithium  ferrite.  We 
2  2 

find  a  =  —  *  52  5  Koe  and  a  -  —  •  542  Koe.  Taking  n  and  n  charac- 

A  b  B  b  A  B 

teristic  of  each  sample,  we  can  calculate  the  expected  hyperfine  fields 
for  the  A  and  B  s  ®s.  Table  II  gives  the  calculated  and  measured  field. 
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Table  II 


Sample  Calculated  Field  Measured  Field 


T  . r 

A 

B 

A  ^ 

B 

X  -  0 

51  5  Koe 

520  Koe 

51  5  Koe 

520  Koe 

i 

0,  L 

t"‘  ■  + 

;  520 

498 

527 

486 

1 

i 

i 

"1 . . 

i 

0.4 

[  0. 6 

!  520 

466 

524 

462 

| 

i  520 

3  96 

_ 

404 

! 

J 
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1.4  DEFECTS  IN  CRYSTALS 


1.4.1  Defect  Chemistry  of  CdS 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 


1.4.  i.  1  High  Temperature  Hall  Effect  Measurements 
F.  A.  Kroger  and  G.  Herechman 


Apparatus  tor  high  temperature  Hall  effect  measurements  of 
CdS  in  Cd  vapor  has  been  modified  to  make  possible  the  use  of  graphite 
yarn  as  current  lines.  Measurements  on  three  different  crystals  with 
foreign  donors  (Al)  in  the  concentration  range  5x  10^-  iO^cm  ^  all 
show  the  following  features: 


1)  in  the  range  500  -  700  C,  the  concentration  of  electrons 
«  ,  indicating  precipitation  -  dissolution  equilibrium 

for  the  donor  (Al^S^), 


2)  between  700  -  900°C,  the  concentration  of  electrons  is 

constant  at  low  p  ,  but  increases  with  p  at  higher 
Cd  Cd 

cadmium  pressures.  This  indicates  that  Cd  excess  gives 
.rise  to  donors  (Cd  or  V_).  It  has  so  far  not  been  possible  to 
decide  between  single  and  double  ionization  of  the  native 
donors. 
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1.4.  1.2  Cd  Self  Diffusion 


F.  A.  Kroger  and  V.  Kumar 

Measurements  of  the  self  diffusion  of  cadmium  in  CdS  under 
Wull-defineu  conditions  (P  T)  using  radio  active  Cd  have  been  started. 

The  goal  of  this  project  is  the  determination  of  the  species 
responsible  for  the  transport  of  cadmium  and  its  thermodynamic  forma¬ 
tion  parameters.  Diffusion  constants  found  so  far  in  the  purest  CdS 

,  ,  ,  ,„17  -3  ,  ,  ,  -10  ,  -11  2 

available  (containing  «  10  cm  Al)  are  of  the  order  10  -  10  cm 

sec  * ,  2-3  orders  of  magnitude  smaller  than  those  reported  by  Woodbury 

in  a  similar  study. 


1.4.2  Electrochemistry  of  Solids 


AF-AFOSR-68- 1405,  Air  Force  Office  of  Scientific  Research 


1.4.2.  1  Zirconia  as  an  Oxygen  Pump 
F.  A.  Kroger  and  D.  Yuan 

With  a  cell  consisting  of  a  zirconia  tube  with  inside  and 
outside  f'~*inum  electrodes,  oxygen  can  be  transferred  from  inside  to 

outside  by  passing  a  current  (outside  positve).  Stai.  .vith  nitrogen 

-  5  3 

containing  «  2  x  10  atm.  of  O-,  flowing  at  a  rate  of  300  cm  per 

-36 

minute  through  the  tube,  oxygen  pressures  as  low  as  10  atm  may 
be  reached.  This  low  pressure  limit  is  determined  by  the  onset  of 
electronic  conduction  in  the  zirconia. 
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1.4.  2,  2  The  Conduction  Properties  of  ZrO^  and  ThO^ 

F.  A.  Kroger  and  H.  Yanagida 

The  point  where  electronic  conduction  sets  in  differs  from 
material  to  material.  We  are  investigating  which  factors  are  responsi¬ 
ble  for  these  differences.  The  point  where  the  electronic  transference 
is  1/2  can  be  found  by  studying  the  cell 

Al,  A1203  |  ZrOz|  02  , 

Th02 

and  comparing  the  emf  with  that  corresponding  to  the  free  enthalpy  of  the 
reaction  2A1  +  ©2  -•  A^O^.  In  order  to  be  able  to  make  dense  materials, 

it  was  necessary  to  study  the  sintering  behavior  of  Zr02  and  Th02>  Gas 
tight  disks  can  now  be  made. 

1.4.2.  3  Sulfur  Electrodes 

F.  A.  Kroger  and  D.  Yuan 

Three  types  of  electrochemical  cells  sensitive  to  S  vapor 
are  being  investigated 

1)  PtOz|  Zr02l  C,S{g) 

This  cell  responds  to  the  sulfur  pressure,  the  empf  ncreasing  with  p  . 

o 

This  suggests  that  the  zirconia  acts  as  an  oxygen  electrode  in  the  usual 
way,  a  low  oxygen  activity  being  caused  by  the  reaction 

S(g)  *•  o2  -  so2 
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2)  Na  |  Na  glass  |  C , S 

3)  Na  |  Na2<J.  llAl^  |  C,S 


1.4.  2. 4  Polarization  Studies  of  AgCl 

F.  A.  Kroger  and  Y.  van  der  Meulen 


cells 


Measurements  have  been  carried  out  of  current  as  f(E)  for 


Ag  j  AgCl  |C  and  C,  Cl^  |  AgCl  |  C  or  Pt 


The  former  show  abnormally  large  noise  effects.  In  addition,  they  have 
currents  «  exp  (Ee/akT)  with  a  >  1.  Probe  measurements  are  now 
being  made  to  see  where  the  noise  is  generated  and  to  find  out  whether 
surface  barriers  play  a  part. 

The  Cl2  cells  show  more  normal  behavior;  however,  a  cannot 
be  determined  because  of  the  presence  of  an  emf  at  zero  current.  A 
new  type  of  cell  is  being  made  in  which  the  Pt  electrode  is  inside  the 
crystal;  this  should  eliminate  the  emf. 


1 . 4. 2.  5  Th Conduction  Mechanism  of  NH  .Cl 

4 

F.  A.  Kroger 

In  order  to  complement  the  emf  measurements  carried  out 
previously,  the  conductivity  of  NH^Cl  as  f  (temperature)  was  measured 
in  various  gas  atmospheres  N2  +  x  NH3+  yH^.  No  significant  variation 
of  the  conduction  with  the  gas  composition  was  observed.  This  indicates 


-85- 


that  the  emf's  are  established  by  equilibration  of  a  relatively  thin  surface 
layer. 

Papers  Submitted  for  Publication 

1.  P.  B.  P.  Phipps  and  F.  A.  Kroger,  J.  Phys.  Chem.  Solids.  Long- 
Wave  Length  Absorption  in  Brominated  AgBr. 

2.  F.  A.  Kroger,  J.  Phys.  Chem.  Solids,  Defects  and  Phase  Stability 
of  Transition  Metal  Compounds. 

3.  H.  Yanagida  and  F.  A.  Kroger,  J.  Am.  Ceram.  Soc.  ,  The  System 
A1  -  O. 

4.  D.  Yuan  and  F.  A.  Kroger,  J.  Electrochem.  Soc.,  Stabilized 
Zirconia  as  an  Oxygen  Pump. 

1.4.3  Defect  Chemistry  of  A^Oj 

AF-AFOSR-68-1405,  Air  Force  Office  of  Scientific 
Research 

F.  A.  Kroger,  R.  J.  Brook,  J.  Yee 

1.  In  a  program  investigating  the  electrical  properties  and 
defect  structure  of  Al^O^,  the  emf  of  the  cell 

Pt(02)|  Al203|  Pt(O^) 

is  being  measured  on  crystals  grown  from  a  PbF^  flux  and  doped  with 
MgC  and  SiO?.  Ionic  conduction  is  indicated  for  undopcd  A1  O  up  to  the 

•  L*  J 

highest  temperature  used  (1400°C). 

2.  The  conductivity  of  the  crystals  is  being  measured  using 
a  volume  guard  ring  to  eliminte  leakage  currents.  A  constant  activation 
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energy  of  2.5  eV  is  observed  over  the  range  900°  to  1400°C. 

3.  An  attempt  is  being  made  to  observe  the  mass  transfer 

+++ 

expected  during  conduction  by  A1  migration  with  metallographic 
techniques;  the  usual  weight  change  method  is  susceptible  to  errors 
introduced  by  surface  conduction  paths. 

4.  Measurements  of  the  conduction  parameters  at  very  low 
oxygen  pressures  are  now  being  prepared  using  the  zirconia  oxygen  pump 
described  above  (section  1 . 4.  2.  1 ). 

The  combination  of  these  approaches  should  afford  a  descrip¬ 
tion  of  the  defect  structure  of  the  oxide. 

1.4.4  Grain  Growth  in  Ceramics 

GK  1487,  National  Science  Foundation 

R.  J.  Brook 

Grain  growth  phenomena  are  of  importance  both  in  the 
processing  and  the  properties  of  ceramic  materials,  and  the  present 
program  is  an  attempt  to  define  the  parameters  that  determine  the  type 
and  the  kinetics  of  grain  growth  for  a  given  system. 

Emphasis  has  been  placed  on  the  effect  of  two  variables, 
namely,  impurity  concentration  and  microstructure. 

Theoretical  aspects  of  impur  ty- influenc  ed  normal  grain 
growth  were  described  in  the  previous  report  ,  present  work  is  con¬ 
cerned  with  defining  the  role  of  porosity  in  the  growth  process  under 
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I 

I 

conditions  of  different  grain  size,  pore  size,  *nd  pore  distribution.  •  , 

Experimental  verification  for  the  impurity  effects  is  being 
sought  in  a  study  of  impurity  distribution  at  different  grain  sizes  in 
MgO  by  means  of  electron  microscopy  and  electron  microprobe  analysis. 

Verification  for  the  microstructural  effects  is  being  studied 
by  comparing  observed  growth  kinetics  in  very  small  grain  size  (thin 
film)  ZnO  with  those  predicted  by  the  analysis. 

Reference 

1.  R.  J.  Brook,  "The  Impurity  Drag  Effect  and  Grain  Growth  Kinetics", 

Scripta  Met.  (2)  375  (1968). 
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1.  5  METALS 


1.5.1  Experimental  Studies  of  Fermi  Surface  Topology  in  Metals 
AF-AFOSR -496 - b7 ,  Joint  Services  Electronics  Program 
M.  H.  Halloran,  J.  Brewer  and  R.  Parker 

The  experimental  determination  of  Fermi  surface  topologies 
in  metals  has  been  an  active  field  of  research  in  solid  state  physics  for 
several  years.  The  importance  of  this  research  derives  from  the  very- 
precise  information  which  can  be  learned  about  the  electronic  band 
structure  at  the  Fermi  energy.  Comparison  of  this  information  with 
theoretical  energy  band  calculations  provides  the  most  accurate  check 
lor  these  calculations  at  the  present  time.  In  addition,  these  studies 
yield  accurate  determinations  of  electron  effective  masses,  which  serve 
to  provide  information  concerning  the  interaction  of  the  electrons  with 
phonons  and/or  magnons  via  the  resulting  mass  enhancement  as  well  as 
additional  further  information  for  checking  the  band  calculations.  Some 
indirect  measurement  of  electron  mean  collision  time  is  also  obtainable, 
as  well  as  carrier  c one ent  rat' ons  ,  magnetic  breakdown  probabilities  and 
open  orbit  directions.  Our  research  is  specifically  directed  to  three 
areas:  Fermi  surface  studies  m  intermctall'r  compounds,  B-H  effects 
in  noble  metals,  and  Fermi  surface  studies  in  transition  metals  and 
rare-earth  metals. 

Fertni  Surface  S tv  lies  in  Intermetallic  Compounds 

Data  are  being  obtained  presently  on  the  intermetallic 
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compound  AuGa,  using  the  deHaas  -  van  Alphen  Effect  as  measured  with 

Cm 

a  sensitive  torsion  balance  in  fields  up  to  28  kOe.  Magnetothe rmal 
oscillations  will  also  be  used  and  a  direct  comparison  of  the  relative 
sensitivities  of  the  two  techniques  should  also  be  possible.  A  crystal 
of  Auln  has  been  zone  refined  and  will  be  checked  soon  to  determine  if 
it  is  sufficiently  pure  to  permit  deHaas-van  Alphen  measurements.  AuGa^, 
Auln  and  AuAF  form  a  family  of  compounds  which  possess  many  simi- 

Cm  Cm 

larities,  but  AuGa^  has  exhibited  an  anomalously  large  Knight  shift  and 

it  has  been  suggested  that  its  Fermi  surface  may  be  substantially  different 

(12) 

from  the  other  two  compounds.  Early  results v  '  '  have  not  been  conclu¬ 
sive,  but  only  a  small  portion  of  *he  Fermi  surface  is  known  thus  far. 
Studies  en  a>l  three  compounds  are  planned,  and  if  low  field  results  remain 
inconclusive,  additional  studies  in  a  .OOkOe  solenoid  will  be  made. 

B-H  Effects  in  the  Noble  Metals 

The  measurement  of  non-linear  effects  in  the  deHaas-van 

A’phen  effect  due  to  differences  between  B  and  H  inside  the  sample  have 

been  extensively  investigated  only  in  beryllium^  \  They  liave  also  been 

(4,  5* 

observed  in  the  noble  metals  (Au.  Ag  and  Cu)  and  perhaps  other 

material^.  The  effect  is  manifested  in  the  line  shape  of  the  oscillations, 
the  appearance  of  large  amplitude  difference  frequencies,  and  anomalous 
oscillation  amplitude  dependence  on  temperature.  A  theoretical  explana¬ 
tion  was  successful  in  Be.  but  certain  qualitative  differences  assoc  iated 
with  the  Fermi  surface  topology  complicate  the  interpretation  in  the 
noble  metals,  and  the  detailed  Fermi  surface  measurements  with  properly 
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shaped  samples  required  for  a  complete  explanation  have  not  been  made 

(6) 

as  yet.  Recent  nm  r  measurements  in  silver  have  indicated  the  exis¬ 
tence  of  domains  due  to  the  B-H  effects,  consistent  with  the  earlier  inter  ■ 
pretations  in  beryllium.  .*•  good  sample  of  silver  is  on  hand  and  when 
higher  fields  are  available  in  the  near  future  we  will  undertake  the 
required  deKaas-van  Alphen  measurements. 

Fermi  Surface  Studies  in  Transitio  .  Metals  and  Rare  Earth  Metals 

These  materials  are  racterized  by  d-bands  or  f-bands  at 
the  Fermi  energy,  resulting  in  large  electronic  masses,  frequently 
with  large  mass  enhancement  due  to  phonon  and/or  magnon  interactions 
with  the  electrons.  Therefore  very  high  fields  and/or  low  temperatures 
are  required  to  study  the  Fermi  surfaces.  To  provide  the  former,  a  100 
kOe  superconducting  soienoid  was  delivered  last  January  but  failed  to 
meet  specifications  and  was  returned.  Subsequent  testing  by  ihe  manu¬ 
facturer  resulted  in  further  performance  degradation,  and  a  second  cod 
was  also  unsuccessful.  A  new  coil  from  a  different  manufacturer  should 
be  delivered  before  January  1,  enabling  research  in  this  area  to  commence. 

In  the  meantime,  in  a  joint  venture  with  Bell  Telephone  Laboratories  and 
the  use  of  their  tacilities,  Fermi  surface  measurements  hav  e  been  made'  ' 
on  Nb  and  Ta  which  indicate  that  the  Fermi  surface  model  determined 
from  APW  band  calculations  and  the  mass  enhancement  is  in  agree¬ 
ment  with  estimates  based  on  superconducting  transition  temperatures. 

References 

1.  M.  H.  Halloran  and  J.  H.  Wernick,  Bull.  'm.  Phy.s.  Sot.  10.  450  (1965). 
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4.  A.  o.  Joseph  and  A.  C.  Thorsen,  Phys,  Rev.  138,  A1159  (1965); 
Phys.  Rev.  HO,  A2046  (1965). 

5.  M.  H.  Halioran,  F.  S.  I.  Hsu  and  J.  E.  Kunr.ler,  in  Proceedings 
of  the  Tenth  International  Conference  on  Low  Temperature  Physics , 
Moscow  (1966)  -  Vimiti  Publ: -hing  House,  Moscow. 
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1 .  5.  2  Basic  Structure  -  Property  Correlations  in  Shock-Loadea 
Metals  and  Alloys  by  Transmission  Electron  Microscopy 
AF-AFQSR-496-67,  Joint  Services  Electronics  Program 
L.  E.  Murr  and  J.  V.  Foltz 

Thia  research  has  as  its  general  purpose  the  study  of  the 
relationship  of  residual  mechanical  properties  to  explosive-shock  induced 
defect  structures  in  metals  and  alloys  of  engineering  interest.  The 
general  plan  of  attack  follows  the  approach  used  in  previous  research, 
portions  of  which  were  supported  in  part  by  this  grant1  \ 

During  the  period  covered  in  this  report,  assemblies  of 
shock-loaded  Inconel  600  alloy  (76%  Ni,  16%  Cr,  7%Fe,  balance  impuri¬ 
ties)  foil  0.001  in.  thic.,.  sandwiched  between  pure  nickel  foils  0.004  in. 

* 

Research  Physicist,  Ballistics  Division,  U.S.  Naval  Weapons  Labora¬ 
tory,  Dahlgren,  Virginia;  on  temporary  duty  assignment  as  a  research 
associate  in  the  Department  of  Materials  Science,  U.S.C. 
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thick, as  shown  in  Fig.  i,  were  disassembled;  and  the  foils  elec tropolished 
for  direct  observation  of  the  residual  microstructures  by  transmission 


electron  microscopy.  The  unique  feature  of  Inconel  alloy  in  the  initial- 
annealed  state  is  the  presence  of  coherent  Ni^(Ti,Al)  precipitates  in  the 
matrix.  It  was  thus  of  fundamental  interest  to  observe  the  influence  of 
such  inclusions  on  the  propagation  of  dislocations  during  the  shock-front 
passage,  and  their  possible  influence  on  the  residual  mechanical  properties. 

As  an  obvious  approach  to  elucidating  the  residual  mechanical 
properties  of  shock-loaded  Inconel,  Vickers  microhardness  mea  arements 
were  made  on  the  recovered  sheet  specimen  using  a  50  gm  load.  The 
results,  when  compared  with  pure  nickel  as  shown  in  Fig.  2,  illustrate 
quite  strikingly  the  fact  that  in  the  pressure  range  0  -  400  Kb,  the 
strengthening  of  the  Inconel  matrix  continues  while  for  nickel,  hardening 
is  effectively  complete  after  200  Kilobars  pressure. 

The  reason  for  the  behavior  shown  in  Fig.  2  is  vested  in  the 
residual  defect  structure;  amenable  to  direct  observation.  Figure  3 
illustrates  the  initial  Inconel  substreture  (before  shock  deformation), 
and  the  typical  appearance  of  the  residual  microstructure  after  shock 
loading  to  100  Kb  pressure.  The  hardening  at  100  Kb  as  observed  in 
Fig.  2  can  now  be  observed  to  be  due  to  dense  dislocation  structures;  and 
possibly  the  interaction  of  dislocations  with  the  precipitates  as  shown  in 
Fig.  3(b).  It  is  also  of  interest  to  point  out  in  Fig.  3(b)  that  the  dense 
arrays  of  dislocations  surrounding  the  precipitates  were  presumably 
generated  at  the  particle-matrix  interface.  This  could  only  occur  if  the 
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Fig.  2.  Residual  microhardness  of  thin  nickel  and  Inconel  alloy  sheets 
following  explosive  shock  deformation. 
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Fig.  3.  (a)  Undeformed  (annealed)  Inconel.  Arrow  indicates  especially  prominent 
coherent  strain  contours  at  a  precipitate.  Other  precipitates  also  display  the  line 
of  no  contrast,  normal  to  the  diffraction  vectors  indicated  (b)  Dislocation  substruc 
ture  in  Inconel  shock-loaded  at  100  Kb.  Note  lack  of  coherency  strains  at  the  preci¬ 
pitates  and  the  interfacial  dislocations.  The  surface  orientation  in  (a)  and  (b)  is 
(110).  (125  KV) 


precipitates  became  strained  from  the  matrix,  i.e.  ,  they  became 
incoherent.  Evidence  for  this  phenomenon  is  readily  apparent  from  the 
fact  that  coherency  strain  contrast  evident  in  Fig.  3(a)  (especially  promi¬ 
nent  at  the  arrow)  is  not  observed  at  the  precipitates  in  Fig.  3(b)  following 
shock  deformation.  In  addition,  the  general  lack  of  prismatic  dislocation 
loops,  and  the  high  density  of  elongated  and  bowed  dislocations;  suggests 
that  the  precipitates  in  many  cases  were  sheared  by  the  passage  of  dis¬ 
locations,  or  acted  as  barriers  to  their  motion. 

Figure  4(a)  shows  the  occurrance  of  twin-faults  in  the  Inconel 
matrix  as  the  mode  of  deformation  at  3  70  Kb  pressure.  Figure  4(b)  shows 
for  comparison  the  corresponding  defect  substructure  in  pure  nickel  at  the 
same  pres-ure.  We  now  observe  that, with  reference  to  Fig.  2,  the 
marked  contrast  in  shock  hardening  between  pure  nickel  and  Inconel 
alloy  at  higher  pressures  is  due  to  the  difference  in  the  residual  defect 
structure. 

In  the  case  of  pure  nickel,  increasing  shock  pressures  simply 
decrease  the  dislocation  cell  sizes,  and  this  character  of  defect  persists. 
On  the  other  hand,  the  planar  dislocation  arrays  in  Inconel  (Fig.  3(b)) 
increase  in  density,  and  transform  to  twin-faults  in  the  pressure  ranee 
above  about  200  Kb. 

Future  work  will  concentrate  on  elucidating  the  mechanism 
of  precipitate  strengthening  in  the  shock-loaded  Inconel,  and  the  effects 
of  the  lower  stacking-fault  energy  of  the  Inconel  as  compared  with  pure 
nickel.  In  addition,  work  is  continuing  on  the  fabrication  of  very  high 
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Fig.  4.  (a)  Deformation  microtwin  structure  in  Inconel  foils  following  explosive 

shock-loading  at  370  Kb.  Note  the  region  of  no  contrast  where  the  twins  intersect. 
The  twin  width  can  be  measured  directly  from  the  width  of  this  region.  The  surface 
orientation  is  close  to  (100).  (b)  Dislocation  cell  structure  in  pure  nickel  foils  fol¬ 
lowing  shock  deformation  at  370  Kb  pressure.  The  grain  surface  orientation  is 
close  to  (110).  (125KV) 
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pressure  shock  assemblies  of  304  stainless  steel.  Three  assemblies  are 

intended  for  explosive  loading  at  pressures  above  1  megabar  and  at  low 

temperatures  in  order  to  investigate  the  shock- recovery  phenomenon 

2 

previously  observed  . 
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1.5.3  Analysis  of  Grain  Boundary  Equilibrium  and  Relative 

Interfacial  Free  Energy  Phenomena  in  Thin  Solid  Films 
AF-AFQSR-496-67,  Joint  Services  Electronics  Program 
L.  E.  Murr  and  R.  J.  Horylcv 

This  research  has  as  its  objective  the  investigation  of  inter¬ 
facial  free  energy  phenomena  and  equilibration  processes  in  solids.  The 
mode  of  attack  is  the  direct  observation  of  interfacial  structure  and 

geometry  in  representative  thin  sections  by  transmission  electron 
1 

microscopy  . 

During  the  current  research  per iod,  preliminary  studies 
concerning  the  characterization  of  low  angle  (energy)  and  high  angle 
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(energy)  grain  boundaries  in  polycrystalline  metals  and  alloys  were 
conducted  on  304  stainless  steel.  The  method  of  analysis  involves  the 
observation  of  grain  boundaries  in  the  electron  microscope  and  the 
recording  of  the  selected  area  electron  diffraction  pattern  in  the  adjoining 
grains.  After  the  assurance  of  the  accuracy  of  the  grain  orientations  on 
opposing  sides  of  the  interface,  the  diffraction  patterns  are  superposed, 
and  the  smallest  angular  difference  between  the  [  110]  directions  of 
grain  orientation  measured.  The  image  contrast  features  are  then 
employed  in  deciding  the  approximate  nature  of  the  interface  (lov  or  high 
energy);  and  the  angles  of  misorientation  are  then  plotted  in  corresponding 
histograms. 

The  preliminary  results  of  this  investigation  have  been  very 
encouraging,  since  there  is  a  marked  difference  in  the  angles  of  misorien¬ 
tation  for  the  two  boundary  types.  Figure  1  illustrates  the  preliminary 
data.  The  real  significance  of  these  results  is  that  for  the  first  time, 
grain  boundary  misorientations  have  been  measured  in  a  real  solid  as 
opposed  to  some  fabrication  of  oriented  single  crystals. 

Encouraged  by  the  preliminary  work,  we  are  now  pursuing 

the  detailed  analysis  of  crystal  misorientation  as  it  relates  to  grain 

boundary  structure  and  energy.  It  now  appears  possible  to  determine 

the  average  values  of  low-angle  grain  boundary  free  energy.  In  order 

to  pursue  this  point  further,  we  are  re-investigating  the  relative  twin 

boundary  grain  boundary  free  energy  ratio  in  stainless  steel  as  previously 
2 

measured  by  Murr  .  However,  in  the  present  extension  of  this  work, 
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(o)  LOW-ANGLE  GRAIN  BOUNDARIES  (b)  HIGH-ANGLE  GRAIN  BOUNDARIES 


the  grain  orientations  in  the  twinned  matrix  ar.d  the  opposing  grain  will 

be  recorded  and  the  misor w  "tat .on  of  [llC]  determ'  ^ed  for  each  junction 

system.  In  this  way,  we  hope  to  relate  the  average  grain  boundary  free 

energy  to  the  angular  misorientation  of  the  competing  grains  assuming  a 

constant  twin  boundary  free  energy,  6  ,  of  18  ergs/cm^  as  determined 

tb 

.  2 

previously  . 
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1 .  5.  4  Vapor  Deposition  of  Metals 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
L.  E.  Murr  and  R.  J.  Horvlev 

We  have  continued  the  fabrication  of  ultra -high  vacuum 
systems  for  evaporation  of  uncontaminated  thin  foils  for  studies  relating 
to  thin  film  energetics.  We  are  also  fabricating  an  evaporation  unit  to 
be  incorporated  into  a  standard  evaporation  chamber  to  accommodate 
some  pilot  studies  on  the  vapor  deposition  of  gadolinium,  with  the  view 
towards  studying  dynamic  recrystallization  and  gram  growth  kinetics  in 
the  pure  foil  or  its  oxide  following  the  previous  work  of  Murr  , 

Initial  designs  have  also  been  drawn  for  the  fabrication  of  a 
special  vacuum  unit  having  the  capability  of  aiding  in  the  investigation 
of  the  influence  of  electric  and  magnetic  fields  on  the  nucleation  and 
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growth  of  vapor  deposited  thin  films. 
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2.  APPLIED  ELECTROMAGNETICS  AND  PLASMAS 


2.1  PLASMAS 

2.1.1  Reflection  and  Transmission  of  Waves  from  Magnetized 
Nonuniform  Plasma  Slabs 
GK  1115,  National  Science  Foundation 
AF-AFOSR-496-67,  _■  mt  Services  Electronics  Program 
H.  H.  Kuehl,  3.  B.  O'Brien 

The  study  of  the  plasma  slab  formed  by  using  a  pulsed  high- 

frequency  transmitter  as  the  energy  source  has  been  continued*.  The 

2 

Buchsbauin-Hasegawa  resonances  have  been  observed  in  several  noble 
gas  afterglow  plasmas.  Probe  measurements  are  presently  being  carried 
out  to  determine  the  density  profile  and  temperature  during  the  afterglow. 
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2.1.2  Photolithography  Laboratory  for  the  Fabrication  of 
Microwave  Devices 

AF-AFOSR -496- 67 ,  Joint  Services  Electronics  Program 
W.V.T.  Rusch,  E.  D.  Garlinger 


Complete  facilities  for  the  fabrication  of  microwave  devices 
using  photolithographic  planar  photoresist  techniques  have  been  established 
and  are  currently  in  operation.  These  facilities  include: 


2. 

3. 

4. 

5. 

6. 

7. 

8. 


A  complete  furnace  system  for  the  deposition  of  SiO^ 
on  semiconductor  substrates. 

A  spinner  for  the  application  of  photoresist  to  substrates. 
Masks,  mask  alignment  and  exposure  system. 

Al!  requisite  chemical  facilities  for  the  photoresist 
processes  and  other  subsequent  procedures. 

Electrolytic  baths  and  a  vacuum  strip-heater  for  the 
preparation  of  back  contacts. 

Electrolytic  baths  and  facility  to  fabricate  Schottky- 
barrier  interfaces  on  planar  structures. 

Electronics  to  test  the  DC  characteristics  of  the  devices. 
Mechanical  equipment  to  mount  the  finished  devices  in 
appropriate  waveguide  configurations. 


Activities  to  the  present  time  have  been  primarily  concerned  with  the 
fabrication  of  one  to  two  micron  diameter  Schottky-barrier  diodes  for  90- 
GIIz  mixer  applications. 
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Additional  project*  have  involved  electrolytic  procedures  for 
the  plating  of  gold  and  tin. 
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2.2  MILLIMETER  WAVE  RADIOMETRY 


2.2.1  Millimeter-Wave  Radiometry  for  Radio  Astronomy 

NASA7-I00,  National  Aeronautics  and  Space  Administration 
952  210,  Jet  Propulsion  Laboratory 
W.V.T.  Rusch,  A.  Cooper,  S.  Slobin 

The  total  integrated  opacity  of  the  atmosphere  can  be  deter¬ 
mined  by  measuring  the  radiation  from  an  extra-terrestrial  microwave 
source  over  a  range  of  zenith  angles.  The  sun  is  a  strong  source  of 
rnillimeter-wave  radiation  that  can  be  used  for  this  purpose.  In  the 
interest  of  acquiring  large  amounts  of  atmospheric  data,  a  semi-automated 
sun  tracking  system  (STS)  is  being  constructed.  The  data  obtained  in  this 
manner  can  be  correlated  with  ground-level  atmospheric  parameter  data. 
Correlations  of  these  two  sets  of  data  should  provide  information  useful 

in  many  types  of  communication  systems  such  as,  for  example,  satellite 

1 

repeater  stations.  Furthermore,  the  present  nodding  subdish  system  , 
because  of  its  beam  switching  capabilities,  should  be  particularly  well 
suited  for  such  a  sun-tracking  system. 

The  STS  will  consist  basically  of  the  existing  90-GHz  radio 
telescope  with  appropriate  modifications.  The  antenna  will  be  tilted  into 
its  equatorial  position.  Accurate  automatic  tracking  of  the  sun  will  be 
provided  by  a  closed-loop  control  ystem  such  as  the  one  proposed  below. 

The  error  sensing  element  will  be  a  refracting  telescope 
consisting  of  a  two-inch  achromatic  objective  having  a  24-inch  focal  length 
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and  a  Barlow  lens  having  a  magnification  of  4X.  The  effective  focal 
length  of  this  optical  system  will  then  be  120  inches.  The  diameter  of 
the  sun's  image  produced  by  this  telescope  will  be  one  inch;  movement 
of  the  sun  through  one  minute  of  arc  will  move  the  image  by  0.3  5  inches. 
Four  photoconductive  cells  will  be  placed  in  the  focal  pl?ne  of  fhis  optical 
system,  as  indicated  in  Figure  1.  When  the  sun  is  centered  evenly  on  the 
four  cells,  the  resistances  of  all  cells  will  be  equal,  and  no  error  signal 
is  produced.  If  the  sur's  image  should  shift,  the  cells  would  be  illumi¬ 
nated  unequally,  thereby  producing  an  error  signal.  Suitable  optical 
filters  will  be  provided  to  obtain  optimum  response  from  the  photocon¬ 
ductive  cells. 

Figure  2  shows  the  electrical  equipment  necessary  to  convert 
error  signals  from  the  photoconductive  cells  into  control  pulses.  There 
control  pulses  will  then  be  fed  into  SLO-SYN  stepping  motors  located  on 
the  right  ascension  and  declination  areas  of  the  telescope.  These  motors 
will  redirect  the  telescope  toward  the  sun.  The  photoconductive  cells 
controlling  the  right  ascension  axis  will  also  control  the  speed  of  a 
sidereal  tracking  motor.  This  system  can  be  converted  into  an  open-loop 
tracking  system  with  appropriate  modifications. 

The  electrical  equipment  shown  in  Figure  2  has  been  con¬ 
structed  and  installed  in  the  telescope  drive  system.  The  teslescope 
and  its  associated  equipment  for  the  photocells  has  been  constructed. 

This  telescope  is  currently  being  installed  on  the  90-GHz  radiometer. 
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Sciences  Laboratory,  School  of  Engineering,  University  of  Southern 
Californir . 


2.2.2  Phase-Center  Relationships  in  an  Asymmetric 
Cassegrainian  Feed-System 

NASA?- 100,  National  Aeronautics  and  Space  Administration 
952  210,  Jet  Propulsion  Laboratory 
W.V.T.  Rusch,  S.  Slobin 


Introduction 

In  an  attempt  to  develop  a  beam- switching,  synchronous 
detection  radio  telescope  for  use  in  radio  astronomy,  the  problem  of 
electromagnetic  wave  scattering  from  a  tilted  hyperboloidal  subreflector 
in  a  Cassegrainian  antenna  system  was  encountered.  Beam  switching  is 
accomplished  by  periodically  nodding  the  subreflector  from  one  asymme- 
tri<-  position  to  another  at  the  synchronous  detection  ^ate.  It  is  necessary 
to  know  the  position  of  the  phase  center  of  the  scattered  field  in  order  to 
accurately  predict  the  displacement  and  position  of  the  resultant  primary 
antenna  beam. 

The  methods  of  phase  center  determination  developed  in  this 
report  utilize  techniques  of  geometrical  optics  and  diffraction  theory,  in 
several  different  forms.  Which  method  is  used  depends  to  a  large  extent 
on  the  accuracy  of  the  result  desired  and  the  availability  of  high-speed 
computing  facilities.  However,  each  method  presented  here  gives  a 
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determination  of  phase  center  position  with  considerable  accuracy.  *  j 

All  phase  center  determinations  are  for  the  specific  geometry 
of  the  nodding  subdish  system.  Various  hyperboloid  D/X  ratios  are 
considered  also. 

Geometrical  Concepts 

Figure  1  indicates  the  geomet.  /  of  the  tilted  hyperboloid 
system.  In  contrast  to  the  symmetn  '  g<  ornetry,  where  the  axes  of  both 
the  hyperboloid  (F'F")  and  the  paraboloid  (  IF)  are  colinear,  the  source 
of  spherical  electromagnetic  waves  (Q)  no  longer  lies  on  the  hyperboloid 
axis  of  symmetry  (F'F").  Three  reflected  rays  R^,  R^,  and  R^  are 
shown.  With  these  reflected  rays  are  associated  virtual  rays  V^,  V^> 
and  V.j.  In  the  tilted  configuration  there  is  no  specific  geometrically  _xact 
point  of  intersection  of  all  the  virtual  rays.  There  is,  rather,  an 
envelope  of  points  or  caus tic  formed  by  the  intersections  of  all  neighbor¬ 
ing  virtual  rays.  The  geometrical  foci  of  the  tilted  hyperboloid  are  the 
points  F'  (virtual)  and  F"  (real).  The  prime  focus  of  the  paraboloid  is 
at  point  F.  Neither  F'  nor  F  may  be  regarded  as  the  hase  center  of  the 
scattered  field  in  the  tilted  configuration. 

Methods  of  Fhase  Center  Determination 

A.  Determination  of  the  Caustic  Curve.  It  can  be  shown 
by  a  simple  geometrical  analysis  that  the  caustic  curve  represents  the 
locus  of  localized  phase  centers  of  sets  of  adjacent  rays  scattered  from 
the  hyperboloid.  Consequently,  it  is  informative  to  investigate  the 
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position  of  the  caustic.  The  equations  of  the  caustic  curve  are  obtained 
from  the  principles  of  differential  geometry.  It  is  found  that  the  caustic 
curve  has  two  branches,  each  of  which  corresponds  to  the  loci  of  the 
two  principal  centers  of  curvature  for  each  point  of  the  reflected  wave- 
front.  Of  special  interest  is  the  point  of  intersection  of  the  two  branches 
of  the  caustic  curve,  because  of  its  nearness  to  the  subsequently  deter¬ 
mined  phase  centers  of  the  total  scattered  field.  The  position  of  the 
caustic  intersection  point  is  listed  in  Table  I  and  shown  in  Figure  2. 

B.  Phase  Center  Determination  by  Geometrical  Optics.  By 
carrying  out  a  ray  tracing  procedure  .a  the  case  of  the  asymmetric  geome¬ 
try,  a  path  length  "phase"  of  the  reflected  field  may  be  calculated  on  a 
ciicle  centered  at  the  origin  of  the  til tt d  (X ' -  Z ' )  coordinate  system.  For 
an  asymmetric  geometry,  the  phase  on  this  circle  will  not  be  constant. 

A  best-fit  phase  center  point  may  be  postulated  to  exist;  and  this  point  is 
found  by  a  least-squares  curve  fitting  technique  in  which  the  rms  differ¬ 
ence  between  the  phase  on  circles  cen  ered  at  point  F’  and  the  phase  on 
circles  centered  at  the  postulated  best-fit  phase  center  is  minimized.  A 
weighting  technique  is  used  to  account  for  the  non-uniform  illumination 
from  the  source  Q.  The  position  of  this  phase  center  point  and  all  other 
phase  center  points  is  shown  in  Figure  2. 

C.  Phase  Center  Determination  by  Fhysical  Optics.  A  more 
accurate  representation  of  the  scattered  field  for  moderately  sized  hyper¬ 
boloids  can  be  obtained  using  the  principles  of  physical  optics.  Using 
results  of  a  diffraction  theory  analysis  of  scattering  from  a  tilted 
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No. 

Analys  is 

Radius  (a), 

wavelengths 

. .  —  .  — .  ,  ■  ■ 

Radius  (a), 
millimete  rs 

I 

Angle  (*) 

!■  ■  -r  ■  a 

- -i 

l 

Geometrical  Optics  Best-Fit 
Weighted,  1 17  pts. 

- 

5.  7  1 4  mm 

87.  15° 

2 

Diffraction  Theory  Best-Fit 
Weighted,  121  pts.,  D/X  = 

59.  54 ,  X  r  3 .  33  mm 

_ 

1. 707  X 

5.  684  mm 

88.  11° 

!  3 

| 

| 

- ? - - - 1 - - - 

|  } 

Diffraction  Theory  Best-Fit 

W e igh ted ,  121  pts.,  D/X-  0.689X  i  5.691  min 

24 ,  X  =  8.  26  mm  1  j 

!  i 

89.  37° 

1  4 

i 

Diffraction  Theory  Best-Fit  i 

Weighted ,  121  pts.,  D/X  -  10  ;  0.289X 

1  -19.83  mm  \ 

j 

L_  _ 

5.  73  1  mm 

94. 04° 

5  Intersection  of  Caustic  Curve  -  5.  32  1  mm 

88. 76° 

■ - - 

6 

Apolication  of  Demagnification' 

Formula  j  1.595X  5.318  mm 

D  /  X  :  59.  54 ,  X  =  3,33  mm 

1  i 

_ i _  _ 

87.  82° 

Table  1 


hyperboloid,  the  phase  of  the  scattered  field  may  bo  specified  on  circles 
centered  at  point  F,  the  origin  of  the  untilted  (X-Z)  coordinate  system. 

A  lea  S  t  -  squa  r  es  curve  fitting  technique,  identical  to  that  described  in 
Section  II,  is  utilized  to  determine  a  best-fit  phase  center  for  the  physmal 
optics  vase.  Hyperboloid  D;  1  ratios  of  M  .  24,  and  iO  are  i  onsidered. 
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unfitted  coordinate 
system 


and  phase  n  rti-r  positions  for  the  thr«.e  cases  are  tound  to  he  within 
approximately  0.7  mm  of  each  other.  A  t  if  teen- fold  deeriaae  in  rir  s 
phase  variation  around  the  best-fit  phase  centers  has  been  determined, 
compared  to  the  phase  variati  is  around  the  paraboloid  prime  focus  (F). 

D  Phase  Center  Determination  from  Demagnification 
Formulas.  The  hype rboloidal  subrefiector  in  a  Ca s s eg ra  1  nia n  antenna 
system  acts  tike  a  demagnifying  optical  system.  Following  optical  rules 
regarding  image  movement  as  a  function  of  source  movement,  an  optically 
equivalent  system  may  be  set  up  for  the  tilted  hyperboloid  geometry.  The 
simple  rules  yield  a  remarkably  accurate  determination  of  phase  center 
position.  Tt  is  to  be  expected,  however,  that  for  more  extreme  tilt  angles 
than  that  in  the  nodding  subdish  system  (2.0b  deg)  the  demagnification 
formulas  will  produce  less  accurate  results  than  those  obtained  by 
previously  discussed  methods. 

Summary  of  Results  and  Conclusions 

The  coordinates  of  all  phase  centers  determined  by  the 
methods  of  the  previous  section  are  shown  in  Figure  2  and  tabulated  in 
Table  1.  Not  shown  in  the  figure  is  the  actual  physical  position  of  the 
hyperboloid  geometric  focus  (F1)  in  its  tilted  position.  The  position  of 
F1  is  almost  3mm  below  the  points  shown  in  Figure  2 

In  general,  the  phase  center  determination  by  the  methods  of 
physical  optics  yields  the  most  accurate  results;  however,  for  D/X  ratios 
of  10  or  larger  and  moderate  tilt  angles,  the  results  of  all  techniques 
agree  tc  within  0.1  to  0.2  wavelengths. 
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The  geometrical  focus  F'  of  the  tilted  hyperboloid  may  not 


be  regarded  as  the  phase  center  of  the  scattered  field. 

The  computed  phase  center  displacement  may  be  used  in 
conjunction  with  standard  beam-deviation-facto»*  formulas  to  estimate 
the  displacement  of  the  primary  antenna  beani  with  considerable  accuracy. 

2.2.3  Pioneer  VI  Solar  Faraday  Rotation  Experiment 
J.  E.  Ohlson,  W.  V.  T.  Rusch,  C.  T.  Stelzried 


Introduction 

In  T  ivember  1966,  the  Pioneer  VI  Spacecraft  will  undergo 
solar  occultation  as  it  passes  approximately  one  earth  orbital  radius 
behind  the  sun.  During  the  immersion  and  emersion  phases  of  the 
occultation  the  spacecraft  radio  signals  to  earth  will  pass  through  the 
solar  corona.  This  offers  a  unique  opportunity  to  measure  the  Faraday 
Rotation  undergone  by  a  linearly  polarized  electromagnetic  wave  upon 
passing  through  the  solar  corona.  The  data  ootamed  will  give  imormation 
on  the  free  electron  density  and  magnetic  field  in  the  corona.  This 
experiment  is  being  performed  jointly  with  personnel  of  the  Communica¬ 
tions  Elements  Research  Section  of  the  Jet  Propulsion  Laboratory,  G.  S, 
Levy,  B.  Seidel  and  T.  Sato,  in  support  of  a  NASA  contract  150-22-13- 
05033. 

Theory 

The  rotation  of  the  plane  of  polarization  of  a  linearly 
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polarized  elec tr omagnetic  wave  upon  pa  ->mg  along  a  path  in  a  magneto- 
ionic  medium  is 

Q  =  — C^S—  j  N(s)  B(s)  •  ds 
f  path 

where  f  is  cyclic  frequency,  N(s)  i&  li<_e  electron  density  and  I  ( s )  is  the 
vector  magnetic  field.  By  obtaining  continuous  data  on  Cl  for  several 
hours  each  day  during  the  occulta  Hon  phases,  a  large  body  of  data  will 
be  available  for  fitting  to  various  models  of  the  corona. 

Measurement  Techniques 

Several  physical  constraints  make  the  measurement  of  Q  ver 

difficult: 

(1)  The  270,000,000  kilometer  spacecraft-earth  distance 
produces  an  extremely  weak  signal  at  earth  (-161  dBm). 

(2)  As  the  spacecraft  nears  the  sun  (in  angle),  intense  solar 
radio  noise  enters  the  antenna  pattern  sidelobes  and  it  appears  that 
tracking  the  spacecraft  closer  than  about  1°  to  the  sun  will  not  be 
possible. 

(3)  Coronal  plasma  turbulence  produces  spectral  spreading 
of  the  signal,  which  makes  signal  tracking  difficult.  This  effect  also 
worsens  as  the  sun  is  reared. 

(4)  The  earth's  ionosphere  contributes  a  varying  rotation  of 
1-5°  to  the  signal.  This  effect  must  be  determined  to  allow  corrections 
of  the  data. 
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(5)  An  apparent  rotation  of  the  spacecraft  is  caused  by  the 
earth's  rotation  and  also  by  the  motion  of  th^  tracking  antenna,  since  it 
has  an  Az-El  instead  of  polar  mounting. 

The  210  foot  diameter  NASA/JPL  antenna  at  Goldstone  with 
maser  receivers  is  the  most  sensitive  facility  in  the  world  for  S-band 
spacecraft  signals.  Extensive  modification  of  the  r.  f.  hardware  and 
receiver  configuration  has  produced  an  antenna  with  a  mechanically 
driven  cross-polarized  feed  system.  Using  monopuls e- rada r  circuitry, 
one  feed  produces  an  error  signal  which  serves  the  primary  feed  to  the 
received  polarization.  The  primary  feed  thus  receives  maximum  signal 
for  a  phase-loc ked-loop  for  frequency  tracking  and  determination  of  a 
detection  reference  phase  for  the  error  channel.  The  mechanical  position 
of  the  feeds  is  the  estimate  of  received  signal  polarization.  This  is 
digitally  recorded  and  later  transformed  {taking  into  account  the  effects 
of  (5)  above)  by  a  computer  to  give  polarization  rotation.  Measurement 
accuracy  is  not  limited  by  mechanical  errors,  but  by  signal  to  noise  ratio, 

A  prediction  of  0  indicates  that  it  will  reach  a  maximum  of 
o 

about  10  before  real-time  tracking  of  the  signal,  as  described  above, will 
not  be  feasible  due  to  sm«iu  SNR.  At  that  point,  open-loop  recording  at 
i.f.  will  be  msde  with  the  polarizer  feeds  in  a  fixed  position  ?n  i  the 
receiver  using  a  fixed  local  oscillator.  Computer  spectral  analysis 
processing  will  allow  following  the  signal  uown  to  a  smaller  SNR  than 
the  real-time  system  will  permit.  The  reason  for  not  using  open-loop 
recording  and  processing  for  all  data  is  that  phase  and  amplitude 
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instabilities  of  the  masers  and  receiver,  while  negligible  tor  most 


applications,  produce  first-order  errors  in  angle  measurements  in  post¬ 
experiment  processing.  These  effects  are  negligible  in  the  real-time 
s  ys  tern . 


Specific  Accomplishments 

(1)  Design  of  antenna  feeds  and  receiver  configuration. 

(2)  Nearly  completed  implementation  of  (1). 

(3)  Analysis  of  system  measurement  capability.  A  good 
approximation  to  the  error  in  measurement  of  polarization  angles  has 
been  shown  to  be 


a.  - 

H 


C  ^  1 


where  is  r.m.s.  error  in  degrees,  T  is  system  operating  tempera¬ 
ture  in  °K,  and  t  is  the  system  response  time  in  seconds. 

(4)  A  computer  program  has  been  implemented  to  determine 
refractive  effects  of  the  corona. 

(5)  A  set  of  computer  programs  to  allow  transformation  of 
pola rization  data  between  terrestrial  and  celestial  coordinates  has  been 
written. 

(6)  A  receiver  with  a  separate  antenna  has  been  installed  to 
measure  the  rotation  produced  by  the  earth's  ionosphere  in  order  to 
remove  this  effect  from  the  data.  This  is  being  done  by  tracking  a 
stationary  s  tellite. 
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(7)  Preliminary  data  indicate  a  system  accuracy  in  agree¬ 
ment  with  (3). 

Future  Plans 

(1)  The  final  experiment  configuration  will  shortly  be 

completed. 

(2)  Data  will  be  taken  through  the  occultation  period. 

(3)  The  second  phase  of  the  experiment  will  begin  as  the  data 
are  evaluated,  coronal  models  are  determined  and  solutions  of  coronal 
parameters  are  found.  This  will  involve  extensive  use  of  data  from 
other  experimenters  who  have  measured  coronal  parameters  in  other 
ways. 


1/ 


/ 
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3.  INFORMATION  SCIENCES 


3.1  CONTROL  SYSTEMS 

3.1.1  Optimal  Control  Problems 

AF-AFOSR-  1029-67A  ,  Air  Force  Office  of  Scientific 

Research 

L.  W.  Neustadt 

Necessary  conditions  for  optimal  control  problems  with 
restricted  phase  coordinates  and  delay-differential  equation  constraints 
were  obtained  and  were  published  in  the  Fh.  D,  dissertation  of  Dr.  S.  C. 
Huang.  These  conditions  extend  the  well-known  Pontryagiu  maximum 
principle.  This  work  was  extended  by  Dr.  Huang  to  obtain  necessary 
conditions  for  optimal  control  problems  where  there  are  constraints  in 
the  form  of  integral  equations.  Investigations  have  begun  on  optimal 
control  problems  with  functional  differential  eq  lation  constraints.  The 
results  that  will  be  obtained  will  hopefully  include  those  of  Huang  as 
spec  ial  cases . 

Three  articles,  describing  the  results  of  previous  research 
which  consider  conventional  optimal  control  problems  from  a  imw  view¬ 
point  and  which  also  give  a  general  abstract  theory  ol  extremality,  have 
been  completed  by  L.  W.  Neustadt. 
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3.1.2  Optimum  Inputs  for  Parameter  Estimation 

NGR-05-018-044,  National  Aeronautics  and  Space 

Administration 

N.  E.  Nahi,  D.  E.  Wallis.  Jr. 

The  quality  of  estimates  of  parameters  of  nonlinear  dynamic 
systems  is  effected  by  the  inputs  to  the  system.  It  is  then  desirable  to 
choose  the  input,  from  an  admissible  class,  in  order  to  enhance  the 
quality  of  the  optimum  estimator. 

A  method  for  finding  optimal  deterministic  inputs  for  estima¬ 
tion  of  system  para  'eters  from  noisy  output  data  is  given.  Analytical 
difficulties  connected  with  the  Bayesian  approach  to  this  joint  optimal 
control  and  optimal  estimation  problem  are  overcome  by  postulating  that 
an  "efficient"  estimator  of  the  parameter  exists.  For  such  estimators, 
the  Cramer-Rao  lower  bound  (CRLB)  is  equal  to  the  optimal  (minimum) 
estimation  variance,  and  its  formulation  is  shown  to  be  a  functional  only 
of  the  system  input.  'Itw-refore,  if  an  input  minimizes  the  CRLB,  it  is 
(together  with  the  postulated  efficient  estimator)  the  optimal  input  for 
joint  minimum- variance  system  parameter  estimation.  Minimization  of 
the  CRLB  is  made  possible  by  a  sensitivity-theoretic  extens.on  of  a 
result  of  Slepian  to  the  case  where  a  signal  is  generated  as  the  output  of 
a  differential  equation  containing  the  parameter.  This  extension  leads 
to  a  deterministic  integral- quad  ratic -c  os  t  optimal  control  problem  which 
may  be  solved  by  establ.shed  techniques.  In  a  certain  sense,  the  method 
presented  enables  determination  ol  an  input  which  optimizes  the  potential 
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for  subsequent  niinimuni-variance  parameter  estmiat'.or., 


Thv  required 


efficient  estimator  is  not  derived. 

3.1.3  Wide  Angle  Optical  Tracking 

NCR  -  05-018-044  ,  National  Aeronautics  and  Space 
A  dm  mist  ra  tion 
N .  E .  N  a  h  1 

The  application  of  linear  (Kalman  -  Rucy)  filtcrir.i;  to  estima¬ 
tion  of  target  trajectory  is  well  known.  Typically  the  problem  is  formu¬ 
lated  as  follows,  "’lu*  unknown  trajectory  is  presented  as  the  output  of  a 
dynamical  system.  The  possible  variations  of  this  trajectory  due  to 
m.predic  ta  ble  variations  of  acceleration  and  other  effects  are  character¬ 
ized  by  random  inputs  to  this  dynamical  system.  Some  function  of  the 
trajectory  is  observed  by  means  of  an  instrument,  such  as  Radar.  This 
observation  is  always  corrupted  by  noise.  An  estimator  (filter)  is  then 
designed  to  operate  on  the  noisv  observation  and  yield  a  best  estimate  of 
the  target  trajectory  (a  present  or  future  value),  recursively. 

What  is  assumed  in  such  an  estimator  as  described  above  is 
that  any  observation  belongs  to  the  trajectory  (which  we  are  estimating) 
with  certainty.  In  other  words,  any  time  an  observation  is  received,  it 
we  exclude  the  contribution  oi  the  observation  noise,  we  assume  that  we 
have  received  a  known  function  of  the  trajectory. 

Let  us  consider  the  problem  of  tracking  a  target  ;n  a  wide- 
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angle  field  of  view.  Now  any  observation  is  either  a  function  of  the 
trajectory  (of  course  corrupted  by  observation  noise)  or  it  may  not 
belong  to  that  trajectory  at  all.  Examples  of  these  appear  in  tracking 
a  target  existing  on  an  optical  or  IR  picture.  The  observation  may  then 
consist  of  detecting,  as  time  proceeds,  bright  points  on  the  picture. 
However,  at  any  time,  false  bright  points  (false  alarms)  may  appear 
which  are  not  connected  with  the  trajectory  being  tracked.  Consequently, 
a  probability  is  associated  with  each  observation  indicating  the  probability 
that  the  obs ervation  belongs  to  the  trajectory.  A  recursive  filter  is  then 
derived  on  the  basis  of  utilizing  this  information. 

Reference 

1.  N.  E.  Nahi,  Estimation  Theory  and  Applications,  John  Wiley  &  Sons, 
1968. 

3.1.4  Optimum  Design  of  Vehicles  with  Constrained  Observation 
NGR  05-018-044,  National  Aeronautics  and  Space 
Administration 

F04(695)-67-C- 109,  Space  Systems  Division,  Air  Force 
Systems  Command 
N.  E.  Nahi,  T.  Cooper 

The  feedback  operation  of  many  systems  is  based  on  the 
information  obtained  by  various  on-board  sensors.  Often  these  sensors 
have  discrete  outputs,  such  as  the  output  of  a  track-scan  rad.*r  system. 
Due  to  various  effects  of  random  noise  there  is  a  non-zero  probability 
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that  a  specific  observation  in  the  lulure  may  not  be  received.  Clearly, 
the  optimum  design  of  the  vehicle  should  utilize  this  critical  information. 
In  other  cases,  the  designer  may  have  influence  over  the  observation 
scheme  used.  For  example,  he  may  increase  the  observation  sampling 
rate.  However,  this,  in  general,  results  in  more  complex  devices 
(such  as  those  yielding  higher  scan  rate)  and  higher  data  processing 
capability  (due  to  the  larger  volume  of  data  to  be  processed).  It  is  then 
logical  to  associate  a  "cost"  to  each  observation  and  to  incorporate  this 
cost  into  the  criterion  cf  system  design. 

The  problem  of  design  of  linear  feedback  systems  with 
quadratic  cost  function, and  when  both  the  system  and  the  observations 
include  additive  random  noise, has  been  satisfactorily  solved  over  the 
last  decade.  It  has  been  shown  that  the  optimum  design  involves  a 
cascade  of  a  Kalman  filter  and  an  optimum  gain  matrix.  Procedures  for 
derivation  of  the  gain  matrix  have  been  derived  leading  to  feasible  com¬ 
puter  algorithms.  The  problems  posed  in  introduction  can  be  considered 
as  extension  and  modification  of  these  results. 

Various  approaches  to  solution  of  the  mentioned  problem  are 
considered.  Basically  the  pro'  fure  of  dynamic  programming  is  utilized 
At  present,  this  attempt  has  yielded  partial  success  toward  solving  the 
first  problem,  namely,  design  of  optimum  linear  discrete  systems  with 
given  probabilities  for  making  future  observations.  The  second  problem 
is  transformed  into  a  linear  system  with  a  new  cost,  by  means  of  a  ?socia 
ting  a  cost  to  making  observations.  It  is  expected  that,  in  conjunction 
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with  the  derivation  of  the  optimum  feedback  control  system,  we  arrive 
at  optimum  observation  schemes,  such  as  how  often  it  is  necessary  to 
make  an  observation, 

3.1,5  Design  of  Random  Circuits 

AF-AFOSR- 1029-67,  Air  Force  Office  of  Scientific 

Research 

D.  F.  Elliott  and  D.  D.  Sworder 

Where  applicable, s tcchas tic  approximation  provides  a  tool 
which  permitB  the  engineer  to  include  component  uncertainties  in  his 
design  work.  The  question  of  convergence  of  this  algorithm  in  circuit 
applications  was  explored  in  [  1 1  using  Dvoretzky's  fundamental  theorem. 
Recently  the  problem  posed  in  [  i  ]  has  been  investigated  and  conditions 
sufficient  to  guarantee  convergence  have  been  given  a  sensitivity  inter¬ 
pretation. 

Suppose  an  engineer  wishes  to  synthesize  a  circuit  to  perform 
a  specified  task  and  suppose  that  the  performance  measure  is  denoted  by 
h.  The  performance  of  the  circuit  depends  upon  its  component  elements, 
some  of  which  are  imperfectly  known.  Let  the  vector  of  system  param¬ 
eters  to  be  selected  by  the  circuit  designer  be  denoted  by  b.  For  example, 
if  the  behavior  of  a  circuit  is  very  sensitive  to  a  pair  of  "10-percent" 
resistors,  b  might  oc  a  two  component  vector  with  the  first  component  equal 
to  the  nominal  value  of  the  first  resistor.  The  performance  measure  can 
be  written  as  h(k>). 


-128- 


The  classical  methods  of  circuit  design  simply  neglect 
component  tolerances  and  the  designer  selects  as  best  that  parameter 
vector  svhich  minimizes  h(b).  In  truth,  however,  choice  of  b  only  fixes 
the  probability  distribution  function  of  the  random  variable  h(b)  since 
the  parameters  actually  realized  in  the  circuit  will  deviate  from  their 
nominal  values.  A  more  appropriate  performance  measure  might  be 
the  expected  value  of  h{b); 

H(b)  =  E  {  h(  b) }  . 

To  be  more  precise,  suppose  that  there  are  q  system  design 
parameters  (resistor  values,  capacitor  values,  amplifier  gains,  etc.). 

Any  specific  circuit  configuration  is  described  by  a  q-dimensional  vector 
b  with  components  b.;  i-i  ,  .  .  .  ,q.  Suppose  that  there  are  constraints  on 

the  allowable  parameter  values: 

b.  e  [  a. ,  g.  ]  ;  i-1 . q  , 

l  1  1 

where  Ij  (3  -  a  jj  <  «,  Denote  the  q-fold  Cartesian  product  of  intervals 
[a  ,  S,.i  by  C(a,6).  The  engineer  then  seeks  that  _G  e  C(Ct  ,  B)  such  that 

H(  0  )  =  inf 

be  C(a  ,  8) 

In  I  Z  ]  it  is  shown  that  under  suitable  restric  tions,  a  modified 
Kiefer-Wolfowitz  procedure  will  generate  a  sequence  of  parameter 
settings  which  converge  to  >n  \sith  probability  one.  The  most  interesting 
restriction  is  one  which  constrains  the  relative  magnitude  ot  the  eigenvalues 
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of  the  Hessian  matrix  of  the  regression  function  at  8.  The  implication 
of  this  result  is  that  convergence  may  not  occur  if  the  system  performance 
is  overly  sensitive  to  a  circuit  parameter. 


References 


1.  K.  B.  Gray,  "The  Application  of  Stochastic  Approximation  to  the 
Optimization  of  Random  Circuits",  Proceedings  of  Symposia  in 
Applied  Mathematics,  Vol.  26,  pp.  178-192,  1964. 
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3.1.6  Control  of  Stochastic  Systems 

AF-AFOSR  - 1029-67,  Air  Force  Office  of  Scientific 

Research 

D.  D.  Sworder 


Recently  it  has  been  shown  that  the  maximum  principle  may 
be  employed  to  design  controllers  for  stochastic  systems".  Indeed,  it 
was  demonstrated  that  the  differential  equations  for  the  adjoint  variables 
are  identical  to  those  of  Pontryagir*  when  the  random  disturbance  is  addi¬ 
tive  white  noise.  This  success  with  the  problem  posed  in  the  indicated 
reference  might  lead  one  to  hope  that  perhaps  the  maximum  principle 
formalism  can  be  extended  to  a  wide  class  of  stochastic  systems. 

Actually,  this  extension  may  nu!  b*>  made  "  ithout  introduc  ing 
additional  complexity.  While  a  technique  similar  to  the  maximum  princi¬ 
ple  may  be  employed  to  find  the  best  control  policy  for  an  important  class? 
of  systems  with  multiplicative  disturbances,  fundamental  changes  must  be 
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mad?  in  the  equations  for  the  adjoint  variable  .  Instead  of  leading  to  a 
set  ox  stochastic  ordinary  differential  equation:-  as  in  Kushner:s  work,  the 
formalism  leads  to  a  set  of  stochastic  partial  differential  equations. 

The  motivation  for  this  work  stems  from  a  need  to  design 
controllers  for  systems  whose  components  are  subject  to  failure.  The 
natural  model  for  plants  of  this  type  n  st  include  multiplicative  distur¬ 
bances  of  the  jump  type.  If  random  variations  occur  in  the  system,  a 
feedback  controller  is  necessary  to  insure  adequate  performance  under 
ali  operating  conditions.  The  question  then  arises  as  to  how  best  to 
process  the  information  available  to  the  controller. 

In  [  3  ]  and  [4]  several  interesting  cases  in  which  the  design 
equations  could  be  solved  explicitly  were  studied.  For  analytical  conven¬ 
ience, attention  was  restricted  to  linear  plants  although  the  contr  Her  and 
hence  the  total  system  was  nonlinear.  It  should  be  explicitly  noted  that 
the  class  of  problems  was  not  restricted  to  Markov  decision  problems. 
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5.1.7  Fuel  Optimum  Spacecraft  Attitude  Control 


NGR  05-018-044,  National  A.eronautics  and  Sj.-.^ce 
Administration 

H.  J.  Payne,  R,  3.  McGhee,  W.  H.  Spuck 

Since  the  mission  life  of  long  duration  space  flights  ts  often 
limited  by  fuel  supply,  it  becomes  important  to  utilize  fuel,  a  'or 
attitude  control,  in  optimal  or  near  optimal  fashion.  It  has  been  the 
objective  of  this  program  to  determine  for  simple  models  of  attitude 
control  systems  (1)  optimal  control  laws  and  he  related  fuel  expenditure 
rate,  and  (2)  the  relative  efficiency  of  practi  al  controllers,  i.  e.  ,  con¬ 
trollers  which  are  specified  in  such  a  way  that  implementation  is  relatively 
easy. 

As  previously  reported,  the  first  portion  of  the  effort  was 
completed  for  a  system  described  by 

x  =  u  +  n(t) 

in  which  x  represents  angular  position,  n  the  control  force  and  n h'  the 

perturbing  forces.  These  results  have  been  reported  internally*  and 

2 

submitted  for  publication  . 

Work  has  been  initiated  which  is  intended  to  provide  the 
basis  for  the  study  of  practical  controllers.  This  effort  includes  the 
refinement  of  numerical  methods  so  that  results  may  be  obtained  more 
efficiently.  It  was  sufficient  in  the  study  of  optimal  control  laws  to  use 
simpler  number  numerical  methods,  i.  e.  ,  explicit  schemes,  since  a 
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limited  number  of  cases  had  to  be  computed.  But  if  several  practical 
controllers  are  to  be  examined,  more  efficient  numerical  schemes  become 
necessary.  Among  the  schemes  being  ''onsidcred  is  the  implicit  Crank- 
Nicoison  procedure. 

R  pfprpnrps 

1.  R.  B.  McGhee,  H.  J.  Payne  and  W.  H.  Spuck,  "A  Numerical  Method 
in  Fuel  Optimum  Stochastic  Attitude  Control",  Department  of 
Electrical  Engineering  Report  No.  USCEE  284,  University  of 
Southern  California. 

2.  R.  B.  McGhee,  H.  J.  Payne  and  W.  H.  Spuck,  "Fuel  Optimum 
Stochastic  Attitude  Control",  submitted  for  publication  to  PG-AC. 

3.1.8  High  Density  Freeway  Traffic  Control 

AF-AFOSR-4%~67,  Joint  Services  Electronics  Program 
H,  J.  Payne 

Bottlenecks  on  the  freeway  of  recurrent  (arising  because  of 
geometry,  e.  g.  ,  points  of  on-ramp  merging  and  freeway  intersections) 
and  nonrecurrent  (e.  g.  ,  due  to  accidents)  nature  add  considerably  to  the 
transit  times  for  freeway  users.  It  has  been  generally  agreed  that  the 
most  efficient  means  of  influencing  traffu  flow  is  to  control  on-ramp 
flow.  The  intelligent  design  of  on-ramp  control  schemes  will  require 
an  understanding  of  high  density  freeway  traffic  in  the  form  of  reasonably 
accurate  mathematical  models. 

Due  to  the  large  number  of  vehicles  involved,  it  is  appropriate 
to  introduce  and  develop  equations  involving  variables  which  represent 


gross  properties  of  the  traffic  flow.  Appropriate  variables  are  p,  density 
(number  of  cars  per  unit  length  of  road);  q„  flow  rate  (number  of  cars 
passing  a  fixed  observer  per  unit  time);  and  u,  mean  speed  (average 
speed  of  vehicles  in  a  fixed  length  of  road).  The  conservation  of  cars 
equation 


and  the  relation 


p  +  q  =  0 

t  x 

q  -  pu 


must  be  supplemented  by  an  equation  describing  the  dynamic  behavior 
of  traffic. 

The  model 

q  '  Q(P)  -  v  P 

x 

has  been  found  to  be  inadequate  in  the  sense  that  perturbations  in  traffic 
always  dissipate,  precluding,  for  example,  a  description  of  the  develop¬ 
ment  of  "stop  and  go"  traffic. 

The  model  presently  being  studied  is  described  by 


u  +  uu  -  -a(q-Q(p)  +  Suo  )/ „ 

t  X  X 

Stationary  solutions,  in  which  all  variables  are  functions  ot  x-Ut  for 
some  fixed  U  have  been  investigated  and  used  to  describe  some  bottleneck 
situations.  One  interesting  ‘eature  of  these  solutions  is  that  traffic 
breaks  down  to  near  crawl  conditions  if  the  rate  of  traffic  flow  into  a 
bottleneck  exceeds  a  critical  value.  Stability  analysis,  using  a  linearized 
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perturoation  equation  with  uniform  flow  the  standard,  shows  that  uniform 


flow  becomes  unstable  above  a  certain  critical  density. 

Since  the  model  depends  on  parameters  0.  and  8  and  the  Road 
Diagram  function  Q(c)  ,  empirical  data  is  necessary  to  complete  the 
mathematical  description. 


3.  1.9  .*  Control  TH'-r'-y  and  Tv.  o- Feint  Doundaiy- 

Value  Problems 


GP8960,  National  Science  Foundation 
R.  bellman 


Work  has  continued  on  the  development  of  a  stability  theory 
for  two-point  bounda  r  y  -  va  lue  problems  paralleling  the  classical  stability 
theory  ot  Poincare  and  Lyapunov.  The  tollowing  three  papers  are  devoted 
to  this : 

1.  R.  Bellman,  On  Analogues  ot  Poincare-1  yapunov  Theory  lor  Multi- 
point  Bounda  ry- V.-.lue  Problems--!,"  .!.  Math.  Anal.  A  p  p  1  ■  Vol.  14, 

1  'Hiti .  pp.  ^22-  ^d(>. 

2.  R.  Bel lm  in.  "On  Analogues  of  Poincare  -  I  yapunov  Theory  tor  Multi¬ 
point  Bounda  ry- Vaiue  Problems- -Correction,  "  J.  Math,  Anal.  Appl.  . 
to  a ppea r . 

1.  R.  bellman,  A  Note  on  Asymptotic  Control  Theory,"  J.  Math.  Anal. 
Appl .  ,  t  o  a  ppea  i  . 


Secondly,  an  analogue  of  classical  iteration 

theory  ior  two-point  bounda  r  v  -  \  a  lue  problems  was  developed.  See 

R.  Bellman.  "Two-Point  Bounda  ry- Value  Problems  and  Iteration. 
A  equationes  Ma  them.it  n  ac  .  to  appear. 


.'.1.10  Dynamic  Programming  and  Partial  Differential  Equations 
GP-8960,  Nath  nal  Science  Foundation 
AT  (1 1  -  1 )-  1 1 3  ,  Atomic  Energy  Commission 
R.  Bellman,  E.  Angel 


Work  has  continued  in  the  application  of  dynamic  programming 
to  the  analytic  computational  solution  of  nonlinear  partial  differential 
equations  of  elliptic  type.  New  methods,  based  on  dynamic  programming, 
have  been  introduced  which  are  far  simpler  than  any  existing  methods 
and  thus  easier  to  teach  and  carry  out,  and  which  seem  to  be  more  effec¬ 
tive  in  a  number  of  cases.  This  work  is  contained  in  Angel’s  thesis, 
Dynamic  Programming  and  Partial  Differential  Equations,  and  ir.  the 
following  references. 

References 

1.  E.  Angel.  "Disc  rete  Inva r  lar.t  Im heeding  and  Elliptic  Boundary- Value 

Problems  Over  Irregular  Regions",  J.  Math,  Anal.  Appl.  ,  Vol.  Zi, 
1968,  pp.  471-484.  . . ' 

2.  E.  Angel.  "A  Building  Block  Technique  for  Elliptic  Bounda  ry  -  Value 
Problems  Over  Irregular  Regions,  "  J.  Math.  Anal.  Appl.  ,  to  appear. 

3.  E.  Angel ,  "Noniterative  Solutions  of  Nonlinear  Elliptic  Equations,  " 


3.2  COMMUNICATION  AND  RADAR  SYSTEMS 


3.2.1  Image  Processing  and  Coding 

NGR  05-018-044,  National  Aeronautics  and  Space 
Admims  tration 

JPL  9523  12,  Jet  Propulsion  Laboratory 
W.  K.  Pratt  and  H.  C.  Andrews 

The  introduction  of  the  fast  Fourier  transform  algorithm  has 
led  to  the  development  of  the  Fourier  transform  image  coding  technique 
whereby  the  two  dimensional  Fourier  transform  of  an  image  is  transmitted 
over  a  channel  rather  than  the  image  itself.  This  development  has  further 
led  to  a  related  image  coding  technique  in  which  an  image  is  transformed 
by  a  Hadamard  matrix  operator.  The  Hadamard  matrix  is  a  square  array 
of  plus  and  minus  ones  whose  rows  and  columns  are  orthogonal  to  one 
another.  A  high  speed  computational  algorithm,  similar  to  the  fast 
Fourier  transform  algorithm,  which  performs  the  Hadamard  transforma¬ 
tion.  has  been  developed.  Since  only  real  number  additions  and  subtrac¬ 
tions  are  required  with  the  Hadamard  transform,  an  order  of  magnitude 
speed  advantage  is  possible  i  ompared  to  U  •>  complex  Rum  Per  Fourier 
transform.  Transmitting  the  Hadamard  transform  oi  an  image  rather 
than  the  spatial  representation  -it  the  image  provides  a  potential  toleration 
to  i  tunnel  errors  and  the  possibility  ot  redm  oil  bandw.dtb.  t  ra  ns  m :  s  s  .on. 

1  ot  the  arras  Six.  v )  represent  the  intinsili  samples  o!  an 
or  1^  :ia!  imago  over  an  arms  o!  N  points .  Then  tho  tv  o  i.m  ns.an.ii 
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Hadamard  transform,  F(u,v),  of  f(x,y)  is  given  by  the  matiix  product 
[  F (u,  v)]  =  [  H(u,  v)  ]  T  H  (u,  v)  ] 

where  [H(u,v)I  is  a  symmetric  Hadamard  matrix  of  order  N.  Pre-  and 
post-multiplication  of  [F(u,v)J  by  the  Hadamard  matrix  gives 

[  H(u.  v;]  [  F(u,  v)l[  H(u,  v)  j  = 

[  H(u,  v)][  H(u,  v)][  f(x,  y)][  H(u,  v)l  [  H(u,  v)"1 

But,  for  a  symmetric  Hacimar  matrix,  FH  =  NI  where  I  is  the  identity 
matrix.  Hence, 

[  f(x,  y)  ]  =  — — r  [  H(u,  v)1  [F{u,  v)]  [  H(u,  v)] 

N 

2 

and,  aside  from  the  constant  scaling  factor,  N  ,  the  arrays  f(x,y)  and 
F(u,v)  are  two  dimensional  transform  pairs. 

Experiments  have  been  performed  to  determine  the  nature  of 
the  two  dimensional  Hadamard  transform  of  an  image,  and  to  assess  the 
effects  of  a  double  Hadamard  transformation.  Figure  la  contains  a 
cathode  ray  tube  display  'f  the  logarithm  ol  the  magnitude  of  the  ordered 
Hadamard  transform  of  a  Surveyor  spacecraft  scene.  In  this  photograph 
the  origin  (zero  frequency)  appears  in  the  lower  left  corner.  As  the 
frequency  increases  the  magnitude  of  the  samples  tend  to  decrease.  This 
is  an  indication  that  there  are  relatively  few  high  amplitude  brightness 
transitions  between  elements  in  the  original  scenes.  Figu-e  lb  shows 


-138- 


(a)  Logarithm  of  the  Magnitude  of  the 
Hadamard  Transform 


(b)  Double  Hadamard  Transform 


Figure  1.  Hadamard  Transforms  of  an  Image. 
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the  double  Hadamard  transform  cf  the  scene.  There  is  no  noticeable 
image  degradation  between  the  originals  and  the  double  Hadamard  trans¬ 
forms. 

3.2.2  Optical  Coherenc e 

AF-AFOSF  - 188-67,  Air  Force  Office  of  Scientific  Research 
R.  S.  Macmillan 


The  fundamental  objective  of  the  research  is  to  develop  a 
technique  for  maximizing  the  information  transfer  through  an  optical 
system  by  adjusting  the  variable  parameters  of  the  system  optimally. 

In  general,  these  parameters  include  the  amplitude  and  pha-  e  shadings 
of  the  system  apertures.  It  has  bean  shown  that  maximization  of  the 
first-order  and  higher-order  transcoherence  functions  accomplishes  this 
objective. 

The  experimental  investigation  of  both  amplitude  and  intensity 
transcoherence  theory*  is  continuing.  The  experimental  results  indicate 
that  a  U-shaped  amplitude  shading  function,  for  intensity,  is  optimum  for 
maximum  information  transfer.  These  results  are  in  good  agreement 
with  theoretical  results  which  show  that  such  shading  maximizes  the 
transcoherence  function  and  consequently  the  information  transfer. 
Experimental  studies  and  theoretical  analyses  are  continuing  into 
variously  shaped  one-  and  two-dimensional  apertures.  The  experimental 
results  not  only  verify  the  theory,  but  suggest  that  a  generalized  adapta¬ 
tion  of  the  basic  theory  can  be  employed  as  an  optimal  lens  design  technique 
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as  well  as  an  optimization  procedure  for  optical  information-carrying 
systems, 

R eference 

1.  Macmillan,  R.  S.  ,  and  G.  O.  Young,  "Optimization  of  Information- 
Processing  Optical  Systems  by  Transcorrelation-Function  Maximiza¬ 
tion",  the  Journal  of  the  Optical  Society  of  America,  volume  58, 
number  3,  page  346,  March,  1968. 

3.2.3  Synchronization  Coding 

DA-ARO-D  31-124-G929,  U.S.  Army  Research  Office 
R.  A.  Scholtz  and  R.  M.  Storwick 

In  a  coded  communication  system,  synchronization  is  of 
utmost  importance,  and  the  ability  to  attain  and  maintain  sync  is  a 
significant  design  factor.  For  the  discrete  noiseless  channel,  the  usual 
approach  to  this  problem  is  to  encode  the  message  so  that,  whatever  its 
initial  point  of  observation,  the  receiver  will  attain  sync  after  observing, 
at  most.some  bounded  number  of  consecutive  code  symbols. 

Relaxation  of  the  bounded  sync  delay  property  results  in 
codes  known  as  statistically  synchronizable  (ss)  codes.  Under  these 
circumstances,  it  is  possible  to  use  many  more  of  the  available 
sequences  as  codewords,  permitting  the  use  of  shorter  codewords  for 
a  given  dictionary  size.  One  approach  to  ss  code  design  is  to  eliminate 
one  n-tuple  as  a  possible  codeword  and  to  encode  messages  with  the 
remaining  n-tuples  so  that  the  forbidden  n-tuple  can  occur  randomly  in 
every  out-of-sync  position  in  a  codeword  sequence.  Hence  the  sync 
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acquisition  delay  is  a  random  variable  with  sync  occuring  when  the 
forbidden  n-tupie  has  been  observed  in  all  but  one  possible  sync  position 
(the  correct  one).  When  messages  are  produced  independently,  an 
expression  for  the  average  sync  acquisition  delay  is  easily  computed. 

A  simple  synchronizer  for  this  scheme  exists. 

Another  method  of  designing  ss  codes  is  to  permit  a  fixed 
synchronizing  sequence  to  appear  in  only  one  position  relative  to  correct 
sync.  Although  this  restricts  the  maximal  dictionary  size  more,  synchron¬ 
ization  is  immediately  attained  upon  observing  the  synchronizing  sequence. 
When  the  messages  are  produced  by  a  memoryiess  source,  the  average 
sync  acquisition  delay  may  be  minimized  through  an  encoding  procedure 
which  is  also  developed.  For  a  given  synchronizing  sequence,  the  maximal 
dictionary  size  has  been  found,  and  the  average  sync  acquisition  delay 
minimized  a3  a  function  of  syr  hronizing  sequence  structure  and  position. 
These  two  minimization  techniques,  ,en  used  together,  provide  a  com¬ 
plete  minimization  of  average  sync  acquisition  delay. 

The  methodology  "»f  these  two  synchronization  techniques  is 
easily  extended  to  more  complex  systems  and  sources. 
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1 .  i .  4  Frequency  Tracking 

DA  -  A  R  0-D3  1  -  124-G92  9,  U.S.  Army  Research  Office 
R.  A.  Sc  Holtz,  W.  Case  ell 

The  problem  under  consideration  is  the  estimation  of  a  time 

varying  center  frequency  of  a  random  process.  The  analytic  representa- 

t  t 

/  ■  N,  N 

tion  of  the  process  is  x(t)  -  N  (i.  cos  (  .  aj(a)da  +  N  (t)  sin  uifalda 

c  v.!  ,/  s  v  y 

where  N  (t),  N  (t)  and  ui(t)  are  Gaussian  processes.  An  inteeral  equation 
c  s 

has  been  derived  for  the  maximum  a  posteriori  estimate  of  the  time  vary¬ 
ing  center  frequency  o> ( t ) .  This  equation  can  be  interpreted  in  terms  of  a 
feed-back  loop  with  non- realizable  filters.  A  linear  model  of  a  similar 
realizable  loop  has  been  derived  and  analyzed.  The  realizable  loop  har 
been  simulated  on  a  digital  computer  and  its  performance  obtained  in 
terms  of  the  normalized  error  variance.  The  linear  model  appears  useful 
in  determining  optimum  loop  parameters  and  in  predicting  trends  in  the 
mean  squared  error. 

3.2.5  Arithmetic  Norms  and  Arithmetic  Error-Detecting  and 
-Correcting  Codes 

DA -A  R  O-D  3  1  -  1 24 -G92  9,  Department  of  Army,  Army- 

Research  Office 

A.C.  L.  Chiang  and  I.  S.  Reed 

Any  integer  N  can  be  expressed  as  a  polynomial  in  radix-r: 
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n 

N  =  E  a.  rl  ,  0  <  |  a.  |  <  r  . 

.  „  1  —  1 

i  0 

The  smallest  number  of  non-zero  coefficients  a.  such  that  this  formula 

i 

holds  is  called  the  r-ary  norm  of  the  integer  N  and  is  denoted  by  ||  N II  . 

r 

In  particular,  if  r  =  2,  |j  N H ^  is  called  the  binary  norm  or  norm  of  the 
integer  N.  Properties  of  the  binary  norm  of  integers  have  been  extensively 
investigated*.  A  programmable  algorithm  for  the  computation  of  the 
binary  norm  of  an  arbitrary  integer  is  developed  and  a  table  of  norms  of 
the  natural  numbers  is  generated. 

Properties  of  the  binary  norm  of  integers  have  been  applied 

2 

to  the  development  of  arithmetic  error-detecting  and  correcting  codes  . 

In  order  to  enumerate  the  total  number  of  possible  error  patterns  of  an 
e-error  arithmetic  error-correcting  code,  one  needs  to  derive  the  number 
of  integers  with  norm  less  than  *e+l).  This  motivates  the  generation  of  a 
certain  abelian  group  which  we  call  a  norm  group.  The  number  of  equiva¬ 
lence  class  of  a  norm  group  is  derived  and  the  distribution  of  int  gers  of 
a  given  norrr  in  a  norm  group  is  computed.  F  rom  this  the  number  of 
possible  error  patternr  of  an  arithmetic  error-correcting  code  can  be 
found. 

LetX  be  a  set  of  consecutive  integers  (0,1,2,....,N-1‘  and 
Y  be  an  m-bit  register  vhich  stores  positive  integers  0  through  2  -1.  An 

arithmetic  error-correcting  code  is  defined  to  be  a  code  whose  coding 
function  0  maps  each  x  in  X  into  y  in  Y  such  that  0(x^  +  x.,)  =  0(Xj )  +  0(x^) 
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holds  for  all  x^,  and  x^+  x^  in  X.  This  definition  can  be  extended  to 

any  r-dirr,ensional  vector  space  Y  -  •  Y  ^ ,  Y^ . Y  .  In  such  a  case 

Y  can  be  considered  to  be  a  register  which  is  a  concatenation  of  all  disjoint 

subregister  Y.'s.  Now  let  the  coding  function  0(x)  -  (0^(x),0^(x) . .  (Mx)) 

be  a  mapping  from  X  into  Y  where  0. :  X  -*  Y  such  that  0.  (x  +  x  }  =  0.  (x  )  + 
rr  l  i  l  1  Z  l  1 

0  (x_)  holds  for  all  x,  ,  x.  and  x  +  x.  in  X.  We  are  particularly  interested 
l  Z  1  Z  1  Z 

in  a  set  of  residue  class  codes  for  which  0.  (x)  -  b,x  modulo  a.  ,  where  b. 

ii  i  l 

and  a  are  relatively  prime.  Emphasis  is  given  to  two  subclasses  of  such 
codes:  the  arithmetic  error-correcting  AN  code  for  which  0(x)  -  ax  and 
the  arithmetic  error-correcting  BN  modulo  A  code  for  which  0(x)  =  (x,  bx 
mod  a). 

We  have  found  bounds  on  the  arithmetic  error-correcting 
codes.  It  is  shown  that  the  upper  bound  on  the  single  error  correcting  AN 
code  can  be  achieved  if  and  only  if  the  constant  a  is  a  prime  and  Z  or  -Z 
is  a  primitive  root  modulo  a.  A  programmable  algorithm  for  computing 
the  best  bounds  on  the  higher  order  error-correcting  AN  codes  is  developed. 

An  arithmetic  error-correcting  BN  modulo  A  code  is  constructed 
in  such  a  way  that  there  is  no  interaction  between  its  information  symbol 
and  check  symbol.  A  necessary  and  sufficient  condition  for  the  construc¬ 
tion  of  a  single  error  correcting  BN  modulo  A  code  is  developed.  Some 
new  results  on  optimal  single  error  and  multiple  error  correcting  BN 
modulo  A  codes  have  been  found. 
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3.2.6  Phase  Locked  Loops  in  Cascade 

DA-ARO-D-3  1  -  124-G929,  U.S.  Army  Research  Office 

J.  J.  Stein  and  C.  I  .  Weber 

Approximations  to  the  joint  density  of  the  phase  errors  of 
two  second  order  phase  locked  loops  in  cascade  have  been  obtained. 

These  approxima^  solutions  to  the  corresponding  Fokker  -  Plane  k  equations 
have  been  shown  to  agree  very  favorably  with  empirical  data  at  low  as  well 
as  intermediate  signal-to-noise  ratios. 

In  addition,  the  approximate  solutions  have  been  shown  to 
compare  favorably  with  a  variety  of  computer  simulations  of  two  second 
order  phase-locked  loops  in  cascade. 

One  technique,  which  has  been  developed  exclusively  for  this 
appl  ication,  appears  to  be  sufficiently  general  to  be  classified  as  a  general 
method  for  simulating  stochastic  nonlinear  differential  equations  on  the 
digital  computer.  The  details  of  this  technique  arc  presently  being  pre¬ 
pared  in  report  form. 
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3.2,7  Signal  Design  for  Coherent  Sys terns 

AF-AFOSR- 1  555-68,  Air  Force  Office  of  Scientific  Research 
S.  Moskowitz,  D.  Sworder,  C.  Weber 

The  signal  design  problem  of  M-ary  coherent  communication 

systems  is  not  in  general  amenable  to  direct  analytical  solution.  Because 

of  the  form  of  the  problem,  gradient  ascent  techniques  suggest  themselves. 

The  gradient  of  the  detection  probability  with  respect  to  the  signal  vectors 

{S  '  can  be  approximated  using  Monte  Carlo  techniques  by  considering 

m 

the  probabilities  o!  detection  o!  a  set  of  M  perturbed  signal  sets  S' 

J ,  k  k  1 

who  re 

•  S  if  i  f  k 

J 

S'  -  - 

i ,  k  S  r  1  ill  k 

J 

The  vet  tor  *  is  such  that  S  *  1  satisfies  the  given  cn*  re\  restrictions 

j 

imposed  on  the  signal  set.  The  gradient  vn  tor  is  then  approximated  b\ 


-147- 


Pn(X,  {s.1-) 

D  J 


6  !i 


PD(X,{S,  })  -  PD(\,{SJ) 

!  p„(\,  f  s.  ‘  1  -  P_(X.  f  s.1) 

D  J .  m  D  j 


Initially,  a  study  was  made  of  sequentially  moving  each  signal 

in  the  signal  set  as  a  function  of  the  appropriate  component  of  the  gradient 

vector.  That  is,  if  (S.  i  is  the  original  signal  set,  then 

J  ° 


'V„+1 


S  }  +  K!  P_(\,  f  S  )  -  P_(\.  f  S  ' 

jn  ID  j ,  k  n  D  jn 


n  >  1 


where  k  =  n  mod  M.  An  important  outcome  of  thi  study  was  that  sequential 
movement  of  the  s  gnals  can  lead  to  false  optimum  points  although  conver¬ 
gence  to  a  true  optimal  point  could  be  attained  when  ■  S  1  was  properly 

j  o 

chosen. 


As  a  result  of  the  above,  the  technique  was  modified  to  calcu¬ 
late  the  total  gradient  vector  and  then  to  move  the  total  signal  set  in  the 
direction  of  the  ascending  path.  The  Mth  signal  vector  is  initially  held 
fixed  and  when  necessary  is  moved  in  the  radial  direction  so  that  the 
energy  constraints  of  the  signal  set  is  satisfied.  Employment  of  this 
technique  has  produced  good  results.  That  is,  for  equal  energy  signals 
restrict  .  to  two  dimensions  all  results  have  been  vcrilied  including  those 
for  unequal  a  priori  signals.  For  this  latter  case,  results  exist  only  for 
signal  sets  composed  of  three  signals.  Hence,  through  use  of  the  pro¬ 
posed  scheme  an  unlimited  number  of  new  results  can  be  attained. 
Furthermore,  for  eq  u-likely,  unity  equal  energy  signal  sets  confined  to 
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three  dimensions,  all  previous  existing  resims,  i.e,  ,  for  b  signals,  9 


signals,  and  4  signals,  have  been  verified  with  good  accuracy.  Additionally, 
the  gradient  approach  has  shown  that, for  7  signals  in  3-space, a  locally 
optimum  signal  s  »t  is  one  composed  of  9  equally  spaced  signals  on  the 
equator  and  2  signals  at  the  north  and  south  poles.  Analytically,  first 
order  condit.ons  for  this  signal  set  have  been  verified  at  all  SNR.  Com¬ 
puter  results  also  indicate  that  another  locally  optimal  set  has  one  signal 
placed  at  the  north  pole,  three  signals  equally  spaced  above  the  equator 
and  at  a  distance  of  1/3  from  the  origin,  and  three  equally  spaced  signals 

below  the  equator  shifted  bO  degrees  from  those  above  the  equator  and  2/3 

7 

of  the  distance  to  the  south  pole.  As  noted,  I  S.  -  0  in  both  of  the  above 

i- 1  ' 

cases. 

The  technique  also  applies  with  a  slight  modification  to  signal 

M  .  2 

sets  constrained  in  average  energy  by  ”  t>.  i|  S  i;  .  And,  through  the  use 

i  1  1 

of  generalized  spherical  coordinates  the  method  can  be  extended  to 
accommodate  signal  sets  confined  to  any  finite  dimensional  space. 
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3.2.8  Optimal  Physically  Realizable  Analog  Demodulation 

JPL  952108,  Jet  Propulsion  Laboratory 

C.  L.  Weber 

Nonlinear  recursive  filtering  techniques  are  being  applied  to 
the  problem  of  obtaining  optimal  analog  demodulation  under  a  variety  of 
modulation  methods ,  signal  assumptions  ,  and  channel  conditions.  The 
emphasis  is  being  placed  on  obtaining  physically  realizable  receivers, 
both  from  the  nonanticipative  point  of  view  and  the  system  complexity 
point  of  view.  Optimal  analog  demodulating  receivers  are  being  obtained 
for  phase,  frequency,  and  amplitude  modulating  transmitters.  Systems 
are  being  considered  in  which  the  received  signal  contains  unknown  param¬ 
eters  such  as  amplitude  (as  in  slow  fading),  frequency  offset  (Doppler 
effect),  and  unknown  phase  refeience. 

The  results  clearly  indicate  that  the  true  optimal  demodulator 
is,  in  most  cas*»s,  prohibitive  from  the  system  complexity  and  implementa¬ 
tion  standpoints.  However,  whenever  the  signal-to-noise  ratio  is  above 
some  nominal  threshold,  significant  approximations  can  be  made  which 
result  in  the  deletion  of  a  substantial  portion  of  the  system.  That  is,  the 
deleted  portion  has  only  second  >rder  effects.  The  end  structures  are 
therefore  approximately  optimal,  physically  realizable,  and  not  excessively 
complex. 

Preliminary  results  are  presently  being  prepared. 
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3.2.9  Generalized  Tracking  Theory 

AF-AFOSR‘496-67,  Joint  Services  Electronics  Program 
W.  C.  Lindsey 

Using  the  theory  of  Markov  processes,  the  response  probability 
density  function  p(_£,  t)  of  a  generalized  tracking  system  is  shown  to  be  the 
solution  to  a  (N+l)  -dimensional  Fokker-Pianck  equation.  The  vector 
=  (t,  y^.  .  .  ,  y  )  is  Markov  in  the  state  variables  and  $  is  the  system 
phase  error. 

It  is  further  shown  that  the  response  distributions  [p($,t), 
p(y  ,  t);  k  =  1 , 2 ,  .  .  .  ,  N ]  of  the  state  variables  satisfy  a  set  of  second- 

K. 

order  partial  differential  equations.  In  the  steady  state,  this  set  of  equa¬ 
tions  becomes  ordinary, first-order  differential  equations  tor  which  the 
exact  solutions  which  specify  the  marginal  probability  densities  [  p($), 
p(y  );  kr  1,2,...  ,  N]  ,  are  determined  by  two  sets  of  conditional  expec- 

K 

tations.  In  particular,  the  marginal  density  p($)  of  the  phase  error  is 
embedded  in  a  knowledge  of  the  conditional  expectations  f  E(y  |$), 
k  =  1,2,.  .  .  ,  N  1  while  the  ma rginal  dens ity  p(y  )  of  the  state  va  r  iable 

K 

y  is  embedded  in  a  knowledge  of  the  conditional  expectation  E(g(4)|y ,  ); 

K 

k  -  1,2,.  .  .  ,  N.  Here  g(^)  belongs  to  that  class  of  nonlinearitics  for 
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which  g(^)  is  an  odd  function.  The  conditional  expectations  are  approxi¬ 
mated  by  two  methods  and,  for  the  case  of  greatest  interest,  i.  e.  , 
g(^)  =  sin  cp ,  the  conditional  expectation  E(y^  j cp)  is  measured  via  com¬ 
puter  simulation  methods.  Agreement  of  the  simulation  results  with  those 
obtained  from  theoretical  considerations  are  within  less  .han  one  percent 
of  each  other. 

Synthesis  procedures  for  effecting  stochastic  optimization  of 
a  generalized  tracker  are  presented.  For  zero  detuning.it  is  shown  that 
the  tracker  which  minimizes  the  mean-squared  value  of  the  phase-error 
is  obtained  when  N  =  0  and  g(4)  is  proportional  to  sgn  (6).  In  the  case  of 
nonzero  detuning, it  i..  shown  that  the  loop  which  minimizes  the  mean-squared 
phase-error  is  second-order,  i,  e.  ,  N  =  l,  and  the  optimum  nonlinearity 
for  this  loop  is  also  proportional  to  sgn{£.)  For  large  loop  signal-to-noise 
ratios  and  equal  mean-squared  loop  phase-errors,  the  second-order  PLL 
system  requires  a  loop  signal-to-noise  ratio,  p,  of  approximately  p/2 
times  larger  than  that  required  in  the  optimum  second- order  tracker. 

At  low  loop  signal-to-noise  ratios  it  is  shown  that  ail  trackers  perform 
approximately  the  same.  For  second-order  trackers  some  results  are 
presented  to  the  case  where  g(^)  must  be  physically  realizable. 

The  field  of  probability  currents  is  derived  and  it  is  shown  that 
this  field  is  rotational  when  there  is  zero  detuning  in  the  loop.  From 
these  currents  the  average  number  of  phase-jumps  (cycles -slipped)  per 
unit  of  time  is  derived.  With  zero  detuning  this  average  is  zero. 
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3.2.  10  Phase  Tracking  in  Shot  Noise  Environments 

F  04701 -68-C-0234 ,  Space  Systems  Division,  Air  Force 
Systems  Command 
M.  Haney,  R.  M.  Gagliardi 

Timing  for  synchronization  and/or  ranging  in  an  optical 
communication  system  can  be  accomplished  by  intensity  modulating  with 
the  timing  signal  and  phase  locking  at  the  receiver  after  optical  demodu¬ 
lation.  When  direct  photodetection  is  used  for  demodulation,  however, 
the  demodulated  intensity  appears  as  a  shot  noise  current  signal  with  the 
desired  timing  signal  immersed  within. 

The  shot  noise  current  can  be  represented  by 

N  (t ) 

i(t)  r  hft-t  ) 

- 1  m 
m  -  1 

where  h (t - 1  )  is  the  current  waveform  due  to  a  single  electron  emitted 

m 

.  n  the  photodetector  at  time  t  .  N(t)  is  a  time  varying  Poisson  distri¬ 
buted  random  variable  such  that  the  probability  of  k  electrons  being 
emitted  in  the  time  interval  f  o.t1  is 

t 

r 

,  t  :  T]  n(T)  dr  , 

j  j  T]  n  (t  )  dr  ,  e  o 
L  -I 

P(N(t)  -  k)  =  - - - 
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where  7]  is  the  photodetector  emission  probability  and  n(t)  is  the  "timing 


signal  +  background  noise"  photon  particle  rate  incident  on  the  photo¬ 
detector. 

Let  the  current,  i(t),  be  the  input  to  a  phase  locked  loop  and 
let  the  photon  particle  rate,  n(t),  be 

n(t)  =  A(l+sin(u>t  +  0  {t)'')  +M 
V  o  1  y 


where  A  is  the  amplitude  of  a  sinusoidal  timing  signal  with  phase  0 ^  (t) , 
and  M  is  the  level  of  the  background  noise.  Also,  let  the  local  oscillator 
generate  a  signal  represented  by  cos  (tu  t  +  9  (t)).  The  phase  error, 

O  M 

4(t),  is  then 

6(t)  =  e^t)  -  ©2(t) 

The  nonlinear  stochastic  integro-differential  equation  describ 
ing  the  behavior  of  the  loop  becomes 


d4(t) 

dt 


dGj(t) 

~dT" 


t 

-  K  f(t-T)  r  h  ( T  -  T  )  cos  {  U)  t  +  0  (t)  dT  j 

m  V  o  2  y 

o  m  =  1 


where  K  is  the  loop  gain  constant  and  f (t - t )  is  the  impulse  response  of 
die  filter  within  the  loop. 

A  study  is  presently  being  conducted  to  determine  the  statistical 
behavior  of  the  phase  error  and  the  overall  loop  performance.  An  analysis 
involving  the  Chapman- Kolm ogorov  method  and  its  equation  in  series  form 
became  mathematically  unsound  because  the  higher  order  moments  of  the 
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incremental  phase  error  divided  by  the  time  increment  At  do  not  go  to 


zero  as  they  do  in  the  conventional  Gaussian  noise  phase  locked  loop.  In 
the  Gaussian  case  the  loop  statistics  are  governed  by  the  Fokker- Planck 
equation  involving  only  the  first  and  second  moments.  Other  methods  of 
loop  analysis  are  presently  being  investigated. 

3.2.11  Hybrid  Processing  of  Complex  Radar  Signals 
GK-3303,  National  Science  Foundation 
J.  E.  Ohlson 

Because  of  peak-power  limitations  in  radar  systems,  the 
pulse  energy  (which  determines  the  system  detection  threshold)  can  be 
increased  only  by  using  long  pulse  lengths.  Range  resolution  can  still 
be  the  same  as  that  of  a  very  short  pulse  if  the  time-bandwidth  product 
is  sufficiently  large.  One  way  of  producing  a  large  time- bandwidth  product 
is  to  linearly  increase  the  carrier  frequency  curing  the  pulse  length. 

This  type  of  signal,  known  as  a  "linear-FM"  or  "chirp"  waveform,  has 
been  used  for  st  eral  years  .  Generation  of  the  transmitted  signal  and 
processing  of  the  returned  signal  are  quite  difficult  procedures  since 
stability  'n  both  time  and  frequency  is  particularly  critical.  Virtually  all 
chirp  procedures  utilize  analog  techniques  and  are  subject  to  the  vagaries 
of  these  methods. 

In  this  research,  effort  is  being  made  to  utilize  digital 
methods  to  contribute  stability  and  precision  to  chirp  techniques.  For 
small  frequency  sweeps  (■  100  kHz)  direct  digital  synthesis  is  possible  , 
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but  for  large  sweeps  (  >  100  MHz)  direct  synthesis  cannot  be  performed 
because  of  state-of-the-art  limitations  on  computer  storage  and  ^--ra¬ 
tional  speed.  Effort  is  thus  being  devoted  to  developing  digitally  assisted 
signal  generation  and  processing.  The  problem  areas  are  (1)  repeatable 
signal  distortion  (e.g.  ,  nonlinearity  of  sweep),  and  (2)  fluctuations  of 
system  gain  and  phase  characteristics.  Work  to  date  has  produced  the 
following  preliminary  results: 

(1)  A  "sampling"  technique  can  be  employed  to  adaptively 
linearize  the  frequency  sweep.  This  will  remove  repeatable  distortion 
(i.e.  ,  constant  from  pulse  to  pulse). 

(2)  System  fluctuations  of  a  slow  nature  (e.g.  ,  temperature 
effects)  can  also  be  compensated  by  (1).  Fast  system  changes  appear  as 
a  noise,  which  because  of  the  large  bandwidths  involved  can  only  be  com¬ 
pensated  by  good  design  practice. 

(3)  A  problem  which  has  application  to  optical  (intensity)  data 
storage  has  been  examined.  It  was  shown  that  a  set  of  bounded,  orthonormai, 
non -negative  definite  basis  functions  cannot  exist  for  [ a, b^ .  The  result 

is  that  an  arbitrary  function  in  L^,  {  a,  b1  (in  particular,  the  non-negative 
subset  of  L^ia,  b1  )  cannot  be  represented  as  a  sum  of  positively  weighted 
non-negative  orthonormal  basis  functions.  A  method  for  constructing  a 
bounded,  linearly  independent  set  of  non-negative  basis  functions  from  a 
tr igonometric  fourier  series  is  found,  which  then  allows  the  desired 
representation  to  be  found  under  somewhat  weaker  conditions. 

(4)  When  the  output  of  a  causal  filter  is  to  be  sampled  for 
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digital  processing,  the  Nyquist  sampling  rate  of  twice  the  bandwidth  is 
undefined  since  the  signal  spectrum  is  not  bandlimited.  A  higher  sampling 
rate  always  produces  greater  fidelity,  but  at  greater  expense.  Take  the 
case  of  a  constant- valued  signal  immersed  in  white  noise  which  is  then 
filtered  by  a  single  pole  filter,  sampled, and  the  samples  summed  to  esti¬ 
mate  the  signal  level.  The  variance  of  this  estimate  is  found  as  a  function 
of  sampling  rate  and  is  compared  to  that  of  the  estimate  obtained  by  a 
finite-time-average  of  the  signal  and  noise  over  the  same  time  interval. 
The  exact  solution  has  also  been  found  for  a  second  order  filter.  For  a 
particular  filter,  the  sampling  rate  necessary  for  a  desired  efficiency  is 
thus  found. 
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Research  Office 

S.  W.  Golomb,  D.  Sherman,  C.  Fuzak,  H.  Trachtenberg 

3.2.  12.  1  A  Model  for  Voicing  Duration 
D.  Sherman 

I.  The  Model 

In  constructing  a  model  for  voicing  duration  based  upon 
phoneme  durations,  two  extremes  are  possible.  The  first  extreme  involves 
simplifying  as  much  as  possible  the  phoneme  duration  distribution.  Essen¬ 
tially,  this  means  assuming  that  all  voiced  phonemes,  in  all  positions  in 
the  speech  stream  (i.e.  ,  all  phonemic  environments),  are  drawn  from  a 
common  population  in  terms  of  duration.  Statistical  tests  (see  part  III) 
indicate  that  so  simple  an  assumption  is  unjustified. 

The  other  extreme  in  modeling  involves  examining  the  duration 
of  each  phoneme  class  individually,  as  a  function  of  phonemic  environment. 
This,  of  court**,  is  an  extremely  tedious  approach,  perhaps  interesting  in 
itself,  but  does  not  necessarily  yield  the  most  useful  voicing  duration 
model.  The  fine  detail  involved  in  a  model  constructed  in  this  manner 
would  be  lost  in  practical  tests  and  applications. 

Asa  compromise  between  these  extremes  a  voiced  segment 
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was  considered  to  be  composed  of  four  types  of  phonemes -- thos e  which 


initiated  voicing,  those  which  terminated  voicing,  those  which  neither 
initiated  nor  terminated  voicing,  and  those  which  did  both.  These  classes 
were  called  Initial,  Terminal,  Medial,  and  Isolated,  respectively.  It  was 
felt  that  these  classes  represented  the  basics  of  phonemic  environment, 
while  ignoring  differences  due  to  the  phonemes  themselves.  If  later  tests 
indicated  need  for  further  refinement,  the  phonemic  differences  might  be 
considered  within  the  classes. 

Given  these  ground  rules,  the  model  for  voicing  duration 
follows  immediately. 

Let:  P.  =  Probability  of  i  phonemes  in  a  voiced  segment. 

Pj ( t )  -  Duration  density  for  an  Initial  voiced  phoneme, 

p  (t)  -  Duration  density  for  a  Medial  voiced  phoneme. 

Pj(0  -  Duration  density  for  a  Terminal  voiced  phoneme, 

p  <0  -  Duration  density  for  an  Isolated  voiced  phoneme. 

p^(t)  -  Duration  density  for  a  voiced  segment. 

Then  the  duration  density  for  a  voiced  segment  is  given  by 


D 


P  p  ,  +  P  (p  *  p  l  +  P  (p  $  p  c-  p  ) 
PS  2  M  HT  3'M  hM  fT 


*  p  (P.  *  P 

n  7 


M 


p  *  p  )  t 

*  M  *T 


—  n  -2  fold  — ■ 

where  the  densities,  p,  are  understood  to  be  functions  of  t.  In  order  to 
make  use  of  this  expression,  each  of  the  terms  on  the  right  must  be 
determined  up  to  a  set  of  parameters.  These  parameters  would  then  be 
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estimated  from  the  voicing  duration  data,  and  fit  tests  performed.  If  a 
good  fit  results,  it  is  then  expected  that  the  parameters  should  prove 
useful  in  determining  such  factoru  as  language,  speaker,  sex,  or  speech 

rate. 

Note  that  in  deriving  this  model  it  was  assumed  that  the 
duration  of  consecutive  phonemes  were  independent;  that  is,  the  duration 
is  a  function  only  of  the  class  to  which  the  phoneme  belongs. 

II.  The  Data 

A  two  minute  continuous  speech  sample  was  obtained,  and 
attempts  were  made  to  phonemically  segment  it  on  the  basis  of  the  ampli¬ 
tude  waveform.  These  attempts  proved  unsuccessful,  so  spectrograms  of 
the  sample  were  obtained  utilizing  the  Kay  Electric  type  B/65  sonagrams 
produced  by  the  USC  Asian  Studies  Department.  The  resulting  phoneme 
duraf'on  data  is  shown  in  Table  I.  The  entries  in  the  table  are  measure¬ 
ment  units,  equal  to  3  1/8  milliseconds.  Casual  inspection  indicates  that 
the  Isolated  phonemes  appear  to  come  from  a  different  population  than  the 
other  three  categories.  The  remaining  three  categories  were  tested  for 
a  possible  common  population. 

The  tests  oerformed  were  of  the  fit  tyoe.  In  ail  cases,  a 
significance  level  of  10%  (C.  r  .10)  was  selected.  This  allowed  a  proba¬ 
bility  of  .  1  of  a  Type  I  error  (rejection  of  the  hypothesis  when  it  is 
actually  true).  This  fairly  large  value  for  1  was  used  in  an  attempt  to 
assure  that  the  resulting  model  would  be  as  accurate  as  possible  based 
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upon  the  available  data.  Hereafter,  rather  than  simply  reporting 
acceptance  or  failure  at  the  10%  level,  the  minimum  significance  level 
for  acceptance  is  reported,  which  will  give  an  estimate  of  "quality"  of 
acceptance  or  rejection, 

A  three  sample  Kruskal- Wallis  H  test  on  the  Initial-Medial- 
Terminal  data  yielded  a  value  of  H  =  14.  32.  For  this  value,  the  hypothesis 
that  the  three  samples  are  drawn  from  a  common  population  is  rejected 
for  any  significance  level  larger  than  .001,  and  we  may  safely  conclude 
that  more  than  one  population  is  involved.  This  test  was  performed  on 
the  data  auer  excluding  voiced  stops  and  fricatives  in  the  Medial  and 
Terminal  positions.  (Stops  and  fricatives  were  not  treated  separately 
in  the  Initial  case  due  to  their  low  frequency  of  occurrence. ) 

Pairwise  (Mann- Whitney)  tests  on  the  same  data  gave  the 
following  results: 

Initial-Medial  X  =  1.70,  corresponding  to  a  ■  .09 

Initial-Terminal  X  =  1.47,  corresponding  to  a  =  .14 
Medial-Terminal  X  =  4.02,  corresponding  to  a  -  .001. 

Since, at  most, one  hypothesis  may  be  accepted,  we  would  accept  the  one 
which  would  be  accepted  for  the  largest  val.  e  of  a  .  At  this  point  we  might 
reasonably  accept  the  hypothesis  that  the  Initial  and  Terminal  samples 
are  drawn  f-om  populations  with  the  same  distribution  functions.  This 
decision  is  postponed,  however,  in  favor  of  further  tests. 

Earlier  investigations  indicated  that  vowel  durations  might 
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be  well  represented  by  a  log-normal  density;  that  is, 


P(t) 
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It  was  decided  to  investigate  this  possibility,  so  chi-squared  goodness  of 
fit  tests  were  performed  using  maximum -likelihood  estimates  of  p  and  o. 
Again,  a  10%  (a  =  .  10)  significance  level  was  established.  The  results  are 
presented  in  Table  II.  Since  the  fit  was  extremely  poor  in  the  case  of 
the  Terminal  data,  fit  tests  were  also  performed  against  a  normal  density 
function.  As  a  result  of  the  tests  we  can  conclude  that  each  of  the  following 
phoneme  categories  may  be  represented  by  log-normal  duration  density 
functions:  Initial;  Medial  '  vowels";  Medial  Stops  &  Fricatives;  Isolated. 

In  addition,  we  may  conclude  that  the  best  fit  in  the  Terminal  case  occurs 
for  total  Terminal  data  vs.  a  log-normal  density.  Thus  the  best  available 
reasonable  assumption  is  that  the  Terminal  data  is  drawn  from  a  popula¬ 
tion  with  a  log-normal  density.  The  fit  in  the  Terminal  case  was  expected 
to  be  the  least  accurate,  due  to  the  voiced-unvoiced  decision  process  used 
on  the  spectrograms.  (The  decision  was  based  upon  the  low  frequency 
power,  a  value  judgement  based  on  the  degree  of  blackness  of  the  spectro¬ 
gram.  The  power  onset  time  is  much  less  than  the  decay  time,  and  timing 
errors  tend  to  accumulate  in  the  terminal  position.  )  Due  to  the  tenative 
nature  of  the  acceptance  of  log-normal  density  for  the  Terminal  phonemes, 
we  will  reject  the  hypothesis  that  the  Initial  and  Terminal  data  come  from 
populations  with  the  same  densities. 
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Chi-Squared  Goodness  of  Fit  Tests 


2 

Number  of 

Significance 

Phonemes 

X 

Regions 

Level 

Fit  to  Log- 

■Normal  Density 

Initial  (total) 

4.  07 

15 

.  98 

Medial  "vowels" 

14.  79 

20 

,20 

Medial  stops  &  frictatives 

4.29 

6 

.25 

Terminal  "vowels" 

25.24 

14 

.  007 

Terminal  (total) 

23.  79 

16 

.04 

Isolated 

3.20 

5 

.20 

Fit  to  Normal  Density 

Terminal  "vowels" 

33.38 

14 

.  005 

Terminal  (total) 

2  5.26 

16 

.  02 

Table  II 


Note:  "vowels  "  are  vowels  , /r/ ,  /!/,  glides,  nasals,  and  diphthongs. 


Table  III  illustrates  the  measured  frequency  of  occurrence 
for  various  numbers  of  phonemes  per  voiced  segment.  Due  to  the  long 
"tail",  the  fit  to  Posson  is  quite  poor,  failing  the  a  =  .  1  test.  The  log¬ 
normal  is  suggested  by  the  fact  that  the  number  of  letters  per  word  and 

the  number  of  phonemes  per  word  both  appear  to  be  log-normal  for 

2 

English.  Maximum  likelihood  estimates  give  p  -  .881,  0  =  .422,  and  a 

2 

value  of  X  (for  a  chi-squared  fit  test)  of  .8698  for  6  regions.  This 
passes  the  test  for  Cl  =  .10,  giving  in  fact  a  significance  level  of  roughly 
80%.  Thus,  in  the  model  for  p^ ,  we  will  use  a  log-normal  expression  to 
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give  the 


P.  's. 


The  dec  is  ion  t 


categories  means  that  p. 


M 


M 


o  treat  the  medial  phonemes  in  two  separate 
in  the  expression  for  p^  will  have  the  form 

-  P  p  +  P  p 
V  fMV  S  hMS 


The  ratio  of  occurrence  of  the  "vowel"  group  (v)  to  the  stop-fricative 
group  (s)  is  theoretically  about  3.  13.  The  measured  frequencies  (103/20 
from  Table  ILA)  yield  about  3.  56.  Accepting  the  theoretical  values  as 
approximately  correct,  we  have  P  =  .76  and  P  =  .24.  These  values 

V  o 

will  be  used  in  the  model. 


#  of  Phonemes  in 
Voiced  Segment 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


#  of 

Occu  ences 

25 

34 

18 

13 

6 

3 

0 

3 

2 

1 

1 

1 


Table  III 
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III.  Continuation 


The  next  step  in  constructing  the  model  is  a  determination  of 

the  exact  form  of  the  P.'s.  Once  this  is  accomplished,  maximum  likeli- 

1 

hood  estimators  must  be  found  for  the  twelve  parameters  in  the  model. 

(These  parameters  are  M  and  0  for  p  ,  p  ,  p  ,  p  ,  p  and  for  the 

I  MV  Mo  T  S 

expression  for  the  P  's. 

The  model  must  then  be  tested  against  duration  data,  and  if 
fit  tests  prove  successful,  then  experiments  must  be  designed  to  test  the 
parameters  for  dependencies  on  factors  such  as  speaker,  language,  sex, 
and  speech  rate. 
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3.3  SWITCHING,  AUTOMATA  THEORY  ,  AND  COMPUTERS 


3,  3.  1  Minimum  Length  Race-Free  Coiling  lor  Asynchronous 

Sequential  Machines 

AF -A FOSR  -  1 0  i  8 -6 7A  ,  Air  Force  Office  of  Scientific 

Research 

R.  B.  McGhee  and  R.  M.  Kashef 

The  preliminary  work  reported  in  previous  progress  reports 
have  been  gene ral ized  for  augmented  parity  check  (APC)  codes  and 
extended  for  flow  table  columns  whose  states  can  be  partitioned  into 
maximum  distance  pairs  and  published  in  the  form  of  R.  M.  Kashef's" 

Ph.  D.  Dissertation  .  An  aostract  of  the  dissertation  follows; 

One  of  the  basic  problems  in  the  design  of  asynchronous 
sequential  machines  is  the  determination  of  an  assignment  of  binary  codes 
to  the  states  of  a  machine  in  a  wav  which  avoids  critical  race  conditions. 
Solutions  to  this  problem  can  be  found  either  by  means  of  a  coding  algor¬ 
ithm  or  by  the  use  of  general  codes  capable  of  encoding  any  machine. 

This  dissertation  is  concerned  with  the  latter  approach  and  is  aimed  at 
the  specific  problem  of  obtaining  a  race-free  encoding  with  a  minimum 
number  of  bits. 

Up  to  the  present  tune,  the  shortest  length  general  race-free 
code  known  has  been  the  intermeshed  row  set  code  devised  by  D.  A, 

❖ 

Middle  initial  of  F.  M.  Kashef  has  recently  been  changed  to  S  from  M. 
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3 

Huffman  .  This  code  encodes  a  table  of  2  rows  with  2S  -1  bits.  In 

o 

this  dissertation,  a  new  family  of  codes,  called  augmented  parity  check 

S 

codes  (APC  codes),  is  introduced.  An  APC  code,  S,  is  a  set  of  2  n-bit 
code  words  such  that  any  element,  s ,  of  S  contains  an  even  (or  odd) 
number  of  ones  in  SQ  specified  bit  positions,  and  all  but  one  of  the  remain¬ 
ing  bits  of  s  are  equal  to  a  constant  value,  if  each  element  of  an  APC 

S0 

code  is  aseigned  to  one  of  the  states  of  a  2  -row  state  table,  then  the 
n  So 

remaining  2  -2  binary  n-tuples  are  available  for  connecting  the  st  e 

codes  as  required  by  the  state  table.  The  main  problem  to  be  solved  is 

thus  the  determination  of  a  minimum  value  for  n  as  a  function  of  S  . 

o 

A  lower  bound  on  n  can  be  obtained  by  considering  a  column 

S0  So-l 

of  a  2  -row  state  table  which  can  be  partitionc  1  into  2  pairs  of  states 
such  that  number  of  bits  in  which  the  corresponding  code  words  differ  is 
the  maximum  possible  for  any  two  code  words  of  S.  In  this  dissertation, 
it  is  shown  that  it  is  impossible  to  code  such  a  column  so  that  every  pair 
of  states  of  the  column  to  be  connected  is  joined  by  disjoint  race-free 


chains  of  n-tuples  with  less  than  S  +  Log  (  < 

o  2  V 


V  1 


2  >  +  1  j  >  bits , 


where  <  x  >  denotes  the  next  integer  greater  than  or  equal  to  x  and 

S  >  4.  In  addition  to  this  result,  it  is  also  shown  that  no  more 
o  - 

S  $  N 

2  <  —  >  +  <  Log^^  <  —  >  >  bits  are  ever  required  to  encode  such 

a  maximum  distance  partition.  For  large  values  of  S^,  both  of  the  above 

numbers  are  substantially  less  than  2 S  -  1. 

o 


By  the  use  of  the  concept  of  equivalence  classes  of  mappings 
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of  a  finite  set  S  into  a  finite  set  R  =  [0,1],  under  a  permutation  group 

G  on  S,  it  was  possible  to  reduce  the  number  equivalence  classes  of 

stable  state  configurations  for  =  3,  n  =  4,  to  13  by  letting  S  be  a  4- 

bit  APC  code  and  G  be  the  group  of  distance  preserving  transformations 
4 

on  S  .  Furthermore,  it  was  possible  to  code  a  canonical  member  of 

each  equivalence  class  directly,  hence  implying  that  any  4-bit  APC  code 

can  be  used  to  obtain  a  race-free  encoding  of  all  8  state  flow  tables.  The 

above  concept  was  used  for  state  tables  with  S  =  4  and  n  =  6  and  the  number 

o 

of  equivalence  classes  of  stable  state  configurations  reduced  to  271.  By 

5 

the  use  of  a  computer  program  based  on  backtrack  programming  ,  every 
possible  type  of  columns  of  16- row  state  tables  were  coded  with  the 
exception  that  for  about  one  percent  of  the  columns  with  5,  6,  and  7  stable 
states  the  coding  was  not  attempted.  There  is  no  evidence  to  indicate  that 
the  uncoded  columns  could  not  be  coded  by  a  further  expenditure  of  com¬ 
puter  time.  An  isomorphism  between  the  group  of  distance  preserving 
transformations  ^n  the  set  of  APC  codes  and  the  group  of  the  symmetries 

of  the  S  -cube  was  established  for  S  -odd,  and  it  was  therefore  concluded 
o  o 

that  it  is  not  practical  to  carry  out  the  proof  for  Sq  >  5,  using  the  method 
for  S  :  4.  The  result  for  Sq-  4  encourages  the  conjecture  that  the  bounds 
obtained  for  maximum  distance  partitions  are  valid  for  any  type  of  flow 
table  columns.  Further  research  is  required  to  determine  whether  or  not 
this  is  actually  the  case. 

Presently,  a  modified  computer  program  has  been  written 
in  order  that  the  remaining  one  percent  of  1 6-row  tables  with  5,  6  and  7 


-169- 


stable  states  can  be  coded.  It  is  anticipated  that  the  program  run  time 
will  exceed  100  hours  on  the  IBM  360,  Model  44,  and  methods  are  being 
investigated  to  reduce  this  time.  Furthermore,  since  this  research 
indicates  that  the  number  of  bits  required  to  code  the  states  of  an  asyn¬ 
chronous  sequential  machine  using  A  PC  codes  is  the  same  as  the  number 
of  bits  for  a  synchronous  sequential  machine  with  the  same  number  of 
states,  plus  the  bits  required  for  the  clock,  then  the  race-free  coding 
problem  of  an  asynchronous  sequential  machine  would  be  completely  solved, 
if  it  would  be  possible  to  prove  that  the  columns  which  can  be  partitioned 

into  maximum  distance  pairs  are  the  only  columns  which  require  the  use 
n-S„ 

o  n 

of  all  2  unused  code  words  among  the  set  of  all  possible  2  n-bit  code 

S0 

words  for  a  2  -row  state  table,  using  A  PC  codes.  Hence,  further  research 
is  under  way  to  investigate  the  methods  by  which  the  above  can  be  proven. 
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3.3.2  Interactive  Hybrid  Computer  Systems 
GK  2716,  National  Science  Foundation 
M.  J.  Merritt,  G.  A.  Bekey 


1 .  Background 

When  an  analog  computer  and  a  digital  computer  are  tied 
together  by  suitable  linkage  devices,  the  result  is  called  a  hybrid  com¬ 
puter.  The  diverse  nature  of  the  two  basic  computing  devices  makes  it 
difficult  for  a  single  individual  to  master  all  aspects  of  both  analog  and 
digital  computation.  Recent  developments  in  computer  generated  graphic 
displays  have  made  possible  the  development  of  software  systems  to  supply 
the  needed  expertise.  In  addition,  the  graphic  display  devices  allow  the 
user  to  communicate  with  the  hybrid  computer  in  a  single,  task  oriented 
language.  The  resultant  improvement  in  the  man-hybrid  computer  inter¬ 
face  is  necessitated  bv  the  increasing  scope  of  hybrid  computation. 

2,  Research  Objectives 

The  research  objectives  of  this  study  alt  involve  the  improve¬ 
ment  inman-hybrid  computer  communication.  Graphic  sketches,  line 
drawings,  and  schematic  diagrams  are  the  natural  forms  of  expression 
used  by  engineers  and  scientists.  Studies  are  underlay  to  determine 
the  most  eifect’vc  graphic  language  format  for  hybrid  computation. 

Existing,  nongraphic,  simulation  languages  provide  a  base  for  these 
s tudies . 

Additional  a  c  r .  is  m  progress  to  develop  graphically  or-ented 
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application  and  support  programs  to  enhance  existing  hybrid  computer 


software  systems.  M’.r  may  be  put  back  into  the  parameter  adjustment, 
optimization  loop  by  providing  him  with  suitable  graphic  displays  and  on¬ 
line  controls.  A  large  class  of  search  procedures,  trial,  and  error 
procedures,  and  iterative  computational  procedures  may  be  adopted  to 
allow  on-line  supervision  and  control  by  a  human  operator  observing 
graphic  presentations. 

3.3.3  Realization  of  Fault  Detection  Sequences  in  Linear  Cascades 
of  Identical  Machines 

AF-AFOSR  -  196-67  ,  Joint  Services  Electronics  Program 
M.  A.  Breuer 

Our  research  deals  with  the  problem  of  constructing  fault 
detection  tests  for  sequential  machines.  This  is  currently  one  of  the 
most  significant  problems  related  to  the  successful  acceptance  of  the 
new  LSI  technology.  We  consider  a  specific  machine  model,  namely  a 
linear  cascade  of  identical  sequential  machines,  as  illustrated  beiow. 
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M  is  called  a  cascade  machine  and  M  M  M  ...  M  . 
i  1  2  N 

Examples  of  linear  cascade  machines  are  sh-ft  registers  and  parallel 
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adders.  Our  main  interest  is  to  find  test  sequences  for  these  machines 

which  are  much  shorter  than  those  produced  by  the  classical  techniques  of 

1  2 

Hennie  and  Kohavi  and  Lavallee  . 

The  following  prop  rties  were  investigated,  and  results 

obtained: 

a)  Wnen  are  two  linear  cascades  of  N  identical  machines 

equivalent  ? 

b)  What  is  the  reduced  structure  of  a  linear  cascade  of  N 
identical  machines  ? 

c)  What  information  lossless  properties  of  cascade  machines 
are  sufficient  in  order  to  insure  testability  of  the  linear  cascade? 

d)  The  concept  of  a  complete  sequence  generation  is  intro¬ 
duced,  and  its  influence  on  testability  examined. 

(Briefly,  M  is  a  complete  sequence  generator  if, starting  in  the  initial 
state,  there  exists  input  sequences  which  will  produce  any  desired 
sequence  on  the  y'  output.  ) 

e)  If  M  is  information  lossless  of  finite  order  and  is  also  a 

1 

complete  sequence  generation,  then  we  show  how  to  test  the  linear  cascade 
of  N  identical  machines  with  a  sequence  of  N  uniform  expe  i  hnents . 

f)  Ultra  efficient  testing  procedures  are  studied,  which  allow 

for  many  of  the  M  's  to  be  tested  simultaneously. 
i 

g)  Sufficient  conditions  are  derived  on  M  m  order  to  insure 
that  either  every  second  machine  or  every  third  machine  can  be  tested 
simultaneously. 
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As  one  practical  example  of  this  work,  we  have  constructed  ] 

i 

! 

a  sequence  of  length  16  which  will  'test  an  N  bit  parallel  binary  adder.  '  j 

i 

This  sequence  is  independent  of  N.  This  work  will  be  presented  at  the 
1968  IEEE  Conference  on  Switching  and  Automata  Theory, 

' 

Work  on  generating  a  test  sequence  of  the  general  model  of  j 

a  sequential  machine  is  continuing.  , 

j 

i 
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3.3.4  Abstract  Pattern  Recognition 

NGR-05-018-044,  National  Aeronautics  and  Space 

Administration 

W.  S.  Meisel,  W.  Yuan 

In  abstract  pattern  recognition  one  is  presented  with  samples 
(in  a  finite-dimensional  vector  space)  of  varic’s  pattern  classes.  The 
problem  is  to  classify  other  points  in  the  space  on  the  basis  of  this  limited  ’ 
information.  One  might  describe  the  motivation  for  such  studies  as  a 
desire  to  identify  the  contents  of  a  "black  box"  which  does  not  contain 
memory.  Obvious  applications  lie  in  engineering,  economics,  and 
biomathematics . 
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The  problem  has  two  aspects:  feature  selection  and  pattern 
class  if  ication. 

Three  basic  approaches  have  been  applied  to  the  latter: 

a)  Nonlinear  parameterization  -  The-  discriminant  function  is 
assumed  to  take  a  certain  form  specified  completely  by  a  finite  set  of 
parameters;  the  parameters  appear  nonlineariy.  A  common  example  is 
the  assumption  that  the  distribution  is  Gaussian. 

b)  Linear  parameterization  -  The  discriminant  is  assumed  to 
be  of  the  form 


R 

P(x)  =  E  C.0.  (x)  , 


where  the  0  are  given  functions,  and  the  C.  paremters  to  be  determined. 
i  b  i 

This  expression  may  represent,  for  example,  a  general  multivariate  poly¬ 
nomial. 

c)  Composition  -  The  discriminant  is  assumed  to  take  the 

form 


P.(x)  =  — - —  T.  Mx,y)  . 

1  |  S  |  y  e  S. 

i  i  i  i  j 

where  S.  is  the  set  of  samples  of  class  i.  This  latter  equation  may  be 
viewed  as  the  superposition  of  contributions  from  each  of  the  sample 
points  of  class  i. 

Our  work  has  concentrated  on  methods  of  type  (b)  and  (c). 
The  choice  of  linear  parameters  by  least-square  or  least-mean- square 
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approximation  to  the  unknown  probability  function  has  resulted  in  some 

practical  algorithm^.  A  study  of  potential  functions  has  led  to  an  algorithm 

2 

with  distinct  advantages  over  similar  algorithms  . 

We  have  generated  artificial  data  to  test  and  compare  tech- 

.  •  3,4 

niques  f. or  abstract  pattern  recognition 
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3.3.  5  Automata  and  Formal  Language  Theory 

AF-AFOSR-496-67,  Joint  Services  Electronics  Program 
R.  Card,  W.  Chandler,  S.  Ginsburg 


Background 

During  the  past  five  years,  a  multitude  of  different  types  of 

acceptors  and  the  languages  they  recogn.ze  have  been  introduced. 

1 

Recently,  it  was  observed  that  most  of  these  devices  and  languages 
would  be  subsumed  and  unified  within  the  notions  of  an  "abstract  family 
of  acceptors",  abbreviated  AFA ,  and  "abstract  family  of  languages", 
abbreviated  AFL.  This  has  resulted  in  a  flood  of  new  questions  and 
results,  hitherto  unsuspected.  Whereas  before  the  advent  of  AFA,  this 
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was  a  surfeit  of  devices;  after  AFA  there  was  a  sudden  scarcity.  This 
situation  arose  because  there  were  not  enough  different  types  of  extant 
devices  to  answer  the  new  questions.  A  typical  new  question  which  is 
receiving  attention  is:  Find  types  of  AFA  whose  associated  family  of 
languages  is  not  closed  under  substitution. 

Progress 

Roger  Card  is  studying  the  properties  of  "associative  memory 
acceptor"  languages.  In  particular,  he  is  examining  their  elementary 
closure  properties  and  relations  to  well-known  families  of  languages. 

William  Chandler  is  studying  languages  accepted  by  one-way 
deterministic  acceptors.  He  has  obtained  a  characterization  of  these 
languages  in  terms  of  closure  operations  and  is  now  examining  their 
elementary  properties. 

A  paper,  "Substitution  in  Families  of  Languages",  was  written 
and  submitted  for  publication.  In  it,  the  role  of  substitution  to  AFJ„  is 
studied.  For  example,  it  is  shown  that  if  and  L^  are  AFL,  then  so  is 
the  set  of  all  languages  obtained  by  substituting  languages  of  into  those 

°‘L2- 

For  a  number  of  years,  linguists  have  been  searching  for 
faithful  models  of  classes  of  natural  languages.  Among  the  models  pro¬ 
posed,  the  one  receiving  the  most  attention  these  days  is  that  of  a 
"transformational  grammar".  As  yet,  no  mathematically  sound  model  of 
these  grammars  has  been  formulated  which  has  received  the  blessings  of 
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a  large  number  of  linguists.  In  a  report,  "A  Mathematical  Model  of 


Transformational  Grammars",  submitted  for  publication,  an  attempt  is 
made  to  provide  such  a  mathematical  model. 

Two  papers,  items  (3)  and  (4),  were  published. 
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4.  BIOMEDICAL  ENGINEERING  AND  MATHEMATICS 


4.  1  CARDIOVASCULAR  AND  RESPIRATORY  SYSTEMS 

4.1.1  Regulations  of  Pulmonary  Ventilation 

He- 11  915-01,  National  Institutes  of  Health 
NGR-05-018-044 ,  National  Aeronautics  and  Space 
Administration 

S.  Yamashiro,  M.  B.  Wolf  and  F.  S.  Grodir.s 

The  goal  of  the  current  laboratory  investigation  is  to  determine 
the  role  of  the  central  chemosensitive  area  in  the  respiratory  response  to 
acute  metaboloic  disturbances  of  acid-base  balance.  We  plan  to  do  this 
by  simultaneously  monitoring  respiration  and  the  acid-base  state  of 
arterial  blood,  brain  venous  blood,  and  CSF  (cerebrospinal  fluid)  following 
acute  metabolic  disturbances.  Some  preliminary  results  have  been 
obtained  fruin  experiments  in  which  respiration,  arterial  blood,  and  CSF 
were  monitored  following  infusion  of  bicarbonate  into  anesthetized  dogs. 

Acute  metabolic  alkalosis  has  been  induced  in  20  anesthetized 
mongrel  dogs  by  injecting  (.3M-1M)  NaHCO^  intravenously  via  a  femoral 
vein.  Doses  of  iO  mEq/Kg  were  used  and  were  infused  at  a  rate  of  20 
ml/min.  Respiration  was  monitored  by  a  spriometer  filled  with.  100% 
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oxygen  and  connected  to  the  animal  via  a  tracheal  cannula.  Arterial 

blood  samples  were  drawn  from  a  femoral  artery  before  infusion  and 

periodically  following  infusion.  Blood  pH  and  Pco^  ware  directly 

measured  with  appropriate  electrodes  (Corning  Model  16  Blood  Analyzer). 

In  5  dogs,  CSF  was  sampled  via  a  spinal  needle  placed  into  the  cisterna 

magna,  and  its  pH  and  Pco?  were  determined  in  the  same  way  as  in  blood. 

In  general,  our  results  agreed  with  the  findings  of  Robin  et  al*.  There  was 

a  consistent  rise  in  ventilation  following  infusion  of  bicarbonate  in  all 

experiments  not  involving  a  gross  technical  error.  In  3  of  5  dogs  in 

which  CSF  was  sampled,  the  average  pH  change  in  CSF  immediately 

following  bicarbonate  infusion  was  -.008  of  a  pH  unit.  Results  from  th° 

other  two  dogs  were  discarded  because  of  gross  technical  errors.  In  all 

three  dogs  a  net  alkaline  shift  of  CSF  pH  was  seen  on'’  hour  after  infusion. 

Robin  et  al.  *  report  a  depression  of  ventilation  following  bicarbonate 

infusion  at  exactly  the  same  time.  These  results  fully  support  the  impor- 

2 

tance  assigned  to  CSF  pH  by  Pappenheimer  in  the  control  of  ventilation 
during  chronic  metabolic  disturbances.  The  present  results  suggest  that 
CSF  pH  is  also  important  in  acute  disturbances. 
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4.1.2  Blood  Pressure  and  Heart  Rate  Regulation  in  the  Fetus 

GMO1724-02,  National  Institutes  of  Health 

W.  Morrison  and  G.  A.  Bekey 

This  project  has  been  concerned  with  both  experimental  and 
computer  studies  of  the  fetal  cardiovascular  system.  During  the  report 
period,  the  emphasis  of  the  project  was  on  computer  studies. 

The  fetal  cardiovascular  system  was  simulated  both  with  and 
without  blood  pressure  control  (the  so-called  "baroceptor  reflex").  Model 
jen3itivi  y  studies  were  carried  out  to  determine  the  influence  of  various 
parameters  on  the  agreement  between  model  behavior  and  experimental 
behavior.  Some  improvement  in  the  fit  of  model  responses  to  experimental 
data  was  obtained  when  certain  nonlinearities  were  added  to  the  model, 
namely: 

(a)  the  resistance  of  certain  vessels  to  blood  flow  was 
assumed  to  be  a  nonlinear  function  of  vessel  diameter. 

(b)  the  elastance  of  the  ascending  aorta  and  pulmonary 
artery  were  made  piecewise  constant  functions  of  the 
vessel  diameter. 

The  final  portion  of  the  computer  study  was  concerned  with 
mathematical  models  of  the  baroreceptor  reflex.  Optimum  parameter 
values  in  the  controller  equations  were  obtained  by  relaxation  methods. 

The  detailed  results  will  be  published  in  W.  Morrison's  Ph.D.  disserta¬ 
tion. 
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4.2  FLUID-ELECTROLYTE  AND  RENAL  SYSTEMS 


4.2.1  Simulation  of  the  Combined  Artificial  Kidney-PatierU 
System 

GMO1724-02,  National  Institutes  of  Health 
M.  B.  Wolf 

This  research  proje  t  is  directed  towards  the  the  retie  ai 
prediction  of  optimum  parame  ers  for  hemodialysis  in  individual  "uremic" 
patients. 

A  mathematical  model  has  been  developed*  to  simulate  me 
combined  artificial  kidney-patient  chemical  dynamics,  and  has  been 
tested  against  actual  patient  data  for  two  of  the  prime  chemical  sub 
stances  of  interest.  Predictions  of  the  dynamics  of  excha..rje  of  these 
substances  as  a  function  of  patient  and  treatment  parameters  were  made, 
and  these  should  provide  the  clinician  with  a  rational  basis  for  choosing 
the  duration  of  the  dialysis  treatment. 

Further  efforts  are  underway  to  theoretically  investigate  the 
dynamics  of  other  substances  such  as  potas  mm  and  water  and  to  predict 
optimum  treatment  time  with  respect  to  them. 

Reference 

1.  Wolf,  M.  B.  ,  Watson,  P.  D.  ,  and  Barbour,  B.  H.  ,  "A  Theoretical 

Evaluation  of  the  Patient-Artificial  Kidney  System  Using  a  Kiil  Western 
Dialyzer,  (to  be  published  in  M- ‘■hematic?  1  Biosciences). 


-182- 


4.3  NEURAL  SYSTEMS 


4.3.1  Studies  of  Neuronal  Interaction 

NB  08207-01,  National  Institutes  of  Health 

G.  P.  Moore 

The  Neurophysiology  Research  Laboratory  of  the  Biomedical 
Engineering  program  is  pursuing  a  program  of  investigation  into  the  signal 
processing  behavior  of  small  neuron  networks.  In  the  past  we  have 
succeeded  in  developing,  in  conjunction  with  others,  the  mathematical 
techniques  for  describing  neuron  pulse  signals  and  the  computer  programs 
for  data  processing  and  simulation^  \  These  fundamental  studies  have 
enabled  us  to  characterize  nerve  cell  activity  and  to  detect  and  character¬ 
ize  interactions  between  cells. 

In  the  past  period,  therefore,  we  have  concentrated  on  (1) 
applying  the  basic  theory  to  real  neuron  on-line  experiments;  {2)  adapting 
large-scale  simulation  programs  (developed  elsewhere)  to  the  IBM  360/44 
computer  in  the  School  of  Engineering, 

An  extensive  literature  survey  in  the  use  of  models  in  the 
study  of  the  generation  and  processing  of  signals  by  nerve  cells  was 
completed  with  the  publication  of  USCEE  Report  290,  entitled  "Models  of 
the  Generation  and  Processing  of  Signals  by  Nerve  Cells:  A  C ategor ically 
Indexed  Abridged  Bibliography". 
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I ,  General  Background 


In  the  past  twelve  months  a  laboratory  capability  has  been 
developed  which  is  suitable  for  neurological  investigation  upon  both  verte¬ 
brate  and  invertebrate  specimens.  To  accomplish  this  has  required  the 
acquisition  and  integration  of  equipment  as  well  as  the  utilization  Oi 
available  equipment.  The  instrumentation  includes  microelectrodes, 
amplifiers  and  multi-channel  tape  recording  system,  CRT  display  devices, 
and  more  significantly,  a  PDP-8  digital  computer  with  appropriate  peri¬ 
pheral  devices  to  permit  both  on-  and  off-line  evaluation  of  experimental 
data.  A  library  of  programs  for  the  computer  is  being  developed  on  a 
continuing  basis.  An  extensive  capability  already  exis ts,  and.  with  the 
recent  addition  of  a  large  capacity  disc  storage,  these  programs  can  be 
kept  readily  available  for  call-up  on  an  on-line  basis.  Hence,  for  example, 
various  statistical  calculations  can  be  performed  on  experimental  data 
while  they  are  being  gathered,  the  results  displayed  for  evaluation  during 
the  experiment,  and  then  discarded  or  saved  on  disc  for  future  utilization 
as  appropriate. 

II.  Neural  -  Net  Modeling  Program 
Obj  ectives 

The  objective  of  this  study  is  to  formulate,  for  simulation  on 
a  general  purpose  digital  computer,  a  nerve-net  model.  This  model  is 
to  be  used  both  as  an  analytical  tool  in  research  concerning  the  functional 
relationships  of  neurons,  and  as  a  pedagogical  tool  for  use  by  students 
studying  functional  relationships  in  nerve  cells. 
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The  three  principal  areas  covered  in  the  study  are: 

1)  The  development  of  a  nerve-net  model  consisting  of  an  explicit 
statement  of  the  network  configuration,  and  of  the  nerve  cell 
dynamics . 

2)  Formulation  of  data  processing  package,  for  use  with  both 
the  spike-train  data  generated  by  the  simulation,  and 
experimentally  derived  data. 

3)  Development  of  a  digital  computer  program  which  permits 
the  model  to  be  used  as  a  general  simulation  tool  on  the 
USC  3^0/44  IBM  Computer. 

Results  and  Conclusions 

An  explicit  statement  of  the  nerve  model  is  formulated, 
including  a  detailed  description  of  the  network  configuration  and  of  the 
nerve  cell  dynamics.  This  formulation  is  presented  in  a  form  suitable 
for  use  in  further  work  in  the  areas  of  model  and  parameter  identification. 

Summary  of  the  Model 

The  neural  net  is  described  by 

s  -  fc\r  r  [r]1 

where:  S  is  an  n  vector  called  the  "stimulation  vector” 

i  c]  is  an  (nyn)  matrix  called  the  ''configuration  matrix” 

[  7  j  is  an  (nxn)  diagonal  matrix  called  the  ’extremal 
configuration  matrix” 
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t  is  an  m  vector  describing  the  input  to  the  network  from 
sources  exterior  to  the  network  and  is  called  the  "extremal 
cell  firing  vector". 

F  is  an  n  vector  describing  the  set  of  tiring  times  for  each 
cell  in  the  network  and  is  called  the  "internal  cell  firing 
vector". 

The  dynamics  of  the  n  ceils  are  represented  by  a  set  of  2n 
first  order  differential  equations 

x  =  [A] x  +  [BIS  +  [D] 

-1 

where:  x  =  ....  is  a  2n  state  vector.  Xj  representing 
x-> 

the  membrane  threshold  potential,  an  1  x^  representing 
the  trans-men  orane  potential. 

The  cell  firing  vector  T  is  generated  by  solving 
x  for  the  following  condition: 

T.  =  [  t/  |  x  -  x_  .  =  0]  v.,  2,...n. 

J  J 

The  nerve  model  is  shown  in  block  diagram  form: 
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Network  Configuration 


Cell  Dynamics 


S 


r 

x  =  Ax+  BS+  D 

r 

r  =  [i] 

L 

when 

lx  -  x- 1  =  0 

.•*-  — 

1  1  2  1 

Nerve  Net  Model 
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4.4  NEUROMUSCULAR  SYSTEMS 


4.4.1  Mathematical  Models  of  Force  Generation  in  Skeletal 

Muscle 

NB-06196-02,  National  Institutes  of  Health 

NGR-05-018-022,  National  Aeronautics  and  Space 

Administration 

J.  C.  Coggsho.ll  and  G.  A.  Bekey 

This  research  project  is  concerned  with  the  synthesis  of 
mathematical  models  of  the  force  generation  properties  of  skeletal 
muscle.  Experimental  data  have  been  obtained  by  recording  the  electro¬ 
myographic  activity  and  force  produced  by  elbow  flexor  and  extensor 
muscles  (specifically  the  biceps,  triceps,  and  brachioradialis  muscles) 
of  normal  human  subjects. 

Any  particular  muscle  may  be  viewed  as  composed  of  a  set 
of  motor  units.  The  motor  unit  (a  single  motoneuron  fiber  and  the  muscle 
fibers  it  innervates)  is  the  basic  functional  unit  of  muscle.  The  starting 
point  for  this  study  was  the  formulation  of  a  mathematical  model  for  the 
electrical  and  force  generating  activity  of  a  single  motor  unit  on  the  basis 
of  available  physiological  data.  The  motor  unit  model  was  used  to  postu¬ 
late  a  rn^del  for  the  aggregate  features  of  a  group  of  motor  units:  namely, 
the  myoelectric  activity  and  force  produced  by  a  whole  muscle.  In  addi¬ 
tion,  transfer  relationships  between  the  EMG  and  force  produced  by  a 
muscle  were  determined  using  on  line  parameter  tracking  with  an  analog 


computer. 


The  validity  of  both  the  postulated  whole  n.  ’icle  model  and 

the  EMG  to  force  transfer  characteristics  was  tested  by  cor  paring  these 

models  with  each  other  and  with  published  data  on  whole  muscles.  The 

resulting  comparisons  are  shown  in  the  amplitude  vs  frequency  plots  of 

Figure  1.  Curve  A  is  a  plot  of  the  amplitude  ratio  obtained  from  human 

tricepts  muscle  models  during  the  study*.  Curve  B  shows  the  results  of 

an  aggregate  model  obtained  by  summing  the  properties  of  single  motor 

2 

units  of  the  cat  soleus  muscle  .  Curve  C  illustrates  the  changes  of 

amplitude  obtained  by  Partridge  when  whole  cat  soleus  muscle  was 

3 

stimulated  by  sinusoidally  varying  nerve  impulse  frequencies  .  In  spite 
of  the  variety  of  sources,  there  is  general  agreement  among  these  curves 
over  a  wide  range  of  frequencies. 
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4.5  BIOMATHEMATICS 


4.  5.  1  Application  of  Control  Theory  to  Administration  of 

Digitalis  and  Antibiotics 

GM-16197-01,  Department  of  Health,  Education  and  Welfare 

J.  Buell,  R.  Kalaba,  and  R.  Bellman 

Buell  and  Kalaba,  working  closely  with  Jelliffe  of  the  LACGH, 
have  developed  mathematical  models  for  the  administration  of  such  drugs 
as  digitalis  and  antibiotics  using  control  theory  and  dynamic  Programming 
techniques.  These  preliminary  models  have  proved  surprisingly  effective. 
Programs  based  on  them  have  been  written,  and  now  are  actually  being  used 
in  the  hospital.  This  work  opens  the  door  to  the  development  of  more 
complex  and  significant  models  in  the  entire  area  of  drug  administration 
and  pharmacokinetics.  We  have  high  hopes  that  this  work  wili  significantly 
improve  patient  treatment.  Furthermore,  a  number  of  new  and  interesting 
types  of  control  processes  arise  from  this  study. 

An  investigation  has  been  made  by  Kaplan  of  certain  matrix 
theoretic  problems  connected  with  linear  compartmental  models  of  drug 
distribution  in  the  body.  For  several  categories  of  n-compartment 
models,  theorems  have  been  developed  showing  the  character  of  the 
eigenvalues  and  eigenvectors,  and  providing  upper  and  lower  bounds  on 
the  eigenvalues.  Further  theorems  have  been  developed  showing  what 
properties  of  the  system  are  deducible  in  principle  from  the  customary 
type  of  radioisotope  experiment.  Algorithms  for  making  these  deductions 


have  been  provided  and  incorporated  into  a  computer  program. 

Numerical  examples  have  been  worked  out  for  the  cases  of  digoxin  and 
albumin. 

Reference 

1.  J.  Buell  and  R.  Kalaba,  "Modern  Control  Theory  and  Optimal  Drug 
Regimens,  I:  The  Plateau  Effect.",  to  appear. 

4.5.2  Identification  of  Systems 

GI.1-16197-01,  Department  of  Health,  Education  and  Welfare 
J.  Bue1!,  R.  Kalaba  and  R.  Bellman 

0»'**r  the  last  year,  Buell  and  Kalaba,  working  with  a  well- 
known  physical-chemical  process ,  the  Michaelis-Menten  model ,  have 
applied  the  theory  of  quasilinearization  to  the  identification  of  rate 
constants  in  this  process.  The  successful  application^  of  methods  of 
this  type  will  have  a  tremendous  influence  upon  experimental  techniques 
in  chemistry,  biochem.  *  •  J  j  ^.neering.  Buell  and  Kalaba  have 
prepared  a  quasilinearization  program  which  will  make  these  methods 
readily  accessible.  A  great  deal  of  further  work  in  this  area  is  being 
done  with  the  aim  of  introducing  modern  sophisticated  techniques  into 
experimental  work. 

Reference 

1.  J.  Buell  and  R.  Kalaba,  "Quasilinearization  and  the  Filtering  of 

Nonlinear  Models  of  Drug  Metabolism  to  Experimental  Kinetic  Data1 
USC EE  Report  312. 
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